The hyperfine transition of $^3$He$^+$ at 3.5cm has been thought as a probe of the high-z IGM since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high-z quasar to characterize the signal and analyze its prospects of detection. We find that the peak of the signal lies in the range 1-50 $\mu$K for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of HeII is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21cm that is very similar to the one from a large collection of galaxies, its signature in 3.5cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyze the prospects of the signal's detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest S/N ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large.

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium (IGM) transmission at the end of reionisation. We have gathered a sample of $13$ spectroscopically confirmed Lyman-break galaxies (LBGs) and $21$ Lyman-$\alpha$ emitters (LAEs) at angular separations $20'' \lesssim \theta \lesssim 10'$ ($\sim 0.1-4$ pMpc at $z\sim 6$) from the sightlines to $8$ background $z\gtrsim 6$ quasars. We report for the first time the detection of an excess of Lyman-$\alpha$ transmission spikes at $\sim 10-60$ cMpc from LAEs ($3.2\sigma$) and LBGs ($1.9\sigma$). We interpret the data with an improved model of the galaxy-Lyman-$\alpha$ transmission and two-point cross-correlations which includes the enhanced photoionisation due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)-Lyman-$\alpha$ transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionisation to $\langle f_{\rm esc} \rangle_{M_{\rm UV}<-12} = 0.14_{-0.05}^{+0.28}\,(0.23_{-0.12}^{+0.46})$. We investigate if the 2PCCF measurement can determine whether bright or faint galaxies are the dominant contributors to reionisation. Our results show that a contribution from faint galaxies ($M_{\rm UV} > -20 \, (2\sigma)$) is necessary to reproduce the observed 2PCCF and that reionisation might be driven by different sub-populations around LBGs and LAEs at $z\sim 6$.

Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$\alpha)\leq50$\AA\ and we find no correlation or anti-correlation between EW (Ly$\alpha$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$\alpha$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$\alpha$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}<t<20\text{Myr}$ in the past. The second object, \textit{Aerith B} at a distance $D=750$ pkpc from the quasar, displays a bright, broad double-peaked \lal emission line. Based on relations calibrated at $z\leq3$, the peak separation implies a low ionizing $f_{\text{esc}} \lesssim 1\%$, the most direct such constraint on a reionization-era galaxy. We fit the Ly-$\alpha$ line with an outflowing shell model, finding a completely typical central density $\text{log N}_{\text{HI}}/\text{cm}^{-2} = 19.3_{-0.2}^{+0.8}$, outflow velocity $v=16_{-11}^{+4}$ km s$^{-1}$, and gas temperature $\text{log} T/\text{K} = 3.8_{-0.7}^{+0.8}$ compared to $2<z<3$ analogue LAEs. Finally, we detect an emission line at $\lambda=8177$ \AA\ in object \textit{Aerith C} which, if it is \lal at $z=5.726$, would correspond closely with the end of the quasar's proximity zone ($\Delta z<0.02$ from the boundary) and suggests the quasar influences the IGM up to $\sim85$ cMpc away, making it the largest quasar proximity zone. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during the epoch of reionization.

We use hydrodynamics and radiative transfer simulations to study the 21~cm signal around a bright QSO at $z \sim 10$. Due to its powerful UV and X-ray radiation, the QSO quickly increases the extent of the fully ionized bubble produced by the pre-existing stellar type sources, in addition to partially ionize and heat the surrounding gas. As expected, a longer QSO lifetime, $t_{\rm QSO}$, results in a 21~cm signal in emission located at increasingly larger angular radii, $\theta$, and covering a wider range of $\theta$. Similar features can be obtained with a higher galactic emissivity efficiency, $f_{\rm UV}$, so that determining the origin of a large ionized bubble (i.e. QSO vs stars) is not straightforward. Such degeneracy could be reduced by taking advantage of the finite light traveltime effect, which is expected to affect an HII region produced by a QSO differently from one created by stellar type sources. From an observational point of view, we find that the 21 cm signal around a QSO at various $t_{\rm QSO}$ could be detected by SKA1-low with a high signal-noise ratio (S/N). As a reference, for $t_{\rm QSO} = 10\,\rm Myr$, a S/N $\sim 8$ is expected assuming that no pre-heating of the IGM has taken place due to high-$z$ energetic sources, while it can reach value above 10 in case of pre-heating. Observations of the 21~cm signal from the environment of a high-$z$ bright QSO could then be used to set constraints on its lifetime, as well as to reduce the degeneracy between $f_{\rm UV}$ and $t_{\rm QSO}$.

We examine Lyman continuum (LyC) leakage through H II regions regulated by turbulence and radiative feedback in a giant molecular cloud in the context of fully-coupled radiation hydrodynamics (RHD). The physical relations of the LyC escape with H I covering fraction, kinematics, spectral hardness, and the emergent Lyman-$\alpha$ (Ly$\alpha$) line profiles are studied using a series of RHD turbulence simulations performed with RAMSES-RT. The turbulence-regulated mechanism allows ionizing photons to leak out at early times before the onset of supernova feedback. The LyC photons escape through turbulence-generated low column density channels which are evacuated efficiently by radiative feedback via photoheating-induced shocks across the D-type ionization fronts. Ly$\alpha$ photons funnel through the photoionized channels along the paths of LyC escape, resulting in a diverse Ly$\alpha$ spectral morphology including narrow double-peaked profiles. The Ly$\alpha$ peak separation is controlled by the residual H I column density of the channels and the line asymmetry correlates with the porosity and multiphase structure of the H II region. This mechanism through the turbulent H II regions can naturally reproduce the observed Ly$\alpha$ spectral characteristics of some of LyC-leaking galaxies. This RHD turbulence-origin provides an appealing hypothesis to explain high LyC leakage from very young ($\sim3$ Myr) star-forming galaxies found in the local Universe without need of extreme galactic outflows nor supernova feedback. We discuss the implications of the turbulent H II regions on other nebular emission lines and a possible observational test with the Magellanic System and local blue compact dwarf galaxies as analogs of reionization-era systems.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

Galactic outflows are critical to our understanding of galaxy formation and evolution. However the details of the underlying feedback process remain unclear. We compare Ly$\alpha$ observations of the circumgalactic medium (CGM) of Lyman Break Galaxies (LBGs) with mock observations of their simulated CGM. We use cosmological hydrodynamical `zoom-in' simulations of an LBG which contains strong, momentum-driven galactic outflows. Simulation snapshots at $z=2.2$ and $z=2.65$ are used, corresponding to the available observational data. The simulation is post-processed with the radiative transfer code \textsc{crash} to account for the impact of ionising photons on hydrogen gas surrounding the simulated LBG. We generate mock absorption line maps for comparison with data derived from observed close galaxy-galaxy pairs. We perform calculations of Ly$\alpha$ photons scattering through the CGM with our newly developed Monte-Carlo code \textsc{slaf}, and compare to observations of diffuse Ly$\alpha$ halos around LBGs. Our fiducial galactic outflow model comes closer to reproducing currently observed characteristics of the CGM in Ly$\alpha$ than a reference inefficient feedback model used for comparison. Nevertheless, our fiducial model still struggles to reproduce the observed data of the inner CGM (at impact parameter $b<30$kpc). Our results suggest that galactic outflows affect Ly$\alpha$ absorption and emission around galaxies mostly at impact parameters $b<50$ kpc, while cold accretion flows dominate at larger distances. We discuss the implications of this result, and underline the potential constraining power of CGM observations - in emission and absorption - on galactic outflow models.

We present a new method to study the contribution of faint sources to the UV background using the 1D correlation of metal absorbers with the intergalactic medium (IGM) transmission in a Quasi Stellar Object (QSO) sightline. We take advantage of a sample of $25$ high signal-to-noise ratio QSO spectra to retrieve $150$ triply-ionised carbon (\cfour) absorbers at $4.5\lesssim z\lesssim 6.2$, of which $37$ systems whose expected H{~\small I} absorption lie in the Lyman-$\alpha$ forest. We derive improved constraints on the cosmic density of \cfour \,at $4.3< z < 6.2$ and infer from abundance-matching that \cfour \,absorbers trace $M_{\text{UV}}\lesssim -16$ galaxies. Correlation with the Lyman-$\alpha$ forest of the QSOs indicates that these objects are surrounded by a highly opaque region at $r\lesssim 5 $ cMpc/h followed by an excess of transmission at $r\gtrsim 10$ cMpc/h detected at $2.7\sigma$. This is in contrast to equivalent measurements at lower redshifts where only the opaque trough is detected. We interpret this excess as a statistical enhancement of the local photoionisation rate due to clustered faint galaxies around the \cfour \,absorbers. Using the analytical framework described in Paper I of this series, we derive a constraint on the average product of the escape fraction and the Lyman continuum photon production efficiency of the galaxy population clustered around the \cfour \,absorbers, $\log \langle f_{\text{esc}}\xi_{\text{ion}}\rangle /[{\rm erg^{-1}~Hz}] = 25.01^{+0.30}_{-0.19}$. This implies that faint galaxies beyond the reach of current facilities may have harder radiation fields and/or larger escape fractions than currently detected objects at the end of the reionisation epoch.

We present a new perturbative approach to "constrained Ly$\alpha$ radiative transfer'" (RT) through the circum- and inter-galactic medium (CGM and IGM). We constrain the HI content and kinematics of the CGM and IGM in a physically motivated model, using the galaxy-Ly$\alpha$ forest clustering data from spectroscopic galaxy surveys in QSO fields at $z\sim2-3$. This enables us to quantify the impact of the CGM/IGM on Ly$\alpha$ emission in an observationally constrained, realistic cosmological environment. Our model predicts that the CGM and IGM at these redshifts transmit $\approx80~\%$ of Ly$\alpha$ photons after having escaped from galaxies. This implies that while the inter-stellar medium primarily regulates Ly$\alpha$ escape, the CGM has a non-negligible impact on the observed Ly$\alpha$ line properties and the inferred Ly$\alpha$ escape fraction, even at $z\sim 2-3$. Ly$\alpha$ scattering in the CGM and IGM further introduces an environmental dependence in the (apparent) Ly$\alpha$ escape fraction, and the observed population of Ly$\alpha$ emitting galaxies: the CGM/IGM more strongly suppresses direct Ly$\alpha$ emission from galaxies in overdense regions in the Universe, and redistributes this emission into brighter Ly$\alpha$ haloes. The resulting mean surface brightness profile of the Ly$\alpha$ haloes is generally found to be a power-law $\propto r^{-2.4}$. Although our model still contains arbitrariness, our results demonstrate how (integral field) spectroscopic surveys of galaxies in QSO fields constrain circumgalactic Ly$\alpha$ RT, and we discuss the potential of these models for studying CGM physics and cosmology.

We introduce a new method for determining the influence of galaxies and active galactic nuclei (AGN) on the physical state of the intergalactic medium (IGM) at high redshift and illustrate its potential via a first application to the field of the $z=6.42$ QSO J1148+5251. By correlating the spatial positions of spectroscopically-confirmed Lyman break galaxies (LBGs) with fluctuations in the Lyman alpha forest seen in the high signal-to-noise spectrum of a background QSO, we provide a statistical measure of the typical escape fraction of Lyman continuum photons close to the end of cosmic reionisation. Here we use Keck DEIMOS spectroscopy to locate 7 colour-selected LBGs in the redshift range $5.3\lesssim z\lesssim 6.4$ and confirm a faint $z=5.701$ AGN. We then examine the spatial correlation between this sample and Ly$\alpha$/Ly$\beta$ transmission fluctuations in a Keck ESI spectrum of the QSO. Interpreting the statistical HI proximity effect as arising from faint galaxies clustered around the detected LBGs, we translate the observed mean Ly$\alpha$ transmitted flux around an average detected LBG into a constraint on the mean escape fraction $\langle f_{\rm esc}\rangle\geq0.08$ at $z\simeq6$. We also report evidence of the individual transverse HI proximity effect of a $z=6.177$ luminous LBG via a Ly$\beta$ transmission spike and two broad Ly$\alpha$ transmission spikes around the $z=5.701$ AGN. We discuss the possible origin of such associations which suggest that while faint galaxies are primarily driving reionisation, luminous galaxies and AGN may provide important contributions to the UV background or thermal fluctuations of the IGM at $z\simeq6$. Although a limited sample, our results demonstrate the potential of making progress using this method in resolving one of the most challenging aspects of the contribution of galaxies and AGN to cosmic reionisation.

Observations of the 21 cm line from neutral hydrogen indicate that an Epoch of Heating (EoH) might have preceded the later Epoch of Reionization (EoR). Here we study the effects on the ionization state and the thermal history of the Intergalactic Medium (IGM) during the EoH induced by different assumptions on ionizing sources in the high redshift Universe: (i) stars, (ii) X-ray binaries (XRBs), (iii) thermal bremsstrahlung of the hot Interstellar Medium (ISM), and (iv) accreting nuclear black holes (BHs). To this aim, we post-process outputs from the ($100 h^{-1}$ cMpc)$^3$ hydrodynamical simulation MassiveBlack-II with the cosmological 3D radiative transfer code CRASH, which follows the propagation of UV and X-ray photons, computing the thermal and ionization state of hydrogen and helium through the EoH. We find that stars determine the fully ionized morphology of the IGM, while the spectrally hard XRBs pave way for efficient subsequent heating and ionization by the spectrally softer ISM. With the seeding prescription in MassiveBlack-II, BHs do not contribute significantly to either ionization or heating. With only stars, most of the IGM remains in a cold state (with a median $T=11$ K at $z=10$), however, the presence of more energetic sources raises the temperature of regions around the brightest and more clustered sources above that of the CMB, opening the possibility to observing the 21 cm signal in emission.

We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km).

We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $\Delta^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$\sigma$) excess variance and a best 2-$\sigma$ upper limit of $\Delta^2_{\rm 21} < (79.6 {\rm mK})^2$ at $k=0.053$$h$cMpc$^{-1}$ in the range $z=$9.6-10.6. The excess variance decreases when optimizing the smoothness of the direction- and frequency-dependent gain calibration, and with increasing the completeness of the sky model. It is likely caused by (i) residual side-lobe noise on calibration baselines, (ii) leverage due to non-linear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications.

We present spatial correlations of galaxies and IGM HI in the COSMOS/UltraVISTA 1.62 deg$^2$ field. Our data consist of 13,415 photo-$z$ galaxies at $z\sim2-3$ with $K_s<23.4$ and the Ly$\alpha$ forest absorptions in the background quasar spectra selected from SDSS data with no signature of damped Ly$\alpha$ system contamination. We estimate a galaxy overdensity $\delta_{gal}$ in an impact parameter of 2.5 pMpc, and calculate the Ly$\alpha$ forest fluctuations $\delta_{\langle F\rangle}$ whose negative values correspond to the strong Ly$\alpha$ forest absorptions. We identify weak evidence of an anti-correlation between $\delta_{gal}$ and $\delta_{\langle F\rangle}$ with a Spearman's rank correlation coefficient of $-0.39$ suggesting that the galaxy overdensities and the Ly$\alpha$ forest absorptions positively correlate in space at the $\sim90\%$ confidence level. This positive correlation indicates that high-$z$ galaxies exist around an excess of HI gas in the Ly$\alpha$ forest. We find four cosmic volumes, dubbed $A_{obs}$-$D_{obs}$, that have extremely large (small) values of $\delta_{gal} \simeq0.8$ ($-1$) and $\delta_{\langle F\rangle}$ $\simeq0.1$ ($-0.4$), three out of which, $B_{obs}$-$D_{obs}$, significantly depart from the correlation, and weaken the correlation signal. We perform cosmological hydrodynamical simulations, and compare with our observational results. Our simulations reproduce the correlation, agreeing with the observational results. Moreover, our simulations have model counterparts of $A_{obs}$-$D_{obs}$, and suggest that the observations pinpoint, by chance, a galaxy overdensity like a proto-cluster, gas filaments lying on the sightline, a large void, and orthogonal low-density filaments. Our simulations indicate that the significant departures of $B_{obs}$-$D_{obs}$ are produced by the filamentary large-scale structures and the observation sightline effects.

We present a detailed analysis of the ionization and thermal structure of the intergalactic medium (IGM) around a high-redshift QSO using a large suite of cosmological, multi-frequency radiative transfer (RT) simulations, exploring the contribution from galaxies as well as the QSO, and the effect of X-rays and secondary ionization. We show that in high-z QSO environments both the central QSO and the surrounding galaxies concertedly control the reionization morphology of hydrogen and helium and have a non-linear impact on the thermal structure of the IGM. A QSO imprints a distinctive morphology on H II regions if its total ionizing photon budget exceeds that of the surrounding galaxies since the onset of hydrogen reionization; otherwise, the morphology shows little difference from that of H II regions produced only by galaxies. In addition, the spectral shape of the collective radiation field from galaxies and QSOs controls the thickness of the I-fronts. While a UV-obscured QSO can broaden the I-front, the contribution from other UV sources, either galaxies or unobscured QSO, is sufficient to maintain a sharp I-front. X-rays photons from the QSO are responsible for a prominent extended tail of partial ionization ahead of the I-front. QSOs leave a unique imprint on the morphology of He II / He III regions. We suggest that, while the physical state of the IGM is modified by QSOs, the most direct test to understand the role of galaxies and QSOs during reionization is to perform galaxy surveys in a region of sky imaged by 21 cm tomography.

The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged]

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.

The first major star-forming galaxies and Active Galactic Nuclei will produce Balmer and higher order extended haloes during the Epoch of Reionization through the scattering of Lyman resonance line photons off the surrounding neutral intergalactic gas. The optical depth dependence of the scattering rates will produce a signal sensitive to both the density and velocity fluctuations of the gas, offering the possibility of probing the ionization region and flow field surrounding young star-forming galaxies. The requirements for detecting the haloes in the infra-red using a space-based telescope are discussed, along with an assessment of the possibility of detecting the haloes using the Tunable Filter Imager on the James Webb Space Telescope.