Search

We predict the 21-cm global signal and power spectra during the Epoch of Reionisation using the MERAXES semi-analytic galaxy formation and reionisation model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionisation of cosmic hydrogen using semi-analytic models (such as MERAXES) requires N-body simulations within large volumes and high mass resolutions. For this, we use a simulation of side-length $210~h^{-1}$ Mpc with $4320^3$ particles resolving dark matter haloes to masses of $5\times10^8~h^{-1}~M_\odot$. To reach the mass resolution of atomically cooled galaxies, thought to be the dominant population contributing to reionisation, at $z=20$ of $\sim 2\times10^7~h^{-1}~M_\odot$, we augment this simulation using the DARKFOREST Monte-Carlo merger tree algorithm (achieving an effective particle count of $\sim10^{12}$). Using this augmented simulation we explore the impact of mass resolution on the predicted reionisation history as well as the impact of X-ray heating on the 21-cm global signal and the 21-cm power spectra. We also explore the cosmic variance of 21-cm statistics within $70^{3}$ $h^{-3}$ Mpc$^3$ sub-volumes. We find that the midpoint of reionisation varies by $\Delta z\sim0.8$ and that the cosmic variance on the power spectrum is underestimated by a factor of $2-4$ at $k\sim 0.1-0.4$ Mpc$^{-1}$ due to the non-Gaussian nature of the 21-cm signal. To our knowledge, this work represents the first model of both reionisation and galaxy formation which resolves low-mass atomically cooled galaxies while simultaneously sampling sufficiently large scales necessary for exploring the effects of X-rays in the early Universe., Comment: 17 pages, 13 figures, revised and MNRAS accepted version

Modeling the 21cm global signal from the Cosmic Dawn is challenging due to the many poorly constrained physical processes that come into play. We address this problem using the semi-analytical code "Cosmic Archaeology Tool" (CAT). CAT follows the evolution of dark matter halos tracking their merger history and provides an ab initio description of their baryonic evolution, starting from the formation of the first (Pop III) stars and black holes (BHs) in mini-halos at z > 20. The model is anchored to observations of galaxies and AGN at z < 6 and predicts a reionization history consistent with constraints. In this work we compute the evolution of the mean global 21cm signal between $4\leq z \leq 40$ based on the rate of formation and emission properties of stars and accreting black holes. We obtain an absorption profile with a maximum depth $δ{\rm T_b} = -95$ mK at $z \sim 26.5$ (54 MHz). This feature is quickly suppressed turning into an emission signal at $z = 20$ due to the contribution of accreting BHs that efficiently heat the IGM at $z < 27$. The high-$z$ absorption feature is caused by the early coupling between the spin and kinetic temperature of the IGM induced by Pop III star formation episodes in mini-halos. Once we account for an additional radio background from early BHs, we are able to reproduce the timing and the depth of the EDGES signal only if we consider a smaller X-ray background from accreting BHs, but not the shape., 16 pages, 10 figures, Accepted for publication in MNRAS

Correlations between black holes and their host galaxies provide insight into what drives black hole-host co-evolution. We use the Meraxes semi-analytic model to investigate the growth of black holes and their host galaxies from high redshift to the present day. Our modelling finds no significant evolution in the black hole-bulge and black hole-total stellar mass relations out to a redshift of 8. The black hole-total stellar mass relation has similar but slightly larger scatter than the black hole-bulge relation, with the scatter in both decreasing with increasing redshift. In our modelling the growth of galaxies, bulges and black holes are all tightly related, even at the highest redshifts. We find that black hole growth is dominated by instability-driven or secular quasar-mode growth and not by merger-driven growth at all redshifts. Our model also predicts that disc-dominated galaxies lie on the black hole-total stellar mass relation, but lie offset from the black hole-bulge mass relation, in agreement with recent observations and hydrodynamical simulations., Comment: 14 pages, 11 figures. Accepted for publication in MNRAS

The bright end of the rest-frame UV luminosity function (UVLF) of high-redshift galaxies is modified by gravitational lensing magnification bias. Motivated by recent discoveries of very high-z galaxies with JWST, we study the dependence of magnification bias on the finite size of sources at $6, Comment: 5 pages, 4 figures, 1 table. Accepted for publication in MNRAS Letters

The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply Markov Chain Monte Carlo-based PSF modelling to determine the detectability of their host galaxies with the James Webb Space Telescope (JWST). While no statistically significant detections are made with HST, we predict that at the same wavelengths and exposure times JWST NIRCam imaging will detect ~50% of quasar host galaxies. We investigate various observational strategies, and find that NIRCam wide-band imaging in the long-wavelength filters results in the highest fraction of successful quasar host detections, detecting >80% of the hosts of bright quasars in exposure times of 5 ks. Exposure times of ~5 ks are required to detect the majority of host galaxies in the NIRCam wide-band filters, however even 10 ks exposures with MIRI result in, 21 pages, 14 figures. Resubmitted to MNRAS following referee feedback, with most notable improvements in Section 5

The highly neutral inter-galactic medium (IGM) during the Epoch of Reionization (EoR) is expected to suppress Ly$\alpha$ emission with damping-wing absorption, causing nearly no Ly$\alpha$ detection from star-forming galaxies at $z{\sim}8$. However, spectroscopic observations of the 4 brightest galaxies (${\rm H}_{160}{\sim}25$ mag) at these redshifts do reveal prominent Ly$\alpha$ line, suggesting locally ionised IGM. In this paper, we explore the Ly$\alpha$ IGM transmission and environment of bright galaxies during the EoR using the Meraxes semi-analytic model. We find brighter galaxies to be less affected by damping-wing absorption as they are effective at ionizing surrounding neutral hydrogen. Specifically, the brightest sources (${\rm H}_{160}{\lesssim}25.5$ mag) lie in the largest ionized regions in our simulation, and have low attenuation of their Ly$\alpha$ from the IGM (optical depth ${, Comment: 10 pages, 3 figures, submitted to MNRAS, comments are welcome

Resolving faint galaxies in large volumes is critical for accurate cosmic reionisation simulations. While less demanding than hydrodynamical simulations, semi-analytic reionisation models still require very large N-body simulations in order to resolve the atomic cooling limit across the whole reionisation history within box sizes $\gtrsim 100 \, h^{-1} {\rm Mpc}$. To facilitate this, we extend the mass resolution of N-body simulations using a Monte Carlo algorithm. We also propose a method to evolve positions of Monte Carlo halos, which can be an input for semi-analytic reionisation models. To illustrate, we present an extended halo catalogue that reaches a mass resolution of $M_\text{halo} = 3.2 \times 10^7 \, h^{-1} \text{M}_\odot$ in a $105 \, h^{-1} {\rm Mpc}$ box, equivalent to an N-body simulation with $\sim 6800^3$ particles. The resulting halo mass function agrees with smaller volume N-body simulations with higher resolution. Our results also produce consistent two-point correlation functions with analytic halo bias predictions. The extended halo catalogues are applied to the \textsc{meraxes} semi-analytic reionisation model, which improves the predictions on stellar mass functions, star formation rate densities and volume-weighted neutral fractions. Comparison of high resolution large volume simulations with both small volume or low resolution simulations confirms that both low resolution and small volume simulations lead to reionisation ending too rapidly. Lingering discrepancies between the star formation rate functions predicted with and without our extensions can be traced to the uncertain contribution of satellite galaxies., 15 pages, 15 figures. Accepted for publication in MNRAS

arXiv:2007.11600v1, We present a first strong-lensing model for the galaxy cluster RM J121218.5+273255.1 (z = 0.35; hereafter RMJ1212; also known as Abell 1489). This cluster is amongst the top 0.1% richest clusters in the redMaPPer catalog; it is significantly detected in X-ray and through the Sunyaev-Zel’dovich effect in ROSAT and Planck data, respectively; and its optical luminosity distribution implies a very large lens, following mass-to-light scaling relations. Based on these properties it was chosen for the Webb Medium Deep Fields (WMDF) JWST/GTO program. In preparation for this program, RMJ1212 was recently imaged with GMOS on Gemini North and in seven optical and near-infrared bands with the Hubble Space Telescope. We use these data to map the inner mass distribution of the cluster, uncovering various sets of multiple images. We also search for high-redshift candidates in the data, as well as for transient sources. We find over a dozen high-redshift (z & 6) candidates based on both photometric redshift and the dropout technique. No prominent (& 5σ) transients were found in the data between the two HST visits. Our lensing analysis reveals a relatively large lens with an effective Einstein radius of θE ' 32 ± 3 00 (zs = 2), in broad agreement with the scaling-relation expectations. RMJ1212 demonstrates that powerful lensing clusters can be selected in a robust and automated way following the light-traces-mass assumption., M.N. acknowledges INAF 1.05.01.86.20 and PRIN MIUR F.OB. 1.05.01.83.08. J.M.D. acknowledges the support of project PGC2018-101814-B-100 (MCIU/AEI/MINECO/FEDER, UE) Ministerio de Ciencia, Investigacion y Universidades. This project was funded by the Agencia Estatal de Investigacion, Unidad de Excelencia Mar´ıa de Maeztu, ref. MDM2017-0765.

A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21cm line. Using the \textsc{BlueTides} simulation, we investigate the measurement of an \textit{average} ionised region towards the beginning of reionization by stacking redshifted 21cm images around optically identified bright galaxies using mock observations. We find that with an SKA 1000 hour observation, assuming perfect foreground subtraction, a $5\sigma$ detection of a stacked HII region can be made with 30 images around some of the brightest galaxies in \textsc{bluetides} (brighter than $M_{UV} < -22.75$) at $z=9$ (corresponding to a neutral fraction of 90.1 \% in our model). We present simulated relationships between the UV magnitude of galaxies, the sizes of the ionised regions they reside in, and the shape of the stacked profiles. These mock observations can also distinguish between scenarios where the IGM is in net emission or absorption of 21cm photons. Once 21cm foreground contamination is included, we find that even with up to 200 images around these rare, bright galaxies, only a tentative $> 1\sigma$ detection will be possible. However, partial foreground subtraction substantially improves signal-to-noise. For example, we predict that reducing the area of Fourier space dominated by foregrounds by 50 (80) percent will allow $> 3\sigma$ ($> 5\sigma$) detections of ionised regions at $z=9$., Comment: 12 pages, 17 figures, Accepted for publication in MNRAS

We use the latest multi-redshift ($z=6.5-8.7$) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. These upper limits are achieved using 298 h of carefully excised data over four observing seasons. To explore these upper limits in the context of reionisation astrophysics, we use 21CMMC. We then connect the disfavoured regions of parameter space to existing observational constraints on reionisation such as high-$z$ galaxy ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of models disfavoured by the MWA limits are already inconsistent with existing observational constraints. These inconsistent models arise from two classes of models: (i) `cold' reionisation and (ii) pure matter density fluctuations (i.e. no reionisation). However, a small subsample of models are consistent implying the existing MWA limits provide unique information in disfavouring models of reionisation, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding $1\sigma$ lower limits of $\epsilon_{{\rm X},0.5-2~{\rm keV}}\gtrsim10^{34.5}$ erg s$^{-1}$ Mpc$^{-3}$. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, 2.4 K at $z=6.5, 6.8, 7.1, 7.8, 8.2, 8.7$. With this we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on $\epsilon_{{\rm X},0.5-2~{\rm keV}}$ and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction., Comment: 14 pages, 7 figures, 1 table. Submitted to MNRAS, comments welcome. arXiv admin note: text overlap with arXiv:2006.03203

We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the 10 most massive black holes and the 191 quasars in the simulation (with $M_{\textrm{UV,AGN}}10^{6.5}M_\odot$ black hole, indicating that minor mergers play at least some role in growing black holes in the early Universe., Comment: 19 pages, 20 figures. Accepted for publication in MNRAS

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, 11 pages, 7 figures; Updated to match the published version. All results and conclusions remain unchanged

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., 17 pages, 27 figures, accepted for publication in MNRAS, minor updates to align with final published version

We present a new analysis of high-redshift UV observations using a semi-analytic galaxy formation model, and provide self-consistent predictions of the infrared excess (IRX) -- $\beta$ relations and cosmic star formation rate density. We combine the Charlot & Fall dust attenuation model with the Meraxes semi-analytic model, and explore three different parametrisations for the dust optical depths, linked to star formation rate, dust-to-gas ratio and gas column density respectively. A Bayesian approach is employed to statistically calibrate model free parameters including star formation efficiency, mass loading factor, dust optical depths and reddening slope directly against UV luminosity functions and colour-magnitude relations at z ~ 4-7. The best-fit models show excellent agreement with the observations. We calculate IRX using energy balance arguments, and find that the large intrinsic scatter in the IRX -$\beta$ plane is driven by the specific star formation rate. Additionally, the difference among the three dust models suggests a factor of two systematic uncertainty in the dust-corrected star formation rate when using the Meurer IRX - $\beta$ relation at z > 4., Comment: 15 pages, 10 figure. Accepted for publication in MNRAS

We study the sizes, angular momenta and morphologies of high-redshift galaxies using an update of the Meraxes semi-analytic galaxy evolution model. Our model successfully reproduces a range of observations from redshifts $z=0$-$10$. We find that the effective radius of a galaxy disc scales with UV luminosity as $R_e\propto L_{\textrm{UV}}^{0.33}$ at $z=5$-$10$, and with stellar mass as $R_e\propto M_\ast^{0.24}$ at $z=5$ but with a slope that increases at higher redshifts. Our model predicts that the median galaxy size scales with redshift as $R_e \propto (1+z)^{-m}$, where $m=1.98\pm0.07$ for galaxies with $(0.3$-$1)L^\ast_{z=3}$ and $m=2.15\pm0.05$ for galaxies with $(0.12$-$0.3)L^\ast_{z=3}$. We find that the ratio between stellar and halo specific angular momentum is typically less than one and decreases with halo and stellar mass. This relation shows no redshift dependence, while the relation between specific angular momentum and stellar mass decreases by $\sim0.5$ dex from $z=7$ to $z=2$. Our model reproduces the distribution of local galaxy morphologies, with bulges formed predominantly through galaxy mergers for low-mass galaxies, disc-instabilities for galaxies with $M_\ast\simeq10^{10}$-$10^{11.5}M_\odot$, and major mergers for the most massive galaxies. At high redshifts, we find galaxy morphologies that are predominantly bulge-dominated., 21 pages, 12 figures. Accepted for publication in MNRAS

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 1e910$), 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., Comment: 18 pages, 9 figures; submitted to MNRAS, comments are welcome

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., Comment: 17 pages, 17 figures, Accepted for publication in MNRAS

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{}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 present the results of a VLT MUSE/FORS2 and Spitzer survey of a unique compact lensing cluster CLIO at z = 0.42, discovered through the GAMA survey using spectroscopic redshifts. Compact and massive clusters such as this are understudied, but provide a unique prospective on dark matter distributions and for finding background lensed high-z galaxies. The CLIO cluster was identified for follow up observations due to its almost unique combination of high mass and dark matter halo concentration, as well as having observed lensing arcs from ground based images. Using dual band optical and infra-red imaging from FORS2 and Spitzer, in combination with MUSE optical spectroscopy we identify 89 cluster members and find background sources out to z = 6.49. We describe the physical state of this cluster, finding a strong correlation between environment and galaxy spectral type. Under the assumption of a NFW profile, we measure the total mass of CLIO to be M$_{200} = (4.49 \pm 0.25) \times 10^{14}$ M$_\odot$. We build and present an initial strong-lensing model for this cluster, and measure a relatively low intracluster light (ICL) fraction of 7.21 $\pm$ 1.53% through galaxy profile fitting. Due to its strong potential for lensing background galaxies and its low ICL, the CLIO cluster will be a target for our 110 hour JWST 'Webb Medium-Deep Field' (WMDF) GTO program., Accepted at MNRAS, 19 pages, 16 figures

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering., Comment: 10 pages, 6 figures, 2 tables, accepted for publication in MNRAS

We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and increases again for higher-mass haloes. However, the galaxy gas accretion does not flatten at intermediate halo masses when AGN feedback is switched off. To better understand these trends, we develop a physically motivated semi-analytic model of galaxy gas accretion. We show that the flattening is produced by the rate of gas cooling from the hot halo. The ratio of the cooling radius and the virial radius does not decrease continuously with increasing halo mass as generally thought. While it decreases up to ${\sim}10^{13}M_{\odot}$ haloes, it increases for higher halo masses, causing an upturn in the galaxy gas accretion rate. This may indicate that in high-mass haloes AGN feedback is not sufficiently efficient. When there is no AGN feedback, the density of the hot halo is higher, the ratio of the cooling and virial radii does not decrease as much and the cooling rate is higher. Changes in the efficiency of stellar feedback can also increase or decrease the accretion rates onto galaxies. The trends can plausibly be explained by the re-accretion of gas ejected by progenitor galaxies and by the suppression of black hole growth, and hence AGN feedback, by stellar feedback., Comment: 17 pages, 9 figures. Accepted for publication in MNRAS. Comments welcome

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., 14 pages, 10 figures, published in MNRAS

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, 21 pages, 12 figures; Updated to match the published version. All results and conclusions remain unchanged

The formation and evolution of galaxies with low neutral atomic hydrogen (HI) masses, M$_{\rm HI}$$, Comment: 14 pages, 10 figures, accepted for pubication in MNRAS

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., Comment: 14 pages, 11 figures, 2 tables, published in MNRAS. Updated to match published version. Additional results and comments added from previous version. All other results and conclusions remain unchanged

We use high-resolution simulations of cosmological volumes to model galaxy formation at high-redshift, with the goal of studying the photon budget for reionization. We demonstrate that galaxy formation models that include a strong, thermally coupled supernovae scheme reproduce current observations of star formation rates and specific star formation rates, both during and after the reionization era. These models produce enough UV photons to sustain reionization at z, 13 pages, 10 figures, 1 table. MNRAS accepted, changed to reflect published version

We show that a bursty model of high redshift star formation explains several puzzling observations of the high redshift galaxy population. We begin by pointing out that the observed specific star formation rate requires a duty-cycle of ~10%, which is much lower than found in many hydro-dynamical simulations. This value follows directly from the fact that the observed star formation rate in galaxies integrated over a Hubble time exceeds the observed stellar mass by an order of magnitude. We use the large observed specific star formation rate to calibrate the efficiency of feedback in a model for the high redshift star formation rate which includes merger driven star formation regulated by SNe feedback. This model reproduces the star formation rate density function and the stellar mass function of galaxies at 4, Comment: 12 pages, 8 figures. Submitted to MNRAS. arXiv admin note: substantial text overlap with arXiv:1209.2215

Large surveys for Lyα 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, Lyα 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. We investigate the impact of Lyα radiative transfer effects on the observed clustering of LAE galaxies. In particular, we focus on the effects of the intergalactic medium velocity gradients, local density within the environment of an LAE galaxy and ionizing background fluctuations. For example, the effect of the linear redshift-space distortion on the power spectrum of LAE galaxies is potentially degenerate with Lyα 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 Hubble expansion rate are independent of Lyα radiative transfer degeneracies; however, they incur slight reductions in their constraining power resulting from the overall reduction of the signal-to-noise due to the Lyα radiative transfer effects. Combining the power spectrum and bispectrum measurements provides improved constraints on the angular diameter distance and Hubble expansion rate. © 2013 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

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., 10 pages, 4 figures, Accepted for publication in MNRAS

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., Comment: 14 pages, 12 figures, published in MNRAS

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., 30 pages, 18 figures; published in MNRAS; Updated to match published version. Additional results and fitting formulae added from previous version. All other results and conclusions remain unchanged

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

The 21 cm forest – H i absorption features in the spectra of high-redshift radio sources – can potentially provide a unique probe of the largely neutral intergalactic medium (IGM) during the epoch of reionization. We present semi-analytical models of the 21 cm forest due to the large-scale structure of the reionization-era IGM, including a prescription for X-ray heating and the percolation of photoionization bubbles. We explore a range of signal-analysis methods to show that, if detected with future instruments such as the Square Kilometre Array (SKA), the 21 cm forest can provide a significant constraint on the thermal history of the IGM. Detection will be aided by consideration of the sudden increase in signal variance at the onset of 21 cm absorption. If radio foregrounds and the intrinsic source spectra are well understood, the flux decrement over wide bandwidths can also improve detection prospects. Our analysis accounts for the possibility of narrow absorption lines from intervening dense regions, but, unlike previous studies, our results do not depend on their properties. Assuming X-ray heating corresponding to a local stellar population and a simple reionization model, we estimate that a statistically significant detection of 21 cm absorption could be made by the SKA in less than a year of observing against a Cygnus A-type source at z ∼ 9, as opposed to nearly a decade for a significant detection of the detailed forest features. We discuss observational challenges due to uncertainties regarding the abundance of background sources and the strength of the 21 cm absorption signal.

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

Measurements of the intergalactic medium (IGM) temperature provide a potentially powerful constraint on the reionisation history due to the thermal imprint left by the photo-ionisation of neutral hydrogen. However, until recently IGM temperature measurements were limited to redshifts 2 6. In this work, we use recent measurements of the IGM temperature in the near-zones of seven quasars at z ~ 5.8 - 6.4, combined with a semi-numerical model for inhomogeneous reionisation, to establish new constraints on the redshift at which hydrogen reionisation completed. We calibrate the model to reproduce observational constraints on the electron scattering optical depth and the HI photo-ionisation rate, and compute the resulting spatially inhomogeneous temperature distribution at z ~ 6 for a variety of reionisation scenarios. Under standard assumptions for the ionising spectra of population-II sources, the near-zone temperature measurements constrain the redshift by which hydrogen reionisation was complete to be z > 7.9 (6.5) at 68 (95) per cent confidence. We conclude that future temperature measurements around other high redshift quasars will significantly increase the power of this technique, enabling these results to be tightened and generalised.

Narrowband observations have detected no Ly-alpha emission within a 70 pMpc^3 volume centered on the z = 2.168 quasar PKS 0424-131. This is in contrast to surveys of Ly-alpha emitters in the field at similar redshifts and flux limits, which indicate that tens of sources should be visible within the same volume. The observed difference indicates that the quasar environment has a significant influence on the observed density of Ly-alpha emitters. To quantify this effect we have constructed a semi-analytic model to simulate the effect of a luminous quasar on nearby Ly-alpha emitters. We find the null detection around PKS 0424-131 implies that the minimum isothermal temperature of Ly-alpha emitter host halos is greater than 3.4 x 10^6 K (68% level), corresponding to a virial mass of ~1.2 x 10^12 M_solar. This indicates that the intense UV emission of the quasar may be suppressing the star formation in nearby galaxies. Our study illustrates that low redshift quasar environments may serve as a surrogate for studying the radiative suppression of galaxy formation during the epoch of reionization.

We present measurements of the intergalactic medium (IGM) temperature within ~5 proper Mpc of seven luminous quasars at z~6. The constraints are obtained from the Doppler widths of Lyman-alpha absorption lines in the quasar near-zones and build upon our previous measurement for the z=6.02 quasar SDSS J0818+1722. The expanded data set, combined with an improved treatment of systematic uncertainties, yields an average temperature at the mean density of log (T_0/K) = 4.21 \pm 0.03 (\pm^0.06_0.07) at 68 (95) per cent confidence for a flat prior distribution over 3.2 < log (T_0/K) < 4.8. In comparison, temperatures measured from the general IGM at z~5 are ~0.3 dex cooler, implying an additional source of heating around these quasars which is not yet present in the general IGM at slightly lower redshift. This heating is most likely due to the recent reionisation of He-II in vicinity of these quasars, which have hard, non-thermal ionising spectra. The elevated temperatures may therefore represent evidence for the earliest stages of He-II reionisation in the most biased regions of the high-redshift Universe. The temperature as a function of distance from the quasars is consistent with being constant, log(T_0/K)~4.2, with no evidence for a line-of-sight thermal proximity effect. However, the limited extent of the quasar near-zones prevents the detection of He-III regions larger than ~5 proper Mpc. Under the assumption the quasars have reionised the He-II in their vicinity, we infer the data are consistent with an average optically bright phase of duration in excess of 10^6.5 yr. These measurements represent the highest-redshift IGM temperature constraints to date, and thus provide a valuable data set for confronting models of H-I reionisation.

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.

Measurement of redshifted 21-cm emission from neutral hydrogen promises to be the most effective method for studying the reionization history of hydrogen and, indirectly, the first galaxies. These studies will be limited not by raw sensitivity to the signal, but rather, by bright foreground radiation from Galactic and extragalactic radio sources and the Galactic continuum. In addition, leakage due to gain errors and non-ideal feeds conspire to further contaminate low-frequency radio observations. This leakage leads to a portion of the complex linear polarization signal finding its way into Stokes I, and inhibits the detection of the non-polarized cosmological signal from the epoch of reionization. In this work, we show that rotation measure synthesis can be used to recover the signature of cosmic hydrogen reionization in the presence of contamination by polarized foregrounds. To achieve this, we apply the rotation measure synthesis technique to the Stokes I component of a synthetic data cube containing Galactic foreground emission, the effect of instrumental polarization leakage and redshifted 21-cm emission by neutral hydrogen from the epoch of reionization. This produces an effective Stokes I Faraday dispersion function for each line of sight, from which instrumental polarization leakage can be fitted and subtracted. Our results show that it is possible to recover the signature of reionization in its late stages (z≈ 7) by way of the 21-cm power spectrum, as well as through tomographic imaging of ionized cavities in the intergalactic medium. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

The measurement of the light scattered from extrasolar planets informs atmospheric and formation models. With the discovery of many hot Jupiter planets orbiting nearby stars, this motivates the development of robust methods of characterisation from follow up observations. In this paper we discuss two methods for determining the planetary albedo in transiting systems. First, the most widely used method for measuring the light scattered by hot Jupiters (Collier Cameron et al.) is investigated for application for typical echelle spectra of a transiting planet system, showing that detection requires high signal-to-noise ratio data of bright planets. Secondly a new Fourier analysis method is also presented, which is model-independent and utilises the benefits of the reduced number of unknown parameters in transiting systems. This approach involves solving for the planet and stellar spectra in Fourier space by least-squares. The sensitivities of the methods are determined via Monte Carlo simulations for a range of planet-to-star fluxes. We find the Fourier analysis method to be better suited to the ideal case of typical observations of a well constrained transiting system than the Collier Cameron et al. method. We apply the Fourier analysis method for extracting the light scattered by transiting hot Jupiters from high resolution spectra to echelle spectra of HD 209458 and HD 189733. Unfortunately we are unable to improve on the previous upper limit of the planet-to-star flux for HD 209458b set by space-based observations. A 1{\sigma}upper limit on the planet-to-star flux of HD 189733b is measured in the wavelength range of 558.83-599.56 nm yielding {\epsilon} < 4.5 \times 10-4. Improvement in the measurement of the upper limit of the planet-to-star flux of this system, with ground-based capabilities, requires data with a higher signal-to-noise ratio, and increased stability of the telescope.

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

The observed number counts of high-redshift galaxy candidates have been used to build up a statistical description of star-forming activity at redshift z >~ 7, when galaxies reionized the Universe. Standard models predict that a high incidence of gravitational lensing will probably distort measurements of flux and number of these earliest galaxies. The raw probability of this happening has been estimated to be ~ 0.5 percent, but can be larger owing to observational biases. Here we report that gravitational lensing is likely to dominate the observed properties of galaxies with redshifts of z >~ 12, when the instrumental limiting magnitude is expected to be brighter than the characteristic magnitude of the galaxy sample. The number counts could be modified by an order of magnitude, with most galaxies being part of multiply imaged systems, located less than 1 arcsec from brighter foreground galaxies at z ~ 2. This lens-induced association of high-redshift and foreground galaxies has perhaps already been observed among a sample of galaxy candidates identified at z ~ 10.6. Future surveys will need to be designed to account for a significant gravitational lensing bias in high-redshift galaxy samples., Nature, Jan. 13, 2011 issue (in press)

In this paper, we explore the constraints on the reionization history that are provided by current observations of the Lyman-alpha forest and the CMB. Rather than using a particular semi-analytic model, we take the novel approach of parametrizing the ionizing sources with arbitrary functions, and perform likelihood analyses to constrain possible reionization histories. We find model independent conclusions that reionization is likely to be mostly complete by z=8 and that the IGM was 50% ionized at z=9-10. Upcoming low-frequency observations of the redshifted 21 cm line of neutral hydrogen are expected to place significantly better constraints on the hydrogen neutral fraction at 6

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

Measurement of the spatial distribution of neutral hydrogen via the redshifted 21 cm line promises to revolutionize our knowledge of the epoch of reionization and the first galaxies, and may provide a powerful new tool for observational cosmology from redshifts 1, Comment: Invited review for Annual Review of Astronomy and Astrophysics (2010 volume)

The thermal state of the intergalactic medium (IGM) provides an indirect probe of both the HI and HeII reionisation epochs. Current constraints on the IGM temperature from the Lya forest are restricted to the redshift range 2