A new upper limit on the 21-cm signal power spectrum at a redshift of z ≈ 9.1 is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen 'excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of 0.25 − 0.45 MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-σ upper limit of ∆ 2 21 < (73) 2 mK 2 at k = 0.075 h cMpc −1 is found, an improvement by a factor ≈ 8 in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.

The 21 cm brightness temperature δT b fluctuations from reionization promise to provide information on the physical processes during that epoch. We present a formalism for generating the δT b distribution using dark matter simulations and an one-dimensional radiative transfer code. Our analysis is able to account for the spin temperature T S fluctuations arising from inhomogeneous X-ray heating and Lyα coupling during cosmic dawn. The δT b power spectrum amplitude at large scales (k ∼ 0.1 Mpc −1 ) is maximum when ∼ 10% of the gas (by volume) is heated above the cosmic microwave background temperature. The power spectrum shows a "bump"-like feature during cosmic dawn and its location measures the typical sizes of heated regions. We find that the effect of peculiar velocities on the power spectrum is negligible at large scales for most part of the reionization history. During early stages (when the volume averaged ionization fraction 0.2) this is because the signal is dominated by fluctuations in T S . For reionization models that are solely driven by stars within high mass ( 10 9 M ⊙ ) haloes, the peculiar velocity effects are prominent only at smaller scales (k 0.4 Mpc −1 ) where patchiness in the neutral hydrogen density dominates the signal. The conclusions are unaffected by changes in the amplitude or steepness in the X-ray spectra of the sources.

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift ≈ 9.1 using new upper limits on the 21-cm power spectrum measured by the LO-FAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction 0.13 and a distribution of the ionized regions with a characteristic size 8 h −1 comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) 16 h −1 Mpc is ruled out. For an IGM with a uniform spin temperature T S 3 K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction 0.34 of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 h −1 Mpc and FWHM of 110 h −1 Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.

The upcoming radio interferometer Square Kilometre Array (SKA) is expected to directly detect the redshifted 21-cm signal from the neutral hydrogen present during the Cosmic Dawn. Temperature fluctuations from X-ray heating of the neutral intergalactic medium can dominate the fluctuations in the 21-cm signal from this time. This heating depends on the abundance, clustering, and properties of the X-ray sources present, which remain highly uncertain. We present a suite of three new large-volume, 349 Mpc a side, fully numerical radiative transfer simulations including QSO-like sources, extending the work previously presented in Ross et al. (2017). The results show that our QSOs have a modest contribution to the heating budget, yet significantly impact the 21-cm signal. Initially, the power spectrum is boosted on large scales by heating from the biased QSO-like sources, before decreasing on all scales. Fluctuations from images of the 21-cm signal with resolutions corresponding to SKA1-Low at the appropriate redshifts are well above the expected noise for deep integrations, indicating that imaging could be feasible for all the X-ray source models considered. The most notable contribution of the QSOs is a dramatic increase in non-Gaussianity of the signal, as measured by the skewness and kurtosis of the 21-cm probability distribution functions. However, in the case of late Lyman-α saturation, this non-Gaussianity could be dramatically decreased particularly when heating occurs earlier. We conclude that increased non-Gaussianity is a promising signature of rare X-ray sources at this time, provided that Lyman-α saturation occurs before heating dominates the 21-cm signal.

Details of various unknown physical processes during the cosmic dawn and the epoch of reionization can be extracted from observations of the redshifted 21-cm signal. These observations, however, will be affected by the evolution of the signal along the line-ofsight which is known as the "light-cone effect". We model this effect by post-processing a dark matter N −body simulation with an 1-D radiative transfer code. We find that the effect is much stronger and dramatic in presence of inhomogeneous heating and Lyα coupling compared to the case where these processes are not accounted for. One finds increase (decrease) in the spherically averaged power spectrum up to a factor of 3 (0.6) at large scales (k ∼ 0.05 Mpc −1 ) when the light-cone effect is included, though these numbers are highly dependent on the source model. The effect is particularly significant near the peak and dip-like features seen in the power spectrum. The peaks and dips are suppressed and thus the power spectrum can be smoothed out to a large extent if the width of the frequency band used in the experiment is large. We argue that it is important to account for the light-cone effect for any 21-cm signal prediction during cosmic dawn.

The ability of the future low frequency component of the Square Kilometre Array radio telescope (SKA-Low) to produce tomographic images of the redshifted 21-cm signal will enable direct studies of the evolution of the sizes and shapes of ionized regions during the Epoch of Reionization. However, a reliable identification of ionized regions in noisy interferometric data is not trivial. Here, we introduce an image processing method known as superpixels for this purpose. We compare this method with two other previously proposed ones, one relying on a chosen threshold and the other employing automatic threshold determination using the K-Means algorithm. We use a correlation test and compare power spectra and bubble size distributions to show that the superpixels method provides a better identification of ionized regions, especially in the case of noisy data. We also describe some possible additional applications of the superpixel method, namely the derivation of the ionization history and constraints on the source properties in specific regions.

Understanding properties of the first sources in the Universe using the redshifted H i 21-cm signal is one of the major aims of present and upcoming low-frequency experiments. We investigate the possibility of imaging the redshifted 21-cm pattern around the first sources during the cosmic dawn using the SKA1-low. We model the H i 21-cm image maps, appropriate for the SKA1-low, around the first sources consisting of stars and X-ray sources within galaxies. In addition to the system noise, we account also for the astrophysical foregrounds by adding them to the signal maps. We find that after subtracting the foregrounds using a polynomial fit and suppressing the noise by smoothing the maps over 10 ′ − 30 ′ angular scale, the isolated sources at z ∼ 15 are detectable with ∼ 4 − 9 σ confidence level in 2000 h of observation with the SKA1-low. Although the 21-cm profiles around the sources get altered because of the Gaussian smoothing, the images can still be used to extract some of the source properties. We account for overlaps in the patterns of the individual sources by generating realistic H i 21-cm maps of the cosmic dawn that are based on N -body simulations and a onedimensional radiative transfer code. We find that these sources should be detectable in the SKA1-low images at z = 15 with an SNR of ∼ 14(4) in 2000 (200) h of observations. One possible observational strategy thus could be to observe multiple fields for shorter observation times, identify fields with SNR 3 and observe these fields for much longer duration. Such observations are expected to be useful in constraining the parameters related to the first sources.

We compare various foreground removal techniques that are being utilized to remove bright foregrounds in various experiments aiming to detect the redshifted 21 cm signal of neutral hydrogen from the epoch of reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limit on the 21 cm signal. Interestingly, we find that GMCA and FastICA reproduce the most recent 2σ upper limit of $\Delta ^2_{21} \lt $ (73)2 mK2 at k = 0.075 hcMpc−1, which resulted from the application of GPR. We also find that FastICA and GMCA begin to deviate from the noise-limit at k-scales larger than ∼0.1 hcMpc−1. We then replicate the data via simulations to see the source of FastICA and GMCA’s limitations, by testing them against various instrumental effects. We find that no single instrumental effect, such as primary beam effects or mode-mixing, can explain the poorer recovery by FastICA and GMCA at larger k-scales. We then test scale-independence of FastICA and GMCA, and find that lower k-scales can be modelled by a smaller number of independent components. For larger scales (k ≳ 0.1 hcMpc−1), more independent components are needed to fit the foregrounds. We conclude that, the current usage of GPR by the LOFAR collaboration is the appropriate removal technique. It is both robust and less prone to overfitting, with future improvements to GPR’s fitting optimization to yield deeper limits.