Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The Epoch of Reionization (EoR) depends on the complex astrophysics governing the birth and evolution of the first galaxies and structures in the intergalactic medium. EoR models rely on cosmic microwave background (CMB) observations, and in particular the large-scale E-mode polarization power spectra (EE PS), to help constrain their highly uncertain parameters. However, rather than directly forward-modelling the EE PS, most EoR models are constrained using a summary statistic -- the Thompson scattering optical depth, $\tau_e$. Compressing CMB observations to $\tau_e$ requires adopting a basis set for the EoR history. The common choice is the unphysical, redshift-symmetric hyperbolic tangent (Tanh) function, which differs in shape from physical EoR models based on hierarchical structure formation. Combining public EoR and CMB codes, 21cmFAST and CLASS, here we quantify how inference using the $\tau_e$ summary statistic impacts the resulting constraints on galaxy properties and EoR histories. Using the last Planck 2018 data release, we show that the marginalized constraints on the EoR history are more sensitive to the choice of the basis set (Tanh vs physical model) than to the CMB likelihood statistic ($\tau_e$ vs PS). For example, EoR histories implied by the growth of structure show a small tail of partial reionization extending to higher redshifts: the lower limit on the volume-averaged neutral hydrogen fraction at $z=10$ changes from $\bar{x}_{\rm HI} \ge 0.925 (1 \sigma)$ using Tanh to $\bar{x}_{\rm HI} \ge 0.849 (1 \sigma)$ using the physical model. However, biases in inference using $\tau_e$ are negligible for the Planck 2018 data.

Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2\sigma$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The very first galaxies that started the cosmic dawn likely resided in so-called "minihaloes", with masses of $\sim10^5$-$10^8\mathrm{M}_\odot$, accreting their gas from the intergalactic medium through H$_2$ cooling. Such molecularly cooled galaxies (MCGs) mostly formed in pristine environments, hosted massive, metal-free stars, and were eventually sterilized by the build-up of a disassociating (Lyman-Werner; LW) background. Therefore, their properties might be very different from the galaxies we see in the later Universe. Although MCGs are probably too faint to be observed directly, we could nevertheless infer their properties from the imprint they leave in the cosmic 21-cm signal. Here we quantify this imprint by extending the public simulation code 21cmFAST to allow for a distinct population of MCGs. We allow MCGs to have different properties from other galaxies, including unique scaling relations for their stellar-to-halo mass ratios, ionizing escape fractions, and spectral energy distributions. We track inhomogeneous recombinations, disassociative LW feedback, and photoheating from reionization. After demonstrating how MCGs can shape the 21-cm signal, we explore to what extent current observations can already place constraints on their properties. The cosmic microwave background optical depth from Planck sets an upper limit on the product of the ionizing escape fraction and the stellar mass in MCGs. When including also the timing of the putative EDGES absorption signal, we find an additional strong degeneracy between the stellar mass and the X-ray luminosity of MCGs. If proven to be of cosmic origin, the timing of the EDGES signal would have been set by MCGs.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$\alpha$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$\alpha$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$\alpha$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$\alpha$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.

Cross-correlating 21cm with known cosmic signals will be invaluable proof of the cosmic origin of the first 21cm detections. As some of the widest fields available, comprising thousands of sources with reasonably known redshifts, narrow-band Lyman alpha emitter (LAE) surveys are an obvious choice for such cross-correlation. Here we revisit the 21cm -- LAE cross-correlation, relaxing the common assumption of reionization occurring in a pre-heated intergalactic medium (IGM). Using specifications from the Square Kilometre Array and the Subary Hyper Supreme-Cam, we present new forecasts of the 21cm -- LAE cross-correlation function at $z\sim7$. We sample a broad parameter space of the mean IGM neutral fraction and spin temperature, ($\bar{x}_{\rm{HI}}$, $\bar{T}_{\rm S}$). The sign of the cross-correlation roughly follows the sign of the 21cm signal: ionized regions which surround LAEs correspond to relative hot spots in the 21cm signal when the neutral IGM is colder than the CMB, and relative cold spots when the neutral IGM is hotter than the CMB. The amplitude of the cross-correlation function generally increases with increasing $\bar{x}_{\rm{HI}}$, following the increasing bias of the cosmic HII regions. As is the case for 21cm, the strongest cross signal occurs when the IGM is colder than the CMB, providing a large contrast between the neutral regions and the ionized regions which host LAEs. We also vary the topology of reionization and the epoch of X-ray heating. The cross-correlation during the first half of reionization is sensitive to these topologies, and could thus be used to constrain them.

By mapping the large-scale thermal and ionization stucture of the intergalactic medium (IGM), cosmic 21-cm interferometry is set to revolutionize our understanding of the Epoch of Reionization (EoR) and the Cosmic Dawn (CD). However, interpreting upcoming observations relies on our ability to accurately model early galaxies and their cosmic radiation fields which drive the signal. The first galaxies are likely very rare and biased, with their abundances modulated by long wavelength modes of the density field. Moreover, the radiation fields from these galaxies interact with the IGM over a large range of scales. Therefore, the limited volume of 21-cm simulations can underestimate the amount of structure in the signal. In this work we quantify the minimum box size for simulating the power spectrum (PS) of the cosmic 21-cm signal. We perform multiple realizations of the initial conditions, for a range of box sizes. We quantify convergence with respect to a mock observation of box length 1.1 Gpc, with thermal noise computed for a 1000h observation with SKA1-low. We find that simulations of box lengths 200-300 Mpc typically do not show a bias in the PS during the EoR; however, they do tend to underestimate the large-scale power during the earlier epoch of CD by ~7-9 % on average. The absolute error on the PS from small-box realizations goes up to ~7 mK^2, or analogously up to ~4 sigma of the total noise. As expected, this bias and scatter between different realizations decreases with increasing box size. We conclude that box lengths of 250 Mpc are needed to converge at the level of <~1 sigma of the total noise.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

The intergalactic medium is expected to clump on scales down to $10^4-10^8$ M$_{\odot}$ before the onset of reionization. The impact of these small-scale structures on reionization is poorly understood despite the modern understanding that gas clumpiness limits the growth of H II regions. We use a suite of radiation-hydrodynamics simulations that capture the $\sim 10^4$ M$_{\odot}$ Jeans mass of unheated gas to study density fluctuations during reionization. Our simulations track the complex ionization and hydrodynamical response of gas in the wake of ionization fronts. The clumping factor of ionized gas (proportional to the recombination rate) rises to a peak value of $5-20$ approximately $\Delta t = 10$ Myr after ionization front passage, depending on the incident intensity, redshift, and degree to which the gas had been pre-heated by the first X-ray sources. The clumping factor reaches its relaxed value of $\approx 3$ by $\Delta t = 300$ Myr. The mean free path of Lyman-limit photons evolves in unison, being up to several times shorter in un-relaxed, recently reionized regions compared to those that were reionized much earlier. Assessing the impact of this response on the global reionizaton process, we find that un-relaxed gaseous structures boost the total number of recombinations by $\approx 50$ % and lead to spatial fluctuations in the mean free path that persist appreciably for several hundred million years after the completion of reionization.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

Galaxy formation during the first billion years of our Universe remains a challenging problem at the forefront of astrophysical cosmology. Although these $z \geq 6$ galaxies are likely responsible for the last major phase change of our Universe, the epoch of reionization (EoR), detailed studies are possible only for relatively rare, bright objects. Characterizing the fainter galaxies which are more representative of the population as a whole is currently done mainly through their non-ionizing UV luminosity function (LF). Observing the faint end of the UV LFs is nevertheless challenging, and current estimates can differ by orders of magnitude. Here we propose a methodology to combine disparate high-$z$ UV LF data sets in a Bayesian framework: Bayesian Data Averaging (BDA). Using a flexible, physically-motivated galaxy model, we compute the relative evidence of various $z=6$ UV LFs within the magnitude range $-20 \leq M_{\rm UV} \leq -15$ which is common to the data sets. Our model, based primarily on power-law scalings of the halo mass function, naturally penalizes systematically jagged data points as well as mis-estimated errors. We then use the relative evidence to weigh the posteriors obtained from disparate LF observations during the EoR, $6 \leq z \leq 10$. The resulting LFs suggest that the star formation rate density (SFRD) integrated down to a UV magnitude of -17 represent $60.9^{+11.3}_{-9.6}\%$ / $28.2^{+9.3}_{-10.1}\%$ / $5.7^{+4.5}_{-4.7}\%$ of the total SFRD at redshifts 6 / 10 / 15. The BDA framework we introduce enables galaxy models to leverage multiple, analogous observational data sets.

Next generation observatories will enable us to study the first billion years of our Universe in unprecedented detail. Foremost among these are 21-cm interferometry with the HERA and the SKA, and high-$z$ galaxy observations with the James Webb Space Telescope (JWST). Taking a basic galaxy model, in which we allow the star formation rates and ionizing escape fractions to have a power-law dependence on halo mass with an exponential turnover below some threshold, we quantify how observations from these instruments can be used to constrain the astrophysics of high-$z$ galaxies. For this purpose, we generate mock JWST LFs, based on two different hydrodynamical cosmological simulations; these have intrinsic luminosity functions (LFs) which turn over at different scales and yet are fully consistent with present-day observations. We also generate mock 21-cm power spectrum observations, using 1000h observations with SKA1 and a moderate foreground model. Using only JWST data, we predict up to a factor of 2-3 improvement (compared with HST) in the fractional uncertainty of the star formation rate to halo mass relation and the scales at which the LFs peak (i.e. turnover). Most parameters regulating the UV galaxy properties can be constrained at the level of $\sim 10$% or better, if either (i) we are able to better characterize systematic lensing uncertainties than currently possible; or (ii) the intrinsic LFs peak at magnitudes brighter than $M_{\rm UV} \lesssim -13$. Otherwise, improvement over HST-based inference is modest. When combining with upcoming 21-cm observations, we are able to significantly mitigate degeneracies, and constrain all of our astrophysical parameters, even for our most pessimistic assumptions about upcoming JWST LFs. The 21-cm observations also result in an order of magnitude improvement in constraints on the EoR history.

Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how different thermal statistics of the intergalactic medium can be used to constrain reionisation. Temperature is highly correlated with redshift of reionisation for a period of time after the gas is heated. However as the gas cools, thermal memory of reionisation is lost, and a power-law temperature-density relation is formed, $T = T_0(1+\delta)^{1-\gamma}$ with $\gamma \approx 1.5$. Constraining our model against observations of electron optical depth and temperature at mean density, we find that reionisation likely finished at $z_{\rm{reion}} = 6.8 ^{+ 0.5} _{-0.8}$ with a soft spectral slope of $\alpha = 2.8 ^{+ 1.2} _{-1.0}$. By restricting spectral slope to the range $[0.5,2.5]$ motivated by population II synthesis models, reionisation timing is further constrained to $z_{\rm{reion}} = 6.9 ^{+ 0.4} _{-0.5}$. We find that, in the future, the degeneracies between reionisation timing and background spectrum can be broken using the scatter in temperatures and integrated thermal history.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

Interferometry of the cosmic 21-cm signal is set to revolutionise our understanding of the Epoch of Reionisation (EoR) and the Cosmic Dawn (CD). The culmination of ongoing efforts will be the upcoming Square Kilometre Array (SKA), which will provide tomography of the 21-cm signal from the first billion years of our Universe. Using a galaxy formation model informed by high-$z$ luminosity functions, here we forecast the accuracy with which the first phase of SKA-low (SKA1-low) can constrain the properties of the unseen galaxies driving the astrophysics of the EoR and CD. We consider three observing strategies: (i) deep (1000h on a single field); (ii) medium-deep (100hr on 10 independent fields); and (iii) shallow (10hr on 100 independent fields). Using the 21-cm power spectrum as a summary statistic, and conservatively only using the 21-cm signal above the foreground wedge, we predict that all three observing strategies should recover astrophysical parameters to a fractional precision of $\sim 0.1$ -- 10 per cent. The reionisation history is recovered to an uncertainty of $\Delta z \lesssim 0.1$ (1$\sigma$) for the bulk of its duration. The medium-deep strategy, balancing thermal noise against cosmic variance, results in the tightest constraints, slightly outperforming the deep strategy. The shallow observational strategy performs the worst, with up to a $\sim 10$ -- 60 per cent increase in the recovered uncertainty. We note, however, that non-Gaussian summary statistics, tomography, as well as unbiased foreground removal would likely favour the deep strategy.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.

The 21-cm power spectrum (PS) has been shown to be a powerful discriminant of reionization and cosmic dawn astrophysical parameters. However, the 21-cm tomographic signal is highly non-Gaussian. Therefore there is additional information which is wasted if only the PS is used for parameter recovery. Here we showcase astrophysical parameter recovery directly from 21-cm images, using deep learning with convolutional neural networks (CNN). Using a database of 2D images taken from 10,000 21-cm lightcones (each generated from different cosmological initial conditions), we show that a CNN is able to recover parameters describing the first galaxies: (i) Tvir , their minimum host halo virial temperatures (or masses) capable of hosting efficient star formation; (ii) {\zeta} , their typical ionizing efficiencies; (iii) LX/SFR , their typical soft-band X-ray luminosity to star formation rate; and (iv) E0 , the minimum X-ray energy capable of escaping the galaxy into the IGM. For most of their allowed ranges, log Tvir and log LX/SFR are recovered with < 1% uncertainty, while {\zeta} and E0 are recovered within 10% uncertainty. Our results are roughly comparable to the accuracy obtained from Monte Carlo Markov Chain sampling of the PS with 21CMMC for the two mock observations analyzed previously, although we caution that we do not yet include noise and foreground contaminants in this proof-of-concept study.

The properties of the first galaxies, expected to drive the Cosmic Dawn (CD) and the Epoch of Reionization (EoR), are encoded in the 3D structure of the cosmic 21-cm signal. Parameter inference from upcoming 21-cm observations promises to revolutionize our understanding of these unseen galaxies. However, prior inference was done using models with several simplifying assumptions. Here we introduce a flexible, physically-motivated parametrization for high-$z$ galaxy properties, implementing it in the public code 21cmFAST. In particular, we allow their star formation rates and ionizing escape fraction to scale with the masses of their host dark matter halos, and directly compute inhomogeneous, sub-grid recombinations in the intergalactic medium. Combining current Hubble observations of the rest-frame UV luminosity function (UV LFs) at high-$z$ with a mock 1000h 21-cm observation using the Hydrogen Epoch of Reionization Arrays (HERA), we constrain the parameters of our model using a Monte Carlo Markov Chain sampler of 3D simulations, 21CMMC. We show that the amplitude and scaling of the stellar mass with halo mass is strongly constrained by LF observations, while the remaining galaxy properties are constrained mainly by 21-cm observations. The two data sets compliment each other quite well, mitigating degeneracies intrinsic to each observation. All eight of our astrophysical parameters are able to be constrained at the level of $\sim 10\%$ or better. The updated versions of 21cmFAST and 21CMMC used in this work are publicly available.

We directly compare predictions of dwarf galaxy properties in a semi-analytic model (SAM) with those extracted from a high-resolution hydrodynamic simulation. We focus on galaxies with halo masses of 1e9<Mvir/Msol<1e11 at high redshift ($z\ge5$). We find that, with the modifications previously proposed in Qin et al. (2018), including to suppress the halo mass and baryon fraction, as well as to modulate gas cooling and star formation efficiencies, the SAM can reproduce the cosmic evolution of galaxy properties predicted by the hydrodynamic simulation. These include the galaxy stellar mass function, total baryonic mass, star-forming gas mass and star formation rate at $z\sim5-11$. However, this agreement is only possible by reducing the star formation threshold relative to that suggested by local observations. Otherwise, too much star-forming gas is trapped in quenched dwarf galaxies. We further find that dwarf galaxies rapidly build up their star-forming reservoirs in the early universe ($z>10$), with the relevant time-scale becoming significantly longer towards lower redshifts. This indicates efficient accretion in cold mode in these low-mass objects at high redshift. Note that the improved SAM, which has been calibrated against hydrodynamic simulations, can provide more accurate predictions of high-redshift dwarf galaxy properties that are essential for reionization study.

We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

The recent detection of ULASJ1342+0928, a bright QSO at $z=7.54$, provides a powerful probe of the ionisation state of the intervening intergalactic medium, potentially allowing us to set strong constraints on the epoch of reionisation (EoR). Here we quantify the presence of Ly$\alpha$ damping wing absorption from the EoR in the spectrum of ULASJ1342+0928. Our Bayesian framework simultaneously accounts for uncertainties on: (i) the intrinsic QSO emission (obtained from reconstructing the Ly$\alpha$ profile from a covariance matrix of emission lines) and (ii) the distribution of HII regions during reionisation (obtained from three different 1.6$^3$ Gpc$^3$ simulations spanning the range of plausible EoR morphologies). Our analysis is complementary to that in the discovery paper (Ba\~nados et al.) and the accompanying method paper (Davies et al.) as it focuses solely on the damping wing imprint redward of Ly$\alpha$ ($1218 < \lambda < 1230$\AA), and uses a different methodology for (i) and (ii). We recover weak evidence for damping wing absorption. Our intermediate EoR model yields a volume-weighted neutral hydrogen fraction at $z=7.5$ of $\bar{x}_{\rm HI} = 0.21\substack{+0.17 \\ -0.19}$ (68 per cent). The constraints depend weakly on the EoR morphology. Our limits are lower than those presented previously, though they are consistent at ~1-1.5$\sigma$. We attribute this difference to: (i) a lower amplitude intrinsic Ly$\alpha$ profile obtained from our reconstruction pipeline, driven by correlations with other high-ionisation lines in the spectrum which are relatively weak; and (ii) only considering transmission redward of Ly$\alpha$ when computing the likelihood, which reduces the available constraining power but makes the results less model-dependent. Our results are consistent with previous estimates of the EoR history, and support the picture of a moderately extended EoR.

We use the sky-average spectrum measured by EDGES High-Band ($90-190$ MHz) to constrain parameters of early galaxies independent of the absorption feature at $78$~MHz reported by Bowman et al. (2018). These parameters represent traditional models of cosmic dawn and the epoch of reionization produced with the 21cmFAST simulation code (Mesinger & Furlanetto 2007, Mesinger et al. 2011). The parameters considered are: (1) the UV ionizing efficiency ($\zeta$), (2) minimum halo virial temperature hosting efficient star-forming galaxies ($T^{\rm min}_{\rm vir}$), (3) integrated soft-band X-ray luminosity ($L_{\rm X\,<\,2\,keV}/{\rm SFR}$), and (4) minimum X-ray energy escaping the first galaxies ($E_{0}$), corresponding to a typical H${\rm \scriptstyle I}$ column density for attenuation through the interstellar medium. The High-Band spectrum disfavors high values of $T^{\rm min}_{\rm vir}$ and $\zeta$, which correspond to signals with late absorption troughs and sharp reionization transitions. It also disfavors intermediate values of $L_{\rm X\,<\,2\,keV}/{\rm SFR}$, which produce relatively deep and narrow troughs within the band. Specifically, we rule out $39.4<\log_{10}\left(L_{\rm X\,<\,2\,keV}/{\rm SFR}\right)<39.8$ ($95\%$ C.L.). We then combine the EDGES High-Band data with constraints on the electron scattering optical depth from Planck and the hydrogen neutral fraction from high-$z$ quasars. This produces a lower degeneracy between $\zeta$ and $T^{\rm min}_{\rm vir}$ than that reported in Greig & Mesinger (2017a) using the Planck and quasar constraints alone. Our main result in this combined analysis is the estimate $4.5$~$\leq \log_{10}\left(T^{\rm min}_{\rm vir}/\rm K\right)\leq$~$5.7$ ($95\%$ C.L.). We leave for future work the evaluation of $21$~cm models using simultaneously data from EDGES Low- and High-Band.

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

Recent detections of Lyman alpha (Ly$\alpha$) emission from $z>7.5$ galaxies were somewhat unexpected given a dearth of previous non-detections in this era when the intergalactic medium (IGM) is still highly neutral. But these detections were from UV bright galaxies, which preferentially live in overdensities which reionize early, and have significantly Doppler-shifted Ly$\alpha$ line profiles emerging from their interstellar media (ISM), making them less affected by the global IGM state. Using a combination of reionization simulations and empirical ISM models we show, as a result of these two effects, UV bright galaxies in overdensities have $>2\times$ higher transmission through the $z\sim7$ IGM than typical field galaxies, and this boosted transmission is enhanced as the neutral fraction increases. The boosted transmission is not sufficient to explain the observed high Ly$\alpha$ fraction of $M_\mathrm{UV} \lesssim -22$ galaxies (Stark et al. 2017), suggesting Ly$\alpha$ emitted by these galaxies must be stronger than expected due to enhanced production and/or selection effects. Despite the bias of UV bright galaxies to reside in overdensities we show Ly$\alpha$ observations of such galaxies can accurately measure the global neutral hydrogen fraction, particularly when Ly$\alpha$ from UV faint galaxies is extinguished, making them ideal candidates for spectroscopic follow-up into the cosmic Dark Ages.

Using a Bayesian framework, we quantify what current observations imply about the history of the epoch of reionisation (EoR). We use a popular, three-parameter EoR model, flexible enough to accommodate a wide range of physically-plausible reionisation histories. We study the impact of various EoR observations: (i) the optical depth to the CMB measured by Planck 2016; (ii) the dark fraction in the Lyman $\alpha$ and $\beta$ forests; (iii) the redshift evolution of galactic Ly$\alpha$ emission (so-called "Ly$\alpha$ fraction"); (iv) the clustering of Ly$\alpha$ emitters; (v) the IGM damping wing imprint in the spectrum of QSO ULASJ1120+0641; (vi) and the patchy kinetic Sunyaev-Zel'dovich signal. Combined, (i) and (ii) already place interesting constraints on the reionisation history, with the epochs corresponding to an average neutral fraction of (75, 50, 25) per cent, constrained at 1$\sigma$ to $z= (9.21\substack{+1.22 -1.15}, 8.14\substack{+1.08 -1.00}, 7.26\substack{+1.13 -0.96})$. Folding-in more model-dependent EoR observations [(iii--vi)], strengthens these constraints by tens of per cent, at the cost of a decrease in the likelihood of the best-fit model, driven mostly by (iii). The tightest constraints come from (v). Unfortunately, no current observational set is sufficient to break degeneracies and constrain the astrophysical EoR parameters. However, model-dependent priors on the EoR parameters themselves can be used to set tight limits by excluding regions of parameter space with strong degeneracies. Motivated by recent observations of $z\sim7$ faint, lensed galaxies, we show how a conservative upper limit on the virial temperature of haloes which host reionising galaxies can constrain the escape fraction of ionising photons to $f_{\rm esc} = 0.14\substack{+0.26 -0.09}$

The cosmic 21 cm signal is set to revolutionise our understanding of the early Universe, allowing us to probe the 3D temperature and ionisation structure of the intergalactic medium (IGM). It will open a window onto the unseen first galaxies, showing us how their UV and X-ray photons drove the cosmic milestones of the epoch of reionisation (EoR) and epoch of heating (EoH). To facilitate parameter inference from the 21 cm signal, we previously developed 21CMMC: a Monte Carlo Markov Chain sampler of 3D EoR simulations. Here we extend 21CMMC to include simultaneous modelling of the EoH, resulting in a complete Bayesian inference framework for the astrophysics dominating the observable epochs of the cosmic 21 cm signal. We demonstrate that second generation interferometers, the Hydrogen Epoch of Reionisation Array (HERA) and Square Kilometre Array (SKA) will be able to constrain ionising and X-ray source properties of the first galaxies with a fractional precision of order $\sim1$-10 per cent (1$\sigma$). The ionisation history of the Universe can be constrained to within a few percent. Using our extended framework, we quantify the bias in EoR parameter recovery incurred by the common simplification of a saturated spin temperature in the IGM. Depending on the extent of overlap between the EoR and EoH, the recovered astrophysical parameters can be biased by $\sim3-10\sigma$.

We quantify the presence of Ly\alpha\ damping wing absorption from a partially-neutral intergalactic medium (IGM) in the spectrum of the $z=7.08$ QSO, ULASJ1120+0641. Using a Bayesian framework, we simultaneously account for uncertainties in: (i) the intrinsic QSO emission spectrum; and (ii) the distribution of cosmic HI patches during the epoch of reionisation (EoR). For (i) we use a new intrinsic Ly\alpha\ emission line reconstruction method (Greig et al.), sampling a covariance matrix of emission line properties built from a large database of moderate-$z$ QSOs. For (ii), we use the Evolution of 21-cm Structure (EOS; Mesinger et al.) simulations, which span a range of physically-motivated EoR models. We find strong evidence for the presence of damping wing absorption redward of Ly\alpha\ (where there is no contamination from the Ly\alpha\ forest). Our analysis implies that the EoR is not yet complete by $z=7.1$, with the volume-weighted IGM neutral fraction constrained to $\bar{x}_{\rm H\,{\scriptsize I}} = 0.40\substack{+0.21 -0.19}$ at $1\sigma$ ($\bar{x}_{\rm H\,{\scriptsize I}} = 0.40\substack{+0.41 -0.32}$ at $2\sigma$). This result is insensitive to the EoR morphology. Our detection of significant neutral HI in the IGM at $z=7.1$ is consistent with the latest Planck 2016 measurements of the CMB Thompson scattering optical depth (Planck Collaboration XLVII).

We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z~6-10) from detections and non-detections of Lyman Alpha (Ly$\alpha$) emission from Lyman break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Ly$\alpha$. We assert that the Ly$\alpha$ line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Ly$\alpha$ lines from host galaxies' systemic redshifts, which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. (2014) and infer a global neutral fraction at z~7 of $\overline{x}_\mathrm{HI} = 0.59_{-0.15}^{+0.11}$, consistent with other robust probes of the EoR and confirming reionization is on-going ~700 Myr after the Big Bang. We show that using the full distribution of Ly$\alpha$ equivalent width detections and upper limits from LBGs places tighter constraints on the evolving IGM than the standard Ly$\alpha$ emitter fraction, and that larger samples are within reach of deep spectroscopic surveys of gravitationally lensed fields and JWST NIRSpec.

We highlight recent progress in the sophistication and diversification of cosmic dawn and reionization simulations. The application of these modeling tools to current observations has allowed us narrow down the timing of reionization, which we now know to within dz ~ 1 for the bulk of reionization. The strongest constraints come from the optical depth to the CMB measured with the {\it Planck} satellite and the first detection of ongoing reionization from the spectra of the z=7.1 QSOs ULASJ1120+0641. However, we still know virtually nothing about the astrophysical sources during the first billion years. The revolution in our understanding will be led by upcoming interferometric observations of the cosmic 21-cm signal. The properties of the sources and sinks of UV and X-ray photons are encoded in the 3D patterns of the signal. The development of Bayesian parameter recovery techniques, which tap into the wealth of the 21-cm signal, will soon usher in an era of precision astrophysical cosmology.

We extend 21CMMC, a Monte Carlo Markov Chain sampler of 3D reionisation simulations, to perform parameter estimation directly on 3D light-cones of the cosmic 21cm signal. This brings theoretical analysis closer to the tomographic 21-cm observations achievable with next generation interferometers like HERA and the SKA. Parameter recovery can therefore account for modes which evolve with redshift/frequency. Additionally, simulated data can be more easily corrupted to resemble real data. Using the light-cone version of 21CMMC, we quantify the biases in the recovered astrophysical parameters if we use the 21cm power spectrum from the co-evolution approximation to fit a 3D light-cone mock observation. While ignoring the light-cone effect under most assumptions will not significantly bias the recovered astrophysical parameters, it can lead to an underestimation of the associated uncertainty. However significant biases ($\sim$few -- 10 $\sigma$) can occur if the 21cm signal evolves rapidly (i.e. the epochs of reionisation and heating overlap significantly) and: (i) foreground removal is very efficient, allowing large physical scales ($k\lesssim0.1$~Mpc$^{-1}$) to be used in the analysis or (ii) theoretical modelling is accurate to within $\sim10$ per cent in the power spectrum amplitude.

Current and upcoming radio interferometric experiments are aiming to make a statistical characterization of the high-redshift 21cm fluctuation signal spanning the hydrogen reionization and X-ray heating epochs of the universe. However, connecting 21cm statistics to underlying physical parameters is complicated by the theoretical challenge of modeling the relevant physics at computational speeds quick enough to enable exploration of the high dimensional and weakly constrained parameter space. In this work, we use machine learning algorithms to build a fast emulator that mimics expensive simulations of the 21cm signal across a wide parameter space to high precision. We embed our emulator within a Markov-Chain Monte Carlo framework, enabling it to explore the posterior distribution over a large number of model parameters, including those that govern the Epoch of Reionization, the Epoch of X-ray Heating, and cosmology. As a worked example, we use our emulator to present an updated parameter constraint forecast for the Hydrogen Epoch of Reionization Array experiment, showing that its characterization of a fiducial 21cm power spectrum will considerably narrow the allowed parameter space of reionization and heating parameters, and could help strengthen Planck's constraints on $\sigma_8$. We provide both our generalized emulator code and its implementation specifically for 21cm parameter constraints as publicly available software.

Motivated by recent measurements of the number density of faint AGN at high redshift, we investigate the contribution of quasars to reionization by tracking the growth of central supermassive black holes in an update of the Meraxes semi-analytic model. The model is calibrated against the observed stellar mass function at $z\sim0.6-7$, the black hole mass function at $z\lesssim0.5$, the global ionizing emissivity at $z\sim2-5$ and the Thomson scattering optical depth. The model reproduces a Magorrian relation in agreement with observations at $z<0.5$ and predicts a decreasing black hole mass towards higher redshifts at fixed total stellar mass. With the implementation of an opening angle of 80 deg for quasar radiation, corresponding to an observable fraction of ${\sim}23.4$ per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at $z\sim0.6-6$. The stellar light from galaxies hosting faint AGN contributes a significant or dominant fraction of the UV flux. At high redshift, the model is consistent with the bright end quasar luminosity function and suggests that the recent faint $z\sim4$ AGN sample compiled by Giallongo et al. (2015) includes a significant fraction of stellar light. Direct application of this luminosity function to the calculation of AGN ionizing emissivity consequently overestimates the number of ionizing photons produced by quasars by a factor of 3 at $z\sim6$. We conclude that quasars are unlikely to make a significant contribution to reionization.

We investigate the clustering properties of Lyman-break galaxies (LBGs) at $z\sim6$ - $8$. Using the semi-analytical model {\scshape Meraxes} constructed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation (DRAGONS) project, we predict the angular correlation function (ACF) of LBGs at $z\sim6$ - $8$. Overall, we find that the predicted ACFs are in good agreement with recent measurements at $z\sim 6$ and $z\sim 7.2$ from observations consisting of the Hubble eXtreme Deep Field (XDF), the Hubble Ultra-Deep Field (HUDF) and Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. We confirm the dependence of clustering on luminosity, with more massive dark matter haloes hosting brighter galaxies, remains valid at high redshift. The predicted galaxy bias at fixed luminosity is found to increase with redshift, in agreement with observations. We find that LBGs of magnitude $M_{{\rm AB(1600)}} < -19.4$ at $6\lesssim z \lesssim 8$ reside in dark matter haloes of mean mass $\sim 10^{11.0}$- $10^{11.5} M_{\rm \odot}$, and this dark matter halo mass does not evolve significantly during reionisation.

Massive quiescent galaxies (MQGs) are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a MQG, ZF-COSMOS-20115 at $z\sim4$, only 1.5 Gyr after the big bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a MQG at $z=3.717$ with an estimated stellar mass of ${\sim}10^{11}\mathrm{M}_\odot$, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ${\sim}10^{13}\mathrm{M}_\odot$, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. ZF-COSMOS-20115 likely remained quiescent at $z<3.7$. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at $z\sim5$, indicating that ZF-COSMOS-20115 and other MQGs may be the descendants of high-redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early.

Direct detection of regions of ionized hydrogen (HII) has been suggested as a promising probe of cosmic reionization. Observing the redshifted 21-cm signal of hydrogen from the epoch of reionization (EoR) is a key scientific driver behind new-generation, low-frequency radio interferometers. We investigate the feasibility of combining low-frequency observations with the Square Kilometre Array and near infra-red survey data of the Wide-Field Infrared Survey Telescope to detect cosmic reionization by imaging HII bubbles surrounding massive galaxies during the cosmic dawn. While individual bubbles will be too small to be detected, we find that by stacking redshifted 21-cm spectra centred on known galaxies, it will be possible to directly detect the EoR at $z \sim 9-12$, and to place qualitative constraints on the evolution of the spin temperature of the intergalactic medium (IGM) at $z \geq 9$. In particular, given a detection of ionized bubbles using this technique, it is possible to determine if the IGM surrounding them is typically in absorption or emission. Determining the globally-averaged neutral fraction of the IGM using this method will prove more difficult due to degeneracy with the average size of HII regions.

Using a series of high-resolution hydrodynamical simulations we show that during the rapid growth of high-redshift (z > 5) galaxies, reserves of molecular gas are consumed over a time-scale of 300Myr, almost independent of feedback scheme. We find that there exists no such simple relation for the total gas fractions of these galaxies, with little correlation between gas fractions and specific star formation rates. The bottleneck or limiting factor in the growth of early galaxies is in converting infalling gas to cold star-forming gas. Thus, we find that the majority of high redshift dwarf galaxies are effectively in recession, with demand (of star formation) never rising to meet supply (of gas), irrespective of the baryonic feedback physics modelled. We conclude that the basic assumption of self-regulation in galaxies - that they can adjust total gas consumption within a Hubble time - does not apply for the dwarf galaxies thought to be responsible for providing most UV photons to reionize the high redshift Universe. We demonstrate how this rapid molecular time-scale improves agreement between semi-analytic model predictions of the early Universe and observed stellar mass functions.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.

We investigate how the hydrostatic suppression of baryonic accretion affects the growth rate of dark matter halos during the Epoch of Reionization. By comparing halo properties in a simplistic hydrodynamic simulation in which gas only cools adiabatically, with its collisionless equivalent, we find that halo growth is slowed as hydrostatic forces prevent gas from collapsing. In our simulations, at the high redshifts relevant for reionization (between ${\sim}6$ and ${\sim}11$), halos that host dwarf galaxies ($\lesssim 10^{9} \mathrm{M_\odot}$) can be reduced by up to a factor of 2 in mass due to the hydrostatic pressure of baryons. Consequently, the inclusion of baryonic effects reduces the amplitude of the low mass tail of the halo mass function by factors of 2 to 4. In addition, we find that the fraction of baryons in dark matter halos hosting dwarf galaxies at high redshift never exceeds ${\sim}90\%$ of the cosmic baryon fraction. When implementing baryonic processes, including cooling, star formation, supernova feedback and reionization, the suppression effects become more significant with further reductions of ${\sim}30\%$ to 60\%. Although convergence tests suggest that the suppression may become weaker in higher resolution simulations, this suppressed growth will be important for semi-analytic models of galaxy formation, in which the halo mass inherited from an underlying N-body simulation directly determines galaxy properties. Based on the adiabatic simulation, we provide tables to account for these effects in N-body simulations, and present a modification of the halo mass function along with explanatory analytic calculations.

By heating the intergalactic medium (IGM) before reionization, X-rays are expected to play a prominent role in the early Universe. The cosmic 21-cm signal from this "Epoch of Heating" (EoH) could serve as a clean probe of high-energy processes inside the first galaxies. Here we improve on prior estimates of this signal by using high-resolution hydrodynamic simulations to calculate the X-ray absorption due to the interstellar medium (ISM) of the host galaxy. X-rays absorbed inside the host galaxy are unable to escape into the IGM and contribute to the EoH. We find that the X-ray opacity through these galaxies can be approximated by a metal-free ISM with a typical column density of log[N / cm^-2] = 21.4 +0.40-0.65. We compute the resulting 21-cm signal by combining these ISM opacities with public spectra of high-mass X-ray binaries (thought to be important X-ray sources in the early Universe). Our results support "standard scenarios" in which the X-ray heating of the IGM is inhomogeneous, and occurs before the bulk of reionization. The large-scale (k ~ 0.1/Mpc) 21-cm power reaches a peak of ~100 mK^2 at z = 10 - 15, with the redshift depending on the cosmic star formation history. This is in contrast to some recent work, motivated by the much larger X-ray absorption towards local HMXBs inside the Milky Way. Our main results can be reproduced by approximating the X-ray emission from HMXBs with a power-law spectrum with energy index alpha = 1, truncated at energies below 0.5 keV.

Using Hubble data, including new grism spectra, Oesch et al. recently identified GN-z11, an $M_\textrm{UV}$=-21.1 galaxy at $z$=11.1 (just 400Myr after the big bang). With an estimated stellar mass of $\sim$10$^9$M$_{\odot}$, this galaxy is surprisingly bright and massive, raising questions as to how such an extreme object could form so early in the Universe. Using \Meraxes{}, a semi-analytic galaxy-formation model developed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) programme, we investigate the potential formation mechanisms and eventual fate of GN-z11. The volume of our simulation is comparable to that of the discovery observations and possesses two analogue galaxies of similar luminosity to this remarkably bright system. Existing in the two most massive subhaloes at $z$=11.1 ($M_\textrm{vir}$=1.4$\times 10^{11}$M$_{\odot}$ and 6.7$\times 10^{10}$M$_{\odot}$), our model analogues show excellent agreement with all available observationally derived properties of GN-z11. Although they are relatively rare outliers from the full galaxy population at high-$z$, they are no longer the most massive or brightest systems by $z$=5. Furthermore, we find that both objects possess relatively smooth, but extremely rapid mass growth histories with consistently high star formation rates and UV luminosities at $z{>}11$, indicating that their brightness is not a transient, merger-driven feature. Our model results suggest that future wide-field surveys with the \textit{James Webb Space Telescope} may be able to detect the progenitors of GN-z11 analogues out to $z{\sim}$13--14, pushing the frontiers of galaxy-formation observations to the early phases of cosmic reionization and providing a valuable glimpse of the first galaxies to reionize the Universe on large scales.

We introduce an intrinsic Ly\alpha\ emission line profile reconstruction method for high-$z$ quasars (QSOs). This approach utilises a covariance matrix of emission line properties obtained from a large, moderate-$z$ ($2 \leq z \leq 2.5$), high signal to noise (S/N > 15) sample of BOSS QSOs. For each QSO, we complete a Monte Carlo Markov Chain fitting of the continuum and emission line properties and perform a visual quality assessment to construct a large database of robustly fit spectra. With this dataset, we construct a covariance matrix to describe the correlations between the high ionisation emission lines Ly\alpha, C IV, Si IV + O IV] and C III], and find it to be well approximated by an $N$-dimensional Gaussian distribution. This covariance matrix characterises the correlations between the line width, peak height and velocity offset from systemic while also allowing for the existence of broad and narrow line components for Ly\alpha\ and C IV. We illustrate how this covariance matrix allows us to statistically characterise the intrinsic Ly\alpha\ line solely from the observed spectrum redward of 1275\AA. This procedure can be used to reconstruct the intrinsic Ly\alpha\ line emission profile in cases where Ly\alpha\ may otherwise be obscured. Applying this reconstruction method to our sample of QSOs, we recovered the Ly\alpha\ line flux to within 15 per cent of the measured flux at 1205\AA\ (1220\AA) ~85 (90) per cent of the time.

We investigate high-redshift galaxy sizes using a semi-analytic model constructed for the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation project. Our fiducial model, including strong feedback from supernovae and photoionization background, accurately reproduces the evolution of the stellar mass function and UV luminosity function. Using this model, we study the size--luminosity relation of galaxies and find that the effective radius scales with UV luminosity as $R_\mathrm{e}\propto L^{0.25}$ at $z{\sim}5$--$9$. We show that recently discovered very luminous galaxies at $z{\sim}7$ (Bowler et al. 2016) and $z{\sim}11$ (Oesch et al. 2016) lie on our predicted size--luminosity relations. We find that a significant fraction of galaxies at $z>8$ will not be resolved by JWST, but GMT will have the ability to resolve all galaxies in haloes above the atomic cooling limit. We show that our fiducial model successfully reproduces the redshift evolution of average galaxy sizes at $z>5$. We also explore galaxy sizes in models without supernova feedback. The no-supernova feedback models produce galaxy sizes that are smaller than observations. We therefore confirm that supernova feedback plays an important role in determining the size--luminosity relation of galaxies and its redshift evolution during reionization.

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology.

We introduce Meraxes, a new, purpose-built semi-analytic galaxy formation model designed for studying galaxy growth during reionization. Meraxes is the first model of its type to include a temporally and spatially coupled treatment of reionization and is built upon a custom (100 Mpc)$^3$ N-body simulation with high temporal and mass resolution, allowing us to resolve the galaxy and star formation physics relevant to early galaxy formation. Our fiducial model with supernova feedback reproduces the observed optical depth to electron scattering and evolution of the galaxy stellar mass function between $z$=5 and 7, predicting that a broad range of halo masses contribute to reionization. Using a constant escape fraction and global recombination rate, our model is unable to simultaneously match the observed ionizing emissivity at $z{\lesssim}6$. However, the use of an evolving escape fraction of 0.05-0.1 at $z{\sim}6$, increasing towards higher redshift, is able to satisfy these three constraints. We also demonstrate that photoionization suppression of low mass galaxy formation during reionization has only a small effect on the ionization history of the inter-galactic medium. This lack of `self-regulation' arises due to the already efficient quenching of star formation by supernova feedback. It is only in models with gas supply-limited star formation that reionization feedback is effective at regulating galaxy growth. We similarly find that reionization has only a small effect on the stellar mass function, with no observationally detectable imprint at $M_{\rm *}{>}10^{7.5}\,{\rm M_{{\odot}}}$. However, patchy reionization has significant effects on individual galaxy masses, with variations of factors of 2-3 at $z$=5 that correlate with environment.

We use the Dark-ages, Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) framework to investigate the effect of galaxy-formation physics on the morphology and statistics of ionized hydrogen (HII) regions during the Epoch of Reioinization (EoR). DRAGONS self-consistently couples a semi-analytic galaxy-formation model with the inhomogeneous ionizing UV background, and can therefore be used to study the dependence of morphology and statistics of reionization on feedback phenomena of the ionizing source galaxy population. Changes in galaxy-formation physics modify the sizes of HII regions and the amplitude and shape of 21-cm power spectra. Of the galaxy physics investigated, we find that supernova feedback plays the most important role in reionization, with HII regions up to $\approx 20$ per cent smaller and a fractional difference in the amplitude of power spectra of up to $\approx 17$ per cent at fixed ionized fraction in the absence of this feedback. We compare our galaxy-formation-based reionization models with past calculations that assume constant stellar-to-halo mass ratios and find that with the correct choice of minimum halo mass, such models can mimic the predicted reionization morphology. Reionization morphology at fixed neutral fraction is therefore not uniquely determined by the details of galaxy formation, but is sensitive to the mass of the haloes hosting the bulk of the ionizing sources. Simple EoR parametrizations are therefore accurate predictors of reionization statistics. However, a complete understanding of reionization using future 21-cm observations will require interpretation with realistic galaxy-formation models, in combination with other observations.

We present upper limits on the 21 cm power spectrum at $z = 5.9$ calculated from the model-independent limit on the neutral fraction of the intergalactic medium of $x_{\rm H{\small I }} < 0.06 + 0.05\ (1\sigma)$ derived from dark pixel statistics of quasar absorption spectra. Using 21CMMC, a Markov chain Monte Carlo Epoch of Reionization analysis code, we explore the probability distribution of 21 cm power spectra consistent with this constraint on the neutral fraction. We present 99 per cent confidence upper limits of $\Delta^2(k) < 10$ to $20\ {\rm mK}^2$ over a range of $k$ from 0.5 to $2.0\ h{\rm Mpc}^{-1}$, with the exact limit dependent on the sampled $k$ mode. This limit can be used as a null test for 21 cm experiments: a detection of power at $z=5.9$ in excess of this value is highly suggestive of residual foreground contamination or other systematic errors affecting the analysis.

The birth of the first luminous sources and the ensuing epoch of reionization are best studied via the redshifted 21-cm emission line, the signature of the first two imprinting the last. In this work we present a fully-Bayesian method, \textsc{hibayes}, for extracting the faint, global (sky-averaged) 21-cm signal from the much brighter foreground emission. We show that a simplified (but plausible), Gaussian model of the 21-cm emission from the Cosmic Dawn epoch ($15 \lesssim z \lesssim 30$), parameterized by an amplitude $A_{\rm HI}$, a frequency peak $\nu_{\rm HI}$ and a width $\sigma_{\rm HI}$, can be extracted even in the presence of a structured foreground frequency spectrum (parameterized as a $7^{\rm th}$-order polynomial), provided sufficient signal-to-noise (400~hours of observation with a single dipole). We apply our method to an early, 19-minute long observation from the Large aperture Experiment to detect the Dark Ages, constraining the 21-cm signal amplitude and width to be $-890 < A_{\rm HI} < 0$ mK and $\sigma_{\rm HI} > 6.5$ MHz (corresponding to $\Delta z > 1.9$ at redshift $z \simeq 20$) respectively at the 95-per-cent confidence level in the range $13.2 < z < 27.4$ ($100 > \nu > 50$ MHz).

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window.

In this paper we present calculations of the UV luminosity function from the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) project, which combines N-body, semi-analytic and semi-numerical modelling designed to study galaxy formation during the Epoch of Reionization. Using galaxy formation physics including supernova feedback, the model naturally reproduces the UV LFs for high-redshift star-forming galaxies from $z{\sim}5$ through to $z{\sim}10$. We investigate the luminosity--star formation rate (SFR) relation, finding that variable SFR histories of galaxies result in a scatter around the median relation of $0.1$--$0.3$ dex depending on UV luminosity. We find close agreement between the model and observationally derived SFR functions. We use our calculated luminosities to investigate the luminosity function below current detection limits, and the ionizing photon budget for reionization. We predict that the slope of the UV LF remains steep below current detection limits and becomes flat at $M_\mathrm{UV}{\gtrsim}{-14}$. We find that $48$ ($17$) per cent of the total UV flux at $z{\sim}6$ ($10$) has been detected above an observational limit of $M_\mathrm{UV}{\sim}{-17}$, and that galaxies fainter than $M_\mathrm{UV}{\sim}{-17}$ are the main source of ionizing photons for reionization. We investigate the luminosity--stellar mass relation, and find a correlation for galaxies with $M_\mathrm{UV}{<}{-14}$ that has the form $M_*{\propto}10^{-0.47M_\mathrm{UV}}$, in good agreement with observations, but which flattens for fainter galaxies. We determine the luminosity--halo mass relation to be $M_\mathrm{vir}{\propto}10^{-0.35M_\mathrm{UV}}$, finding that galaxies with $M_\mathrm{UV}{=}{-20}$ reside in host dark matter haloes of $10^{11.0\pm 0.1}\mathrm{M_\odot}$ at $z{\sim}6$, and that this mass decreases towards high redshift.

We use the Fisher matrix formalism and semi-numerical simulations to derive quantitative predictions of the constraints that power spectrum measurements on next-generation interferometers, such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA), will place on the characteristics of the X-ray sources that heated the high redshift intergalactic medium. Incorporating observations between $z=5$ and $z=25$, we find that the proposed 331 element HERA and SKA phase 1 will be capable of placing $\lesssim 10\%$ constraints on the spectral properties of these first X-ray sources, even if one is unable to perform measurements within the foreground contaminated "wedge" or the FM band. When accounting for the enhancement in power spectrum amplitude from spin temperature fluctuations, we find that the observable signatures of reionization extend well beyond the peak in the power spectrum usually associated with it. We also find that lower redshift degeneracies between the signatures of heating and reionization physics lead to errors on reionization parameters that are significantly greater than previously predicted. Observations over the heating epoch are able to break these degeneracies and improve our constraints considerably. For these two reasons, 21\,cm observations during the heating epoch significantly enhance our understanding of reionization as well.

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.

We introduce the Evolution of 21-cm Structure (EOS) project: providing periodic, public releases of the latest cosmological 21-cm simulations. 21-cm interferometry is set to revolutionize studies of the Cosmic Dawn (CD) and epoch of reionization (EoR), eventually resulting in 3D maps of the first billion years of our Universe. Progress will depend on sophisticated data analysis pipelines, which are in turn tested on large-scale mock observations. Here we present the 2016 EOS data release, consisting of the largest (1.6 Gpc on side with a 1024^3 grid), public 21-cm simulations of the CD and EoR. We include calibrated, sub-grid prescriptions for inhomogeneous recombinations and photo-heating suppression of star formation in small mass galaxies. We present two simulation runs that approximately bracket the contribution from faint unseen galaxies. From these two extremes, we predict that the duration of reionization (defined as a change in the mean neutral fraction from 0.9 to 0.1) should be between 2.7 < Delta z < 5.7. The large-scale 21-cm power during the advanced EoR stages can be different by up to a factor of ~10, depending on the model. This difference has a comparable contribution from: (i) the typical bias of sources; and (ii) a more efficient negative feedback in models with an extended EoR driven by faint galaxies. We also make detectability forecasts. With a 1000h integration, HERA and SKA1-low should achieve a signal-to-noise of ~few-hundreds throughout the EoR/CD, while in the maximally optimistic scenario of perfect foreground cleaning, all instruments should make a statistical detection of the cosmic signal. We also caution that our ability to clean foregrounds determines the relative performance of narrow/deep vs. wide/shallow surveys expected with SKA1. Our 21-cm power spectra, simulation outputs and visualizations are publicly available.

Interferometry of the cosmic 21-cm signal is set to revolutionize our understanding of the Epoch of Reionization (EoR), eventually providing 3D maps of the early Universe. Initial detections however will be low signal-to-noise, limited by systematics. To confirm a putative 21-cm detection, and check the accuracy of 21-cm data analysis pipelines, it would be very useful to cross-correlate against a genuine cosmological signal. The most promising cosmological signals are wide-field maps of Lyman alpha emitting galaxies (LAEs), expected from the Subaru Hyper-Suprime Cam (HSC) Ultra-Deep field. Here we present estimates of the correlation between LAE maps at z~7 and the 21-cm signal observed by both the Low Frequency Array (LOFAR) and the planned Square Kilometer Array Phase 1 (SKA1). We adopt a systematic approach, varying both: (i) the prescription of assigning LAEs to host halos; and (ii) the large-scale structure of neutral and ionized regions (i.e. EoR morphology). We find that the LAE-21cm cross-correlation is insensitive to (i), thus making it a robust probe of the EoR. A 1000h observation with LOFAR would be sufficient to discriminate at >1 standard deviation a fully ionized Universe from one with a mean neutral fraction of xHI~0.50, using the LAE-21cm cross-correlation function on scales of R~3-10 Mpc. Unlike LOFAR, whose detection of the LAE-21cm cross-correlation is limited by noise, SKA1 is mostly limited by ignorance of the EoR morphology. However, the planned 100h wide-field SKA1-Low survey will be sufficient to discriminate an ionized Universe from one with xHI~0.25, even with maximally pessimistic assumptions.

We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically-averaged density profiles of individual haloes. After removing haloes that are out-of-equilibrium, we find a $z{>}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $\Delta c_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $\alpha$, depends on peak height, $\nu$, in a manner that is accurately described by $\alpha {=}0.0070\nu^2{+}0.1839$. The distribution of the spin parameter, $\lambda$, has a weak dependence on equilibrium state; $\lambda$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift.

We compute robust lower limits on the spin temperature, $T_{\rm S}$, of the $z=8.4$ intergalactic medium (IGM), implied by the upper limits on the 21-cm power spectrum recently measured by PAPER-64. Unlike previous studies which used a single epoch of reionization (EoR) model, our approach samples a large parameter space of EoR models: the dominant uncertainty when estimating constraints on $T_{\rm S}$. Allowing $T_{\rm S}$ to be a free parameter and marginalizing over EoR parameters in our Markov Chain Monte Carlo code 21CMMC, we infer $T_{\rm S}\ge3 {\rm K}$ (corresponding approximately to $1\sigma$) for a mean IGM neutral fraction of $\bar{x}_{\rm H{\scriptsize I}}\gtrsim0.1$. We further improve on these limits by folding-in additional EoR constraints based on: (i) the dark fraction in QSO spectra, which implies a strict upper limit of $\bar{x}_{\rm H{\scriptsize I}}[z=5.9]\leq 0.06+0.05 \,(1\sigma)$; and (ii) the electron scattering optical depth, $\tau_{\rm e}=0.066\pm0.016\,(1\sigma)$ measured by the Planck satellite. By restricting the allowed EoR models, these additional observations tighten the approximate $1\sigma$ lower limits on the spin temperature to $T_{\rm S} \ge 6$ K. Thus, even such preliminary 21-cm observations begin to rule out extreme scenarios such as `cold reionization', implying at least some prior heating of the IGM. The analysis framework developed here can be applied to upcoming 21-cm observations, thereby providing unique insights into the sources which heated and subsequently reionized the very early Universe.

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

With the first phase of the Square Kilometre Array (SKA1) entering into its final pre-construction phase, we investigate how best to maximise its scientific return. Specifically, we focus on the statistical measurement of the 21 cm power spectrum (PS) from the epoch of reionization (EoR) using the low frequency array, SKA1-low. To facilitate this investigation we use the recently developed MCMC based EoR analysis tool 21CMMC (Greig & Mesinger). In light of the recent 50 per cent cost reduction, we consider several different SKA core baseline designs, changing: (i) the number of antenna stations; (ii) the number of dipoles per station; and also (iii) the distribution of baseline lengths. We find that a design with a reduced number of dipoles per core station (increased field of view and total number of core stations), together with shortened baselines, maximises the recovered EoR signal. With this optimal baseline design, we investigate three observing strategies, analysing the trade-off between lowering the instrumental thermal noise against increasing the field of view. SKA1-low intends to perform a three tiered observing approach, including a deep 100 deg$^{2}$ at 1000 h, a medium-deep 1000 deg$^{2}$ at 100 h and a shallow 10,000 deg$^{2}$ at 10 h survey. We find that the three observing strategies result in comparable ($\lesssim$ per cent) constraints on our EoR astrophysical parameters. This is contrary to naive predictions based purely on the total signal-to-noise, thus highlighting the need to use EoR parameter constraints as a figure of merit, in order to maximise scientific returns with next generation interferometers.

The Ly$\alpha$ line of high-redshift galaxies has emerged as a powerful probe of both early galaxy evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-field survey with the Subaru Hyper Suprime-Cam (HSC), we study the angular correlation function (ACF) of narrow-band selected, $z\approx7$ Ly$\alpha$ emitting galaxies (LAEs). The clustering of LAEs is determined by both: (i) their typical host halo masses, $\bar{M}_{\rm h}$; and (ii) the absorption due to a patchy EoR, characterized by an average neutral fraction of the intergalactic medium, $\bar{x}_{\rm HI}$. We bracket the allowed LAE ACFs by exploring extreme scenarios for both the intrinsic Ly$\alpha$ emission and the large-scale pattern (i.e. morphology) of cosmic ionized patches in physical EoR models. Independent of the EoR morphology, current $z\approx 7$ ACF measurements constrain $\bar{x}_{\rm HI}\lesssim 0.5$ (1-$\sigma$). We also find that the low values of the currently-observed ACF imply that LAEs are hosted by relatively small dark matter halos: $\bar{M}_{\rm h} \lesssim 10^{10} M_{\odot}$, with corresponding duty cycles of $\lesssim$ few per cent. These values are over an order of magnitude lower than the analogous ones for color-selected, Lyman break galaxies, suggesting that $z\approx 7$ narrow-band LAEs searches are preferentially selecting young, star-burst galaxies, residing in less massive halos. The upcoming Ultra Deep campaign with the HSC will significantly improve constraints on both the EoR and LAE host halos.

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 investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic HII patches; and (ii) as reservoirs of neutral gas they contribute a non-zero 21-cm signal in otherwise ionized regions. Both effects damp the contrast between neutral and ionized patches during reionization, making detection of the epoch of reionization with 21-cm interferometry more challenging. Here we systematically investigate these effects using the latest semi-numerical simulations. We find that sinks dramatically suppress the peak in the redshift evolution of the variance, corresponding to the midpoint of reionization. As previously predicted, skewness changes sign at midpoint, but the fluctuations in the residual HI suppress a late-time rise. Furthermore, large levels of residual HI dramatically alter the evolution of the variance, skewness and power spectrum from that seen at lower levels. In general, the evolution of the large-scale modes provides a better, cleaner, higher signal-to-noise probe of reionization.

We present a framework for high-redshift ($z \geq 7$) galaxy formation that traces their dark matter (DM) and baryonic assembly in four cosmologies: Cold Dark Matter (CDM) and Warm Dark Matter (WDM) with particle masses of $m_x =$ 1.5, 3 and 5 ${\rm keV}$. We use the same astrophysical parameters regulating star formation and feedback, chosen to match current observations of the evolving ultra violet luminosity function (UV LF). We find that the assembly of observable (with current and upcoming instruments) galaxies in CDM and $m_x \geq 3 {\rm keV}$ WDM results in similar halo mass to light ratios (M/L), stellar mass densities (SMDs) and UV LFs. However the suppression of small-scale structure leads to a notably delayed and subsequently more rapid stellar assembly in the $1.5 {\rm keV}$ WDM model. Thus galaxy assembly in $m_x \leq 2 {\rm keV}$ WDM cosmologies is characterized by: (i) a dearth of small-mass halos hosting faint galaxies; and (ii) a younger, more UV bright stellar population, for a given stellar mass. The higher M/L ratio (effect ii) partially compensates for the dearth of small-mass halos (effect i), making the resulting UV LFs closer to CDM than expected from simple estimates of halo abundances. We find that the redshift evolution of the SMD is a powerful probe of the nature of DM. Integrating down to a limit of $M_{UV} =-16.5$ for the James Webb Space Telescope (JWST), the SMD evolves as $\log$(SMD)$\propto -0.63 (1+z)$ in $m_x = 1.5 {\rm keV}$ WDM, as compared to $\log$(SMD)$\propto -0.44 (1+z)$ in CDM. Thus high-redshift stellar assembly provides a powerful testbed for WDM models, accessible with the upcoming JWST.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

We introduce 21CMMC: a parallelized, Monte Carlo Markov Chain analysis tool, incorporating the epoch of reionization (EoR) seminumerical simulation 21CMFAST. 21CMMC estimates astrophysical parameter constraints from 21 cm EoR experiments, accommodating a variety of EoR models, as well as priors on model parameters and the reionization history. To illustrate its utility, we consider two different EoR scenarios, one with a single population of galaxies (with a mass-independent ionizing efficiency) and a second, more general model with two different, feedback-regulated populations (each with mass-dependent ionizing efficiencies). As an example, combining three observations (z=8, 9 and 10) of the 21 cm power spectrum with a conservative noise estimate and uniform model priors, we find that interferometers with specifications like the Low Frequency Array/Hydrogen Epoch of Reionization Array (HERA)/Square Kilometre Array 1 (SKA1) can constrain common reionization parameters: the ionizing efficiency (or similarly the escape fraction), the mean free path of ionizing photons and the log of the minimum virial temperature of star-forming haloes to within 45.3/22.0/16.7, 33.5/18.4/17.8 and 6.3/3.3/2.4 per cent, ~$1\sigma$ fractional uncertainty, respectively. Instead, if we optimistically assume that we can perfectly characterize the EoR modelling uncertainties, we can improve on these constraints by up to a factor of ~few. Similarly, the fractional uncertainty on the average neutral fraction can be constrained to within $\lesssim10$ per cent for HERA and SKA1. By studying the resulting impact on astrophysical constraints, 21CMMC can be used to optimize (i) interferometer designs; (ii) foreground cleaning algorithms; (iii) observing strategies; (iv) alternative statistics characterizing the 21 cm signal; and (v) synergies with other observational programs.

We present an overview of the theory of high-redshift star and X-ray source formation, and how they affect the 21-cm background. Primary focus is given to Lyman alpha pumping and X-ray heating mechanisms at cosmic dawn, opening a new observational window for high-redshift astrophysics by generating sizable fluctuations in the 21-cm background. We describe observational prospects for power spectrum analysis and 3D tomography (imaging) of the signature of these early astrophysical sources by SKA1-LOW and SKA2.

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

In this chapter we provide an overview of the current status of the simulations and modelling of the Cosmic Dawn and Epoch of Reionization. We discuss the modelling requirements as dictated by the characteristic scales of the problem and the SKA instrumental properties and the planned survey parameters. Current simulations include most of the relevant physical processes. They can follow the full nonlinear dynamics and are now reaching the required scale and dynamic range, although small-scale physics still needs to be included at sub-grid level. However, despite a significant progress in developing novel numerical methods for efficient utilization of current hardware they remain quite computationally expensive. In response, a number of alternative approaches, particularly semi-analytical/semi-numerical methods, have been developed. While necessarily more approximate, if appropriately constructed and calibrated on simulations they could be used to quickly explore the vast parameter space available. Further work is still required on including some physical processes in both simulations and semi-analytical modelling. This hybrid approach of fast, approximate modelling calibrated on numerical simulations can then be used to construct large libraries of reionization models for reliable interpretation of the observational data.

We provide an overview of 21cm tomography of the Cosmic Dawn and Epoch of Reionization as possible with SKA-Low. We show why tomography is essential for studying CD/EoR and present the scales which can be imaged at different frequencies for the different phases of SKA- Low. Next we discuss the different ways in which tomographic data can be analyzed. We end with an overview of science questions which can only be answered by tomography, ranging from the characterization of individual objects to understanding the global processes shaping the Universe during the CD/EoR

The Square Kilometre Array (SKA) will offer an unprecedented view onto the early Universe, using interferometric observations of the redshifted 21cm line. The 21cm line probes the thermal and ionization state of the cosmic gas, which is governed by the birth and evolution of the first structures in our Universe. Here we show how the evolution of the 21cm signal will allow us to study when the first generations of galaxies appeared, what were their properties, and what was the structure of the intergalactic medium. We highlight qualitative trends which will offer robust insights into the early Universe.

Observations of the redshifted 21 cm line from neutral hydrogen will open a new window on the early Universe. By influencing the thermal and ionization history of the intergalactic medium (IGM), annihilating dark matter (DM) can leave a detectable imprint in the 21 cm signal. Building on the publicly available 21cmFAST code, we compute the 21 cm signal for a 10 GeV WIMP DM candidate. The most pronounced role of DM annihilations is in heating the IGM earlier and more uniformly than astrophysical sources of X-rays. This leaves several unambiguous, qualitative signatures in the redshift evolution of the large-scale ($k\approx0.1$ Mpc$^{-1}$) 21 cm power amplitude: (i) the local maximum (peak) associated with IGM heating can be lower than the other maxima; (ii) the heating peak can occur while the IGM is in emission against the cosmic microwave background (CMB); (iii) there can be a dramatic drop in power (a global minimum) corresponding to the epoch when the IGM temperature is comparable to the CMB temperature. These signatures are robust to astrophysical uncertainties, and will be easily detectable with second generation interferometers. We also briefly show that decaying warm dark matter has a negligible role in heating the IGM.

We present new upper limits on the volume-weighted neutral hydrogen fraction, <xHI>, at z~5-6 derived from spectroscopy of bright quasars. The fraction of the Lyman-alpha and Lyman-beta forests that is "dark" (with zero flux) provides the only model-independent upper limit on <xHI>, requiring no assumptions about the physical conditions in the intergalactic medium or the quasar's unabsorbed UV continuum. In this work we update our previous results using a larger sample (22 objects) of medium-depth (~ few hours) spectra of high-redshift quasars obtained with the Magellan, MMT, and VLT. This significantly improves the upper bound on <xHI> derived from dark pixel analysis to <xHI> <= 0.06 + 0.05 (1{\sigma}) at z=5.9, and <xHI> <= 0.04 + 0.05 (1{\sigma}) at z=5.6. These results provide robust constraints for theoretical models of reionization, and provide the strongest available evidence that reionization has completed (or is very nearly complete) by z~6.

The large cross-section of the Lyman alpha (Lya) line makes it a sensitive probe of the ionization state of the intergalactic medium (IGM). Here we present the most complete study to date of the IGM Lya opacity, and its application to the redshift evolution of the 'Lya fraction', i.e. the fraction of color-selected galaxies with a detectable Lya emission line. We use a tiered approach, which combines large-scale semi-numeric simulations of reionization with moderate-scale hydrodynamic simulations of the ionized IGM. This allows us to simultaneously account for evolution in both: (i) the opacity from an incomplete (patchy) reionization, parameterized by the filling factor of ionized regions, Q_HII; and (ii) the opacity from self-shielded systems in the ionized IGM, parameterized by the average photo-ionization rate inside HII regions, \Gamma. In contrast to recent empirical models, attenuation from patchy reionization has a unimodal distribution along different sightlines, while attenuation from self-shielded systems is more bimodal. We quantify the average IGM transmission in our (Q_HII, \Gamma) parameter space, which can easily be used to interpret new data sets. Using current observations, we predict that the Lya fraction cannot drop by more than a factor of ~2 with IGM attenuation alone, even for HII filling factors as low as Q_HII>0.1. Larger changes in the Lya fraction could result from a co-evolution with galaxy properties. Marginalizing over \Gamma, we find that current observations constrain Q_HII < 0.6 at z=7 [68% confidence level (C.L.)]. However, all of our parameter space is consistent with observations at 95% C.L., highlighting the need for larger observational samples at z >= 6.

The cosmological 21cm signal is a physics-rich probe of the early Universe, encoding information about both the ionization and the thermal history of the intergalactic medium (IGM). The latter is likely governed by X-rays from star-formation processes inside very high redshift (z > 15) galaxies. Due to the strong dependence of the mean free path on the photon energy, the X-ray SED can have a significant impact on the interferometric signal from the cosmic dawn. Recent Chandra observations of nearby, star-forming galaxies show that their SEDs are more complicated than is usually assumed in 21cm studies. In particular, these galaxies have ubiquitous, sub-keV thermal emission from the hot interstellar medium (ISM), which generally dominates the soft X-ray luminosity (with energies < 1 keV, sufficiently low to significantly interact with the IGM). Using illustrative soft and hard SEDs, we show that the IGM temperature fluctuations in the early Universe would be substantially increased if the X-ray spectra of the first galaxies were dominated by the hot ISM, compared with X-ray binaries with harder spectra. The associated large-scale power of the 21cm signal would be higher by roughly a factor of three. More generally, we show that the peak in the redshift evolution of the large-scale (k = 0.2 1/Mpc) 21cm power is a robust probe of the soft-band SED of the first galaxies, and importantly, is not degenerate with their bolometric luminosities. On the other hand, the redshift of the peak constrains the X-ray luminosity and halo masses which host the first galaxies.

We describe a new technique for constraining the radio loud population of active galactic nuclei at high redshift by measuring the imprint of 21 cm spectral absorption features (the 21 cm forest) on the 21 cm power spectrum. Using semi-numerical simulations of the intergalactic medium and a semi-empirical source population we show that the 21 cm forest dominates a distinctive region of $k$-space, $k \gtrsim 0.5 \text{Mpc}^{-1}$. By simulating foregrounds and noise for current and potential radio arrays, we find that a next generation instrument with a collecting area on the order of $\sim 0.1\text{km}^2$ (such as the Hydrogen Epoch of Reionization Array) may separately constrain the X-ray heating history at large spatial scales and radio loud active galactic nuclei of the model we study at small ones. We extrapolate our detectability predictions for a single radio loud active galactic nuclei population to arbitrary source scenarios by analytically relating the 21 cm forest power spectrum to the optical depth power spectrum and an integral over the radio luminosity function.

The nature of the seeds of high-redshift supermassive black holes (SMBHs) is a key question in cosmology. Direct collapse black holes (DCBH) that form in pristine, atomic-line cooling halos, illuminated by a Lyman-Werner (LW) UV flux exceeding a critical threshold J_crit, represent an attractive possibility. We investigate when and where these conditions are met during cosmic evolution. For the LW intensity, J_LW, we account for departures from the background value in close proximity to star forming galaxies. For the pristine halo fraction, we account for both (i) supernova driven outflows, and (ii) the inherent pollution from progenitor halos. We estimate the abundance of DCBH formation sites, n_DCBH(z), and find that it increases with cosmic time from n_DCBH(z=20) ~ 1e-12 -1e-7 cMpc^-3 to n_DCBH(z=10) ~ 1e-10 - 1e-5 cMpc^-3. Our analysis shows the possible importance of galactic winds, which can suppress the predicted n_DCBH by several orders of magnitude, and cause DCBH formation to preferentially occur around the UV-brightest (M_UV ~ -22 to -20) star forming galaxies. Our analysis further highlights the dependence of these predictions on (i) the escape fraction of LW photons, (ii) J_crit, and (iii) the galactic outflow prescription.

[ABRIDGED] The cosmological 21cm signal is set to become the most powerful probe of the early Universe, with first generation interferometers aiming to make statistical detections of reionization. There is increasing interest also in the pre-reionization epoch when the intergalactic medium was heated by an early X-ray background. Here we perform parameter studies varying the halo masses hosting galaxies, and their X-ray production efficiencies. We also relate these to popular models of Warm Dark Matter cosmologies. For each parameter combination we compute the signal-to-noise (S/N) of the large-scale (k~0.1/Mpc) 21cm power for both reionization and X-ray heating for a 2000h observation with several instruments: 128 tile Murchison Wide Field Array (MWA128T), a 256 tile extension (MWA256T), the Low Frequency Array (LOFAR), the 128 element Precision Array for Probing the Epoch of Reionization (PAPER), and the second generation Square Kilometre Array (SKA). We show that X-ray heating and reionization in many cases are of comparable detectability. For fiducial astrophysical parameters, MWA128T might detect X-ray heating thanks to its extended bandpass. When it comes to reionization, both MWA128T and PAPER will also only achieve marginal detections, unless foregrounds on larger scales can be mitigated. On the other hand, LOFAR should detect plausible models of reionization at S/N > 10. The SKA will easily detect both X-ray heating and reionization.

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.

By depleting the ionizing photon budget available to expand cosmic HII regions, recombining systems (or Lyman limit systems) can have a large impact during (and following) cosmic reionization. Unfortunately, directly resolving such structures in large-scale reionization simulations is computationally impractical. Instead, here we implement a sub-grid prescription for tracking inhomogeneous recombinations in the intergalactic medium. Building on previous work parameterizing photo-heating feedback on star-formation, we present large-scale, semi-numeric reionization simulations which self-consistently track the local (sub-grid) evolution of both sources and sinks of ionizing photons. Our simple, single-parameter model naturally results in both an extended reionization and a modest, slowly-evolving emissivity, consistent with observations. Recombinations are instrumental in slowing the growth of large HII regions, and damping the rapid rise of the ionizing background in the late stages of (and following) reionization. As a result, typical HII regions are smaller by factors of $\sim 2-3$ throughout reionization. The large-scale ($k\lesssim 0.2\text{ Mpc$^{-1}$}$) ionization power spectrum is suppressed by factors of $\gtrsim 2-3$ in the second half of reionization. Therefore properly modeling recombinations is important in interpreting virtually all reionization observables, including upcoming interferometry with the redshifted 21 cm line. Consistent with previous works, we find the clumping factor of ionized gas to be $C_{\rm HII}\sim 4$ at the end of reionization.

We investigate the effects of warm dark matter (WDM) on the cosmic 21-cm signal. If dark matter exists as WDM instead of cold dark matter (CDM), its non-negligible velocities can inhibit the formation of low-mass halos that normally form first in CDM models, therefore delaying star-formation. The absence of early sources delays the build-up of UV and X-ray backgrounds that affect the 21-cm radiation signal produced by neutral hydrogen. With use of the 21CMFAST, code, we demonstrate that the pre-reionization 21-cm signal can be changed significantly in WDM models with a free-streaming length equivalent to that of a thermal relic with mass mx of up to ~10-20 keV. In such a WDM cosmology, the 21-cm signal traces the growth of more massive halos, resulting in a delay of the 21-cm absorption signature and followed by accelerated X-ray heating. CDM models where astrophysical sources have a suppressed photon-production efficiency can delay the 21-cm signal as well, although its subsequent evolution is not as rapid as compared to WDM. This motivates using the gradient of the global 21-cm signal to differentiate between some CDM and WDM models. Finally, we show that the degeneracy between the astrophysics and mx can be broken with the 21-cm power spectrum, as WDM models should have a bias-induced excess of power on large scales. This boost in power should be detectable with current interferometers for models with mx < 3 keV, while next generation instruments will easily be able to measure this difference for all relevant WDM models.

We present a fully self-consistent simulation of a synthetic survey of the furthermost cosmic explosions. The appearance of the first generation of stars (Population III) in the Universe represents a critical point during cosmic evolution, signaling the end of the dark ages, a period of absence of light sources. Despite their importance, there is no confirmed detection of Population III stars so far. A fraction of these primordial stars are expected to die as pair-instability supernovae (PISNe), and should be bright enough to be observed up to a few hundred million years after the big bang. While the quest for Population III stars continues, detailed theoretical models and computer simulations serve as a testbed for their observability. With the upcoming near-infrared missions, estimates of the feasibility of detecting PISNe are not only timely but imperative. To address this problem, we combine state-of-the-art cosmological and radiative simulations into a complete and self-consistent framework, which includes detailed features of the observational process. We show that a dedicated observational strategy using $\lesssim 8$ per cent of total allocation time of the James Webb Space Telescope mission can provide us up to $\sim 9-15$ detectable PISNe per year.

We propose a novel use of high-redshift galaxies, discovered in deep Hubble Space Telescope (HST) fields around strong lensing clusters. These fields probe small comoving volumes (about 1000 cubic Mpc) at high magnification ({\mu} > 10), and can detect otherwise inaccessible ultra-faint galaxies. Even a few galaxies found in such small volumes require a very high number density of collapsed dark matter (DM) halos. This implies significant primordial power on small scales, allowing these observations to rule out popular alternatives to standard cold dark matter (CDM) models, such as warm dark matter (WDM). In this work, we analytically compute WDM halo mass functions at z = 10, including the effects of both particle free-streaming and residual velocity dispersion. The preliminary number density corresponding to the two galaxies at z about 10 already detected by the Cluster Lensing And Supernova survey with Hubble (CLASH) constrains the WDM particle mass to m_x > 1 (0.9) keV at 84% (99.9%) confidence limit (for a thermal relic relativistic at decoupling). This limit depends only on the WDM halo mass function and, unlike previous constraints on m_x, is independent of any astrophysical modeling. The forthcoming HST Frontier Fields can significantly tighten these constraints.

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

Recent He II Lyman-alpha forest observations from 2.0 < z < 3.2 show large fluctuations in the optical depth at z > 2.7. These results point to a fluctuating He-ionizing background, which may be due to the end of helium reionization of this era. We present a fast, semi-numeric procedure to approximate detailed cosmological simulations. We compute the distribution of dark matter halos, ionization state of helium, and density field at z = 3 in broad agreement with recent simulations. Given our speed and flexibility, we investigate a range of ionizing source and active quasar prescriptions. Spanning a large area of parameter space, we find order-of-magnitude fluctuations in the He II ionization rate in the post-reionization regime. During reionization, the fluctuations are even stronger and develop a bimodal distribution, in contrast to semi-analytic models and the hydrogen equivalent. These distributions indicate a low-level ionizing background even at significant He II fractions.

[abridged] With their long mean free paths and efficient heating of the intergalactic medium (IGM), X-rays could have a dramatic impact on the thermal and ionization history of the Universe. We explore this in various signals: (i) Reionization history: including X-rays results in an earlier, more extended reionization. Efficient thermal feedback from X-ray heating could yield an extended, ~10% ionized epoch. (ii) Reionization morphology: a sizable (~10%) contribution of X-rays to reionization results in a more uniform morphology, though the impact is modest when compared at the same global neutral fraction, xH. However, changes in morphology cannot be countered by increasing the bias of the ionizing sources, making them a robust signature. (iii) The kinetic Sunyaev-Zel'dovich (kSZ) effect: at a fixed reionization history, X-rays decrease the kSZ power at l=3000 by ~0.5 microK^2. Our extreme model in which X-rays dominate reionization is the only one that is marginally consistent with upper limits from the South Pole Telescope, assuming no thermal Sunyaev-Zel'dovich (tSZ) - dusty galaxy correlation. Since this extreme model is unlikely, we conclude that there should be a sizable tSZ-dusty galaxy signal. (iv) The cosmic 21cm signal: the impact of X-rays on the 21cm power spectrum during the advanced stages of reionization (xH<0.7) is modest, except in extreme, X-ray dominated models. The largest impact of X-rays is to govern IGM heating. In fact, unless thermal feedback is efficient, the epoch of X-ray heating likely overlaps with the beginning of reionization (xH>0.9). This results in a 21cm power spectrum which is ~ 10-100 times higher than obtained from naive estimates ignoring this overlap. However, if thermal feedback is efficient, the resulting extended epoch between X-ray heating and reionization could provide a clean probe of the matter power spectrum in emission.

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^\alpha$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $\alpha$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

An ionizing UV background (UVB) inhibits gas accretion and photo-evaporates gas from the shallow potential wells of small, dwarf galaxies. During cosmological reionization, this effect can result in negative feedback: suppressing star-formation inside HII regions, thus impeding their continued growth. It is difficult to model this process, given the enormous range of scales involved. We tackle this problem using a tiered approach: combining parameterized results from single-halo collapse simulations with large-scale models of reionization. In the resulting reionization models, the ionizing emissivity of galaxies depends on the local values of the reionization redshift and the UVB intensity. We present a physically-motivated analytic expression for the average minimum mass of star-forming galaxies, which can be readily used in modeling galaxy formation. We find that UVB feedback: (i) delays the end stages of reionization by less than 0.5 in redshift; (ii) results in a more uniform distribution of HII regions, peaked on smaller-scales (with large-scale ionization power suppressed by tens of percent); and (iii) suppresses the global photoionization rate per baryon by a factor of < 2 towards the end of reionization. However, the impact is modest, since the hydrodynamic response of the gas to the UVB occurs on a time-scale comparable to reionization. In particular, the popular approach of modeling UVB feedback with an instantaneous transition in the minimum mass of star-forming galaxies, dramatically overestimates its importance. UVB feedback does not significantly affect reionization unless: (i) molecularly-cooled galaxies contribute significantly to reionization; or (ii) internal feedback processes strongly couple with UVB feedback in the early Universe. Since both are considered unlikely, we conclude that there is no significant self-regulation of reionization by UVB feedback.

The reionization of the intergalactic medium (IGM) was likely inhomogeneous and extended. By heating the IGM and photo-evaporating gas from the outskirts of galaxies, this process can have a dramatic impact on the growth of structures. Using a suite of spherically-symmetric collapse simulations spanning a large parameter space, we study the impact of an ionizing ultraviolet background (UVB) on the condensation of baryons onto dark matter halos. We present an expression for the halo baryon fraction, which is an explicit function of: (i) halo mass; (ii) UVB intensity; (iii) redshift; (iv) redshift at which the halo was exposed to a UVB. We also present a corresponding expression for the characteristic or critical mass, defined as the halo mass which retains half of its baryons compared to the global value. Since our results are general and physically-motivated, they can be broadly applied to inhomogeneous reionization models.

We study the imprint of dark matter (DM) annihilation on the global 21 cm signal from the Dark Ages to Cosmic Reionization. Motivated by recent observations, we focus on three DM candidates: (i) a 10 GeV Bino-like neutralino (ii) a 200 GeV Wino and (iii) a 1 TeV heavier particle annihilating into leptons. For each DM candidate we assume two values for the thermally averaged annihilation cross section \sigma v, the standard thermal value \sigma v_th = 3 x 10^-26 cm^3 s^-1 and the maximum value allowed by WMAP7 data, \sigma v_max. We include the enhancement of DM annihilations due to collapsed structures, detailed estimates of energy deposition into the intergalactic medium (IGM), as well realistic prescriptions for astrophysical sources of UV and X-ray radiation. In these models, the additional heat input from DM annihilation suppresses the mean 21cm brightness temperature offset by \delta Tb ~ few-100 mK. In particular, the very deep \delta Tb ~ -150 mK absorption feature at ~20 < z < 25 predicted by popular models of the first galaxies is considerably reduced or totally erased by some of the considered DM candidates. Such an enhancement in IGM heating could come from either DM annihilations or a stronger-than-expected astrophysical component (i.e. abundant early X-ray sources). However, we find that the two signatures are not degenerate, since the DM heating is dominated by halos several orders of magnitude smaller than those hosting galaxies, whose fractional abundance evolves more slowly resulting in a smaller gradient: d \delta Tb/d \nu < 4 mK/MHz in the range \nu ~60-80 MHz. The detection of such signals by future radio telescopes would be clear evidence of DM energy injection at high-redshifts.

The spectra of several high-redshift (z>6) quasars have shown evidence for a Gunn-Peterson (GP) damping wing, indicating a substantial mean neutral hydrogen fraction (x_HI > 0.03) in the z ~ 6 intergalactic medium (IGM). However, previous analyses assumed that the IGM was uniformly ionized outside of the quasar's HII region. Here we relax this assumption and model patchy reionization scenarios for a range of IGM and quasar parameters. We quantify the impact of these differences on the inferred x_HI, by fitting the spectra of three quasars: SDSS J1148+5251 (z=6.419), J1030+0524 (z=6.308), and J1623+3112 (z=6.247). We find that the best-fit values of x_HI in the patchy models agree well with the uniform case. More importantly, we confirm that the observed spectra favor the presence of a GP damping wing, with peak likelihoods decreasing by factors of > few - 10 when the spectra are modeled without a damping wing. We also find that the Ly alpha absorption spectra, by themselves, cannot distinguish the damping wing in a relatively neutral IGM from a damping wing in a highly ionized IGM, caused either by an isolated neutral patch, or by a damped Ly alpha absorber (DLA). However, neutral patches in a highly ionized universe (x_HI < 0.01), and DLAs with the large required column densities (N_HI > few x 10^{20} cm^{-2}) are both rare. As a result, when we include reasonable prior probabilities for the line of sight (LOS) to intercept either a neutral patch or a DLA at the required distance of ~ 40-60 comoving Mpc away from the quasar, we find strong lower limits on the neutral fraction in the IGM, x_HI > 0.1 (at 95% confidence). This strengthens earlier claims that a substantial global fraction of hydrogen in the z~6 IGM is in neutral form.

[ABRIDGED] Inhomogeneous reionization acts as a source of arcminute-scale anisotropies in the cosmic microwave background (CMB), the most important of which is the kinetic Sunyaev-Zel'dovich (kSZ) effect. Observational efforts with the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) are poised to detect this signal for the first time. Indeed, recent SPT measurements place a bound on the dimensionless kSZ power spectrum around a multipole of l=3000 of P_tot < 2.8 (6) micro K^2 at 95% C.L., by ignoring (allowing) correlations between the thermal Sunyaev-Zel'dovich (tSZ) effect and the cosmic infrared background (CIB). To interpret these and upcoming observations, we compute the kSZ signal from a suite of ~ 100 reionization models using the publicly available code 21cmFAST. Our physically motivated reionization models are parameterized by the ionizing efficiency of high-redshift galaxies, the minimum virial temperature of halos capable of hosting stars, and the ionizing photon mean free path. We predict the contribution of patchy reionization to be P_patchy = 1.5-3.5 micro K^2. Therefore, even when conservatively adopting a low estimate of the post-reionization signal, P_OV ~ 2 micro K^2, none of our models are consistent with the aggressive 2sigma SPT bound that does not include correlations. This implies that either: (i) the early stages of reionization occurred in a much more homogeneous manner than suggested by the stellar-driven scenarios we explore, such as would be the case if, e.g., very high energy X-rays or exotic particles contributed significantly; and/or (ii) that there is a significant correlation between the CIB and the tSZ. On the other hand, the conservative SPT bound is compatible with all of our models, and is on the boarder of constraining reionization.

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.

The first generations of astrophysical objects made a substantial impact on our Universe with their radiation. X-rays from the first sources, with large mean free paths, likely quickly heated the intergalactic medium (IGM). The second generation of 21cm instruments can provide a unique view into this early epoch. The early stages of reionization likely followed, driven by so-called "minihalos", i.e. molecularly-cooled halos. These small halos were susceptible to complex feedback mechanisms, especially from the soft-UV background which preceded reionization, resulting in complex and possibly extended early stages of reionization. When atomically-cooled galaxies emerged as the dominant ionizers, reionization could proceed more rapidly, with these being less sensitive to radiative feedback than previously thought. Reionization could have slowed in the final stages when the ionized bubbles grew larger than the separation of Lyman limit absorption systems (LLSs). The final stages likely involved the photo-evaporation of LLSs, which by then regulated the rise of the UV background. I discuss the theoretical underpinnings of this narrative, as well as how future 21cm observations may help shed light on the outstanding uncertainties.

Cosmic reionization is expected to be complex, extended and very inhomogeneous. Existing constraints at z~6 on the volume-averaged neutral hydrogen fraction, <x_HI>, are highly model-dependent and controversial. Constraints at z<6, suggesting that the Universe is highly ionized, are also model-dependent, but more fundamentally are invalid in the context of inhomogeneous reionization. As such, it has recently been pointed out that there is no conclusive evidence that reionization has completed by z~5--6, a fact that has important ramifications on the interpretation of high-redshift observations and theoretical models. We present the first direct upper limits on <x_HI> at z~5--6 using the simple and robust statistic of the covering fraction of dark pixels in the Ly alpha/beta forests of high redshift quasars. With a sample of 13 Keck ESI spectra we constrain <x_HI> < 0.2 at 5 < z < 5.5, rising to <x_HI> < 0.8 at z~6.1. We also find tentative evidence for a break in the redshift evolution of the dark covering fraction at z~5.5. A subsample of two deep spectra provides a more stringent constraint of <x_HI>(z=6.1) < 0.5 when combined with conservative estimates of cosmic variance. This upper limit is comparable to existing results at z~6 but is more robust. The results presented here do not rely on assumptions about the quasar continuum, IGM density, HII morphology or ionizing background fields, and thus are a good starting point for future interpretation of high redshift observations.

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

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc/h simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative application of a smoothing filter to the underlying source and density fields. We compare these algorithms' resulting ionization and 21 cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement with each other (with the cross-correlation coefficient of the ionization fields >0.8 for k < 10 h/Mpc and in good agreement with the analytic schemes (>0.6 for k < 1 h/Mpc). When used to predict the 21cm power spectrum at different times during reionization, all ionization algorithms agree with one another at the 10s of percent level. This agreement suggests that the different approximations involved in the ray tracing algorithms are sensible and that semi-numerical schemes provide a numerically-inexpensive, yet fairly accurate, description of the reionization process.

It is generally taken for granted that reionization has completed by z=6, due to the detection of flux in the Lyman alpha forest at redshifts z<6. However, since reionization is expected to be highly inhomogeneous, much of the spectra pass through just the ionized component of the intergalactic medium (IGM) even for non-negligible values of the volume-weighted mean neutral hydrogen fraction, x_HI. We study the expected signature of an incomplete reionization at z ~ 5--6, using very large-scale (2 Gpc) seminumeric simulations. We find that ruling out an incomplete reionization is difficult at these redshifts since: (1) quasars reside in biased regions of the ionization field, with fewer surrounding HI patches than implied by the global mean, x_H; (2) absorption from the residual neutral hydrogen inside the ionized IGM generally dominates over the absorption from the remaining HI regions; (3) modeling the Lyman alpha forest and its redshift evolution even in just the ionized IGM is very difficult, and nearly impossible to do a priori. We propose using the fraction of pixels which are dark as a simple, nearly model-independent upper limit on x_HI. Alternately, the size distribution of regions with no detectable flux (dark gaps) can be used to place a more model dependent constraint. Either way, the current sample of quasars is statistically insufficient to constrain x_HI at z~6 to even the 10 per cent level. At z~5, where there are more available sightlines and the forest is less dark, constraining x_HI < 0.1 might be possible given a large dynamic range from very deep spectra and/or the Lyman beta forest. We conclude with the caution against over-interpreting the observations. There is currently no direct evidence that reionization was complete by z ~ 5-6.

UV radiation from early astrophysical sources could have a large impact on subsequent star formation in nearby protogalaxies. Here we study the radiative feedback from the first, short-lived stars using hydrodynamical simulations with transient UV backgrounds (UVBs) and persistent Lyman-Werner backgrounds (LWBs) of varying intensity. We extend our prior work in Mesinger et al. (2006), by studying a more typical region whose proto-galaxies form at lower redshifts, z~13-20, in the epoch likely preceding the bulk of reionization. We confirm our previous results that feedback in the relic HII regions resulting from such transient radiation, is itself transient. Feedback effects dwindle away after ~30% of the Hubble time, and the same critical specific intensity of J_UV~0.1 x 10^{-21} ergs/s/cm^2/Hz/sr separates positive and negative feedback regimes. Additionally, we discover a second episode of eventual positive feedback in halos which have not yet collapsed when their progenitor regions were exposed to the transient UVB. This eventual positive feedback appears in all runs, regardless of the strength of the UVB. However, this feedback regime is very sensitive to the presence of Lyman-Werner radiation, and notable effects disappear under fairly modest background intensities of J_LW>10^{-3} x 10^{-21} ergs/s/cm^2/Hz/sr. We conclude that UV radiative feedback in relic HII regions, although a complicated process, seems unlikely to have a major impact on the progress of cosmological reionization, provided that present estimates of the lifetime and luminosity of a PopIII star are accurate. More likely is that the build-up of the LWB ultimately governs the feedback strength until a persistent UV background can be established. [abridged]

We study the spatial fluctuations in the hydrogen ionizing background in the epoch following reionization (z ~ 5--6). The rapid decrease with redshift in the photon mean free path (m.f.p.), combined with the clustering of increasingly rare ionizing sources, can result in a very inhomogenous ionizing background during this epoch. We systematically investigate the probability density functions (PDFs) and power spectra of ionizing flux, by varying several parameters such as the m.f.p., minimum halo mass capable of hosting stars, and halo duty cycle. In order to be versatile, we make use of analytic, semi-numeric and numeric approaches. Our models show that the ionizing background indeed has sizable fluctuations during this epoch sourced by the clustering of sources, with the PDFs being a factor of few wide at half of the maximum likelihood. The distributions also show marked asymmetries, with a high-value tail set by clustering on small scales, and a shorter low-value tail which is set by the mean free path. The power spectrum of the ionizing background is much more sensitive to source properties than the PDF and can be well-understood analytically with a framework similar to the halo model (usually used to describe dark matter clustering). Nevertheless, we find that Lya forest spectra are extremely insensitive to the details of the UVB, despite marked differences in the PDFs and power spectra of our various ionizing backgrounds. Assuming a uniform ionizing background only underestimates the value of the mean ionization rate inferred from the Lya forest by a few percent. Instead, analysis of the Lya forest is dominated by the uncertainties in the density field. Thus, our results justify the common assumption of a uniform ionizing background in Lya forest analysis even during this epoch.

We study the gas content of halos in the early universe using high resolution hydrodynamical simulations. We extract from the simulations and also predict based on linear theory the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was $\sim 3 \times 10^4$ M$_\odot$. We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20% accuracy the halo must be resolved into at least 500 dark matter particles.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

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.

One of the most sought-after signatures of reionization is a rapid increase in the ionizing background (usually measured through the Lyman-alpha optical depth toward distant quasars). Conventional wisdom associates this with the "overlap" phase when ionized bubbles merge, allowing each source to affect a much larger volume. We argue that this picture fails to describe the transition to the post-overlap Universe, where Lyman-limit systems absorb ionizing photons over moderate lengthscales (20-100 Mpc). Using an analytic model, we compute the probability distribution of the amplitude of the ionizing background throughout reionization, including both discrete ionized bubbles and Lyman-limit systems (parameterized by an attenuation length). We show that overlap does not by itself cause a rapid increase in the ionizing background or a rapid decrease in the mean Lyman-alpha transmission toward distant quasars. More detailed semi-numeric models support these conclusions. We argue that rapid changes should instead be interpreted as evolution in the attenuation length itself, which may or may not be directly related to overlap.

The ionizing ultraviolet background (UVB) during reionization can suppress the gas content of low-mass galaxies, even those capable of efficient atomic cooling, and thus lead to an extended reionization epoch. In this work, we explore the importance of negative UV radiative feedback on Tvir > 10^4 K halos during the middle and late stages of reionization. We do not try to self-consistently model reionization; instead, we explore a large parameter space in an attempt to draw general, robust conclusions. We do this using a tiered approach. Using 1-D hydrodynamical simulations, we model the collapse of gas onto halos of various masses under UVBs of various intensities. We then generate realistic, parametrized maps of the inhomogeneous UVB, using large-scale semi-numeric simulations. By combining these results, we find that under all reasonably conservative scenarios, UV feedback on atomically-cooled halos is not strong enough to notably delay the bulk of reionization. Such a delay is only likely if ionizing efficiencies of z > 10 sources are much higher (~ two orders of magnitude) than z ~ 6 data seem to imply. We also find that feedback is very strongly dependent on halo mass. Our results suggest that the natural time-scale for the bulk of reionization is the growth of the global collapsed fraction contained in Tvir > 10^4 K halos. Finally, our results underscore the importance of taking into account extended dynamical ranges when modeling reionization.

We investigate the potential of exploiting Lya Emitters (LAEs) to constrain the volume-weighted mean neutral hydrogen fraction of the IGM, x_H, at high redshifts (specifically z~9). We use "semi-numerical'' simulations to efficiently generate density, velocity, and halo fields at z=9 in a 250 Mpc box, resolving halos with masses M>2.2e8 solar masses. We construct ionization fields corresponding to various values of x_H. With these, we generate LAE luminosity functions and "counts-in-cell'' statistics. As in previous studies, we find that LAEs begin to disappear rapidly when x_H > 0.5. Constraining x_H(z=9) with luminosity functions is difficult due to the many uncertainties inherent in the host halo mass <--> Lya luminosity mapping. However, using a very conservative mapping, we show that the number densities derived using the six z~9 LAEs recently discovered by Stark et al. (2007) imply x_H < 0.7. On a more fundamental level, these LAE number densities, if genuine, require substantial star formation in halos with M < 10^9 solar masses, making them unique among the current sample of observed high-z objects. Furthermore, reionization increases the apparent clustering of the observed LAEs. We show that a ``counts-in-cell'' statistic is a powerful probe of this effect, especially in the early stages of reionization. Specifically, we show that a field of view (typical of upcoming IR instruments) containing LAEs has >10% higher probability of containing more than one LAE in a x_H>0.5 universe than a x_H=0 universe with the same overall number density. With this statistic, a fully ionized universe can be robustly distinguished from one with x_H > 0.5 using a survey containing only ~ 20--100 galaxies.

One well-known way to constrain the hydrogen neutral fraction, x_H, of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lya absorption line. We examine this method's effectiveness in light of recent models showing that the IGM neutral fraction is highly inhomogeneous on large scales during reionization. Using both analytic models and "semi-numeric" simulations, we show that the "picket-fence" absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of x_H. In particular, we show that simple fits to the damping wing tend to overestimate the true neutral fraction in a partially ionized universe, with a fractional error of ~ 30% near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate x_H if the observational sample only probes a subset of the entire halo population, such as quasars with large HII regions. We also find that the damping wing absorption profile is generally steeper than one would naively expect in a homogeneously ionized universe. The profile steepens and the sightline-to-sightline scatter increases as reionization progresses. Of course, the bias and scatter also depend on x_H and so can, at least in principle, be used to constrain it. Damping wing constraints must therefore be interpreted by comparison to theoretical models of inhomogeneous reionization.

We present a method to construct semi-numerical ``simulations'', which can efficiently generate realizations of halo distributions and ionization maps at high redshifts. Our procedure combines an excursion-set approach with first-order Lagrangian perturbation theory and operates directly on the linear density and velocity fields. As such, the achievable dynamic range with our algorithm surpasses the current practical limit of N-body codes by orders of magnitude. This is particularly significant in studies of reionization, where the dynamic range is the principal limiting factor. We test our halo-finding and HII bubble-finding algorithms independently against N-body simulations with radiative transfer and obtain excellent agreement. We compute the size distributions of ionized and neutral regions in our maps. We find even larger ionized bubbles than do purely analytic models at the same volume-weighted mean hydrogen neutral fraction. We also generate maps and power spectra of 21-cm brightness temperature fluctuations, which for the first time include corrections due to gas bulk velocities. We find that velocities widen the tails of the temperature distributions and increase small-scale power, though these effects quickly diminish as reionization progresses. We also include some preliminary results from a simulation run with the largest dynamic range to date: a 250 Mpc box that resolves halos with masses M >~ 2.2 x10^8 M_sun. We show that accurately modeling the late stages of reionization requires such large scales. The speed and dynamic range provided by our semi-numerical approach will be extremely useful in the modeling of early structure formation and reionization.

We make use of hydrodynamical simulations of the intergalactic medium (IGM) to create model quasar absorption spectra. We compare these model spectra with the observed Keck spectra of three z>6.2 quasars with full Gunn-Peterson troughs: SDSS J1148+5251 (z=6.42), SDSS J1030+0524 (z=6.28), and SDSS J1623+3112 (z=6.22). We fit the probability density distributions (PDFs) of the observed Ly alpha optical depths with those generated from the simulation, by exploring a range of values for the size of the quasar's surrounding HII region, R_S, the volume-weighted mean neutral hydrogen fraction in the ambient IGM, x_H, and the quasar's ionizing photon emissivity, N_Q. In order to avoid averaging over possibly large sightline-to-sightline fluctuations in IGM properties, we analyze each observed quasar independently. We find the following results for J1148+5251, J1030+0524, and J1623+3112: The best-fit sizes R_S are 40, 41, and 29 (comoving) Mpc, respectively. These constraints are tight, with only ~ 10% uncertainties, comparable to those caused by redshift--determination errors. The best-fit values of N_Q are 2.1, 1.3, and 0.9 x 10^{57} s^{-1}, respectively, with a factor of ~ 2 uncertainty in each case. Finally, the best--fit values of x_H are 0.16, 1.0, and 1.0, respectively. The uncertainty in the case of J1148+5251 is large, and x_H is not well constrained. However, for both J1030+0524 and J1623+3112, we find a significant lower limit of x_H >~ 0.033. Our method is different from previous analyses of the GP absorption spectra of these quasars, and our results strengthen the evidence that the rapid end--stage of reionization is occurring near z ~ 6.

We use three-dimensional hydrodynamic simulations to investigate the effects of a transient photoionizing ultraviolet (UV) flux on the collapse and cooling of pregalactic clouds. These clouds have masses in the range 10^5 -10^7 M_sun, form at high redshifts (z>18), are assumed to lie within the short-lived cosmological HII regions around the first generation of stars. In addition, we study the combined effects of this transient UV flux and a persistent Lyman-Werner (LW) background from distant sources. In the absence of a LW background, we find that a critical specific intensity of J_UV ~ 0.1 x 10^-21 ergs s^-1 cm^-2 Hz^-1 sr^-1 demarcates a transition from net negative to positive feedback for the halo population. A weaker UV flux stimulates subsequent star formation inside the fossil HII regions, by enhancing the H_2 molecule abundance. A stronger UV flux significantly delays star-formation by reducing the gas density, and increasing the cooling time, at the centers of collapsing halos. At a fixed J_UV, the sign of the feedback also depends strongly on the density of the gas at the time of UV illumination. Regardless of the whether the feedback is positive or negative, we find that once the UV flux is turned off, its impact stars to diminish after ~30% of the Hubble time. In the more realistic case when a LW background is present, with J_LW > 0.01 x 10^-21 ergs s^-1 cm^-2 Hz^-1 sr^-1, strong suppression persists down to the lowest redshift (z=18) in our simulations. Finally, we find evidence that heating and photoevaporation by the transient UV flux renders the ~10^6 M_sun halos inside fossil HII regions more vulnerable to subsequent H_2 photo-dissociation by a LW background.

We develop and assess the potential of several powerful techniques, designed to investigate the details of reionization. First, we present a procedure to probe the neutral fraction, x_HI, using the Lyman alpha transmission statistics of high-redshift (z > 6) sources. We find that only tens of bright quasar spectra could distinguish between x_HI ~ 1 and x_HI < 0.01. A rudimentary application of such a technique on quasar SDSS J1030+0524 has yielded compelling evidence of a large neutral fraction (x_HI > 0.2) at z ~ 6. We also generate the observable, high-z supernovae (SNe) rates and quantify the prospects of detecting the suppression of star-formation in low-mass galaxies at reionization from such SNe rates, specifically from those obtainable from the James Webb Space Telescope (JWST). Our analysis suggests that searches for SNe could yield thousands of SNe per unit redshift at z ~ 6, and be a valuable tool at studying reionization features and feedback effects out to z < 13.

We model the number of detectable supernovae (SNe) as a function of redshift at different flux thresholds, making use of the observed properties of local SNe, such as their lightcurves, fiducial spectra, and peak magnitude distributions. We obtain a rate of 0.4-2.3 SNe/arcmin^2/yr at z > 5 at the near infrared (4.5 micron) flux density threshold of 3 nJy (achievable with the James Webb Space Telescope (JWST) in a 10^5 s integration). In a hypothetical one-year survey, it should be possible to detect up to several thousand SNe per unit redshift at z > 6. We discuss the possible application of such a large sample of distant SNe as a probe of the epoch of reionization. By heating the intergalactic medium and raising the cosmological Jeans mass, the process of reionization can suppress star formation in low-mass galaxies. This could have produced a relatively sharp drop in the SNR around the redshift of reionization. We quantify the detectability of this feature in future surveys of distant SNe by varying the redshift and duration of the feature, as well as its impact on the SFR in low-mass halos, which results in different redshifts, widths, and sizes of the drop in the expected SFR. We find that the drop can be detected out to z_re ~ 13, as long as (i) the reionization history contains a relatively rapid feature which is synchronized over different regions to within \Delta z < 1 - 3, (ii) the star-formation efficiency in halos that dominate reionization is ~ 10%, and (iii) reionization significantly suppresses the star formation in low-mass halos. Depending on the details of (i) - (iii), this could be achieved with a survey lasting less than two weeks. Detecting this signature would also help elucidate the feedback mechanism that regulates reionization.

Cosmological models that include suppression of the power spectrum of density fluctuations on small scales exhibit an exponential reduction of high-redshift, non-linear structures, including a reduction in the rate of gamma ray bursts (GRBs). Here we quantify the constraints that the detection of distant GRBs would place on structure formation models with reduced small-scale power. We compute the number of GRBs that could be detectable by the Swift satellite at high redshifts (z > 6), assuming that the GRBs trace the cosmic star formation history, which itself traces the formation of non-linear structures. We calibrate simple models of the intrinsic luminosity function of the bursts to the number and flux distribution of GRBs observed by the Burst And Transient Source Experiment (BATSE). We find that a discovery of high-z GRBs would imply strong constraints on models with reduced small-scale power. For example, a single GRB at z > 10, or 10 GRBs at z > 5, discovered by Swift during its scheduled two-year mission, would rule out an exponential suppression of the power spectrum on scales below R_c=0.09 Mpc (exemplified by warm dark matter models with a particle mass of m_x=2 keV). Models with a less sharp suppression of small-scale power, such as those with a red tilt or a running scalar index, n_s, are more difficult to constrain, because they are more degenerate with an increase in the power spectrum normalization, sigma_8, and with models in which star-formation is allowed in low-mass minihalos. We find that a tilt of \delta n_s ~ 0.1 is difficult to detect; however, an observed rate of 1 GRB/yr at z > 12 would yield an upper limit on the running of the spectral index, alpha = d(n_s)/d(ln k) > -0.05.

We discuss currently available observational constraints on the reionization history of the intergalactic medium (IGM), and the extent to which accreting black holes (BHs) can help explain these observations. We show new evidence, based on the combined statistics of Lyman alpha and beta absorption in quasar spectra, that the IGM contains a significant amount of neutral hydrogen, and is experiencing rapid ionization around redshift z=6. However, we argue that quasar BHs, even faint ones that are below the detection thresholds of existing optical surveys, are unlikely to drive the evolution of the neutral fraction around this epoch, because they would over-produce the present-day soft X-ray background. On the other hand, the seeds of the quasar BHs around z=6 likely appeared at much earlier epochs (around z=20), and produced hard ionizing radiation by accretion. These early BHs are promising candidates to account for the high redshift ionization (around z=15) implied by the recent cosmic microwave anisotropy data from WMAP. Using a model for the growth of BHs by accretion and mergers in a hierarchical cosmology, we suggest that the early growth of quasars must include a super-Eddington growth phase, and that, although not yet optically identified, the FIRST radio survey may have already detected several thousand BHs at z>6 with Mbh>10^8 Msun.

A bright quasar residing in a dense and largely neutral intergalactic medium (IGM) at high redshifts (z > 6) will be surrounded by a large cosmological Stromgren sphere. The quasar's spectrum will then show a sharp increase in resonant Lyman line absorption at wavelengths approaching and shorter than that corresponding to the Stromgren sphere's boundary along the line of sight. We show here that simultaneously considering the measured absorption in two or more hydrogen Lyman lines can provide the dynamical range required to detect this feature. We model broad and robust features of the Lyman alpha and Lyman beta regions of the spectrum of the z=6.28 quasar SDSS J1030+0524, using a hydrodynamical simulation. From the steep wavelength-dependence of the inferred absorption opacity, we detect the boundary of the Stromgren sphere at a proper distance of 6.0 +/- 0.2 Mpc away from the source redshift. From the spectrum alone, we also find that beyond this distance, cosmic hydrogen turns nearly neutral, with a neutral fraction of x_HI > 0.2, and that the ionizing luminosity of this quasar is in the range (5.2 +/- 2.5) times 10^{56} photons/sec. The method presented here, when applied to future quasars, can probe the complex topology of overlapping ionized regions, and can be used to study the details of the reionization process.

We quantify and discuss the footprints of neutral hydrogen in the intergalactic medium (IGM) on the spectra of high-redshift (z ~ 6) sources, using mock spectra generated from hydrodynamical simulations of the IGM. We show that it should be possible to extract relevant parameters, including the mean neutral fraction in the IGM, and the radius of the local cosmological Stromgren region, from the flux distribution in the observed spectra of distant sources. We focus on quasars, but a similar analysis is applicable to galaxies and gamma ray burst (GRB) afterglows. We explicitly include uncertainties in the spectral shape of the assumed source template near the Lyman alpha line. Our results suggest that a mean neutral hydrogen fraction, x(HI) of unity can be statistically distinguished from x(HI)<0.01, by combining the spectra of tens of bright (M = -27) quasars. Alternatively, the same distinction can be achieved using the spectra of several hundred sources that are ~100 times fainter. Furthermore, if the radius of the Stromgren sphere can be independently constrained to within ~10 percent, this distinction can be achieved using a single source. The information derived from such spectra will help in settling the current debate as to what extent the universe was reionized at redshifts near z=6.

The recently discovered z=10 galaxy (Pello et al 2004) has a strong Lyman alpha emission line that is consistent with being surprisingly symmetric, even given the relatively poor quality of its spectrum. The blue wing of a Lyman alpha line originating at high redshift should be strongly suppressed by resonant hydrogen absorption along the line of sight, an expectation borne out by the observed asymmetric shapes of the existing sample of Lyman alpha emitting sources at lower redshifts (3< z < 6.7). Absorption on the blue side of the line of the Pello et al. source could be reduced if the intergalactic medium (IGM) in the vicinity of the galaxy is highly ionized, but we show that this requires an unrealistically high ionizing emissivity. We suggest instead that the Lyman alpha emitting gas be receding relative to the surrounding gas with a velocity of >= 35km/s, a large velocity that is plausible only if the galaxy is part of a larger system (group of galaxies) with a velocity dispersion >= 35km/s. We find that with this velocity shift, the observed strength and shape of the line is still consistent with the galaxy being surrounded by its own Stromgren sphere embedded in a fully neutral IGM. More generally, we predict that at any given redshift, the bright Lyman alpha emitters with broader lines would exhibit stronger asymmetry than fainter ones. Bright galaxies with symmetric Lyman alpha lines may be signposts for groups and clusters of galaxies, within which they can acquire random velocities comparable to or larger than their linewidths.