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.

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.

Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($\zeta$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $\zeta$ to be biased by up to $5\sigma$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1\sigma$.

We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

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.

Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $\Delta^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude.

Once we have a measurement, how do we extract this information from the signal? This chapter focusses on the inference of the interesting astrophysics and cosmology once we obtain a detection of the 21-cm signal. Essentially, inference of the astrophysics can be broken down into three parts: 1. Characterisation of the observed data: The observed 21-cm signal varies spatially as well as along the line-of-sight (frequency or redshift dimension) to provide a full three dimensional movie of the intergalactic medium in the early Universe. However, we cannot perform a full pixel-by-pixel comparison between theoretical models and the observed signal. Instead, we require a variety of statistical methods to average the observational data in order to be able to better characterise and compare the behaviour of the faint signal. 2. An efficient method to model the 21-cm signal: In order to interpret the observations and understand the astrophysical processes responsible, we must be able to produce physically motivated models capable of replicating the signal. Further, these must be as computationally efficient as possible in order to be able to realistically investigate the 21-cm signal. 3. A robust probabilistic framework to extract the physics: The observed 21-cm signal is dependent on numerous physical processes, which within our models or simulations are described by many unknown parameters. Further, these contain approximations in order to deal with the requisite dynamic range. We must be able to characterise our ignorance in a meaningful way in order to be truly able to infer the astrophysical processes of the epoch of reionisation and cosmic dawn. In this chapter we will focus on each separately, discussing the current state-of-the-art in inferring astrophysical and cosmological information from the 21cm signal.

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.

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.

In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the systematics confronting 21 cm cosmology into a plan for overcoming them, enabling next-generation, mid-scale 21 cm arrays to be proposed late in the decade. Submitted for consideration by the Astro2020 Decadal Survey Program Panel for Radio, Millimeter, and Submillimeter Observations from the Ground as a Medium-Sized Project.

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.

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.

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.

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}$

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).

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 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.

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 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 study the redshift evolution of the quasar UV Luminosity Function (LF) for 0.5 < z < 6.5, by collecting the most up to date observational data and, in particular, the recently discovered population of faint AGNs. We fit the QSO LF using either a double power-law or a Schechter function, finding that both forms provide good fits to the data. We derive empirical relations for the LF parameters as a function of redshift and, based on these results, predict the quasar UV LF at z=8. From the inferred LF evolution, we compute the redshift evolution of the QSO/AGN comoving ionizing emissivity and hydrogen photoionization rate. If faint AGNs are included, the contribution of quasars to reionization increases substantially. However, their level of contribution critically depends on the detailed shape of the QSO LF, which can be constrained by efficient searches of high-z quasars. To this aim, we predict the expected (i) number of z>6 quasars detectable by ongoing and future NIR surveys (as EUCLID and WFIRST), and (ii) number counts for a single radio-recombination line observation with SKA-MID (FoV = 0.49 deg^2) as a function of the Hnalpha flux density, at 0<z<8. These surveys (even at z<6) will be fundamental to better constrain the role of quasars as reionization sources.

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 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.

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 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.

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.

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 develop a semi-analytic method for assessing the impact of the large-scale IGM temperature fluctuations expected following He${\rm\,{\scriptstyle II}}$ reionization on three-dimensional clustering measurements of the Ly$\alpha$ forest. Our methodology builds upon the existing large volume, mock Ly$\alpha$ forest survey simulations presented by Greig et al. by including a prescription for a spatially inhomogeneous ionizing background, temperature fluctuations induced by patchy He${\rm\,{\scriptstyle II}}$ photoheating and the clustering of quasars. This approach enables us to achieve a dynamic range within our semi-analytic model substantially larger than currently feasible with computationally expensive, fully numerical simulations. The results agree well with existing numerical simulations, with large-scale temperature fluctuations introducing a scale-dependent increase in the spherically averaged 3D Ly$\alpha$ forest power spectrum of up to 20-30 per cent at wavenumbers $k\sim0.02$ Mpc$^{-1}$. Although these large-scale thermal fluctuations will not substantially impact upon the recovery of the baryon acoustic oscillation scale from existing and forthcoming dark energy spectroscopic surveys, any complete forward modelling of the broad-band term in the Ly$\alpha$ correlation function will none the less require their inclusion.

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.

Large surveys for Lyman-alpha emitting (LAE) galaxies have been proposed as a new method for measuring clustering of the galaxy population at high redshift with the goal of determining cosmological parameters. However, Lyman-alpha radiative transfer effects may modify the observed clustering of LAE galaxies in a way that mimics gravitational effects, potentially reducing the precision of cosmological constraints. For example, the effect of the linear redshift-space distortion on the power spectrum of LAE galaxies is potentially degenerate with Lyman-alpha radiative transfer effects owing to the dependence of observed flux on intergalactic medium velocity gradients. In this paper, we show that the three-point function (bispectrum) can distinguish between gravitational and non-gravitational effects, and thus breaks these degeneracies, making it possible to recover cosmological parameters from LAE galaxy surveys. Constraints on the angular diameter distance and the Hubble expansion rate can also be improved by combining power spectrum and bispectrum measurements.

High redshift measurements of the baryonic acoustic oscillation scale (BAO) from large Ly-alpha forest surveys represent the next frontier of dark energy studies. As part of this effort, efficient simulations of the BAO signature from the Ly-alpha forest will be required. We construct a model for producing fast, large volume simulations of the Ly-alpha forest for this purpose. Utilising a calibrated semi-analytic approach, we are able to run very large simulations in 1 Gpc^3 volumes which fully resolve the Jeans scale in less than a day on a desktop PC using a GPU enabled version of our code. The Ly-alpha forest spectra extracted from our semi-analytical simulations are in excellent agreement with those obtained from a fully hydrodynamical reference simulation. Furthermore, we find our simulated data are in broad agreement with observational measurements of the flux probability distribution and 1D flux power spectrum. We are able to correctly recover the input BAO scale from the 3D Ly-alpha flux power spectrum measured from our simulated data, and estimate that a BOSS-like 10^4 deg^2 survey with ~15 background sources per square degree and a signal-to-noise of ~5 per pixel should achieve a measurement of the BAO scale to within ~1.4 per cent. We also use our simulations to provide simple power-law expressions for estimating the fractional error on the BAO scale on varying the signal-to-noise and the number density of background sources. The speed and flexibility of our approach is well suited for exploring parameter space and the impact of observational and astrophysical systematics on the recovery of the BAO signature from forthcoming large scale spectroscopic surveys.