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

Measurements of the Epoch of Reionization (EoR) 21-cm signal hold the potential to constrain models of reionization. In this paper we consider a reionization model with three astrophysical parameters namely (1) the minimum halo mass which can host ionizing sources, $M_{\rm min}$, (2) the number of ionizing photons escaping into the IGM per baryon within the halo, $N_{\rm ion}$ and (3) the mean free path of the ionizing photons within the IGM, $R_{\rm mfp}$. We predict the accuracy with which these parameters can be measured from future observations of the 21-cm power spectrum (PS) using the upcoming SKA-Low. Unlike several earlier works, we account for the non-Gaussianity of the inherent EoR 21-cm signal. Considering cosmic variance only and assuming that foregrounds are completely removed, we find that non-Gaussianity increases the volume of the $1 \sigma$ error ellipsoid of the parameters by a factor of $133$ relative to the Gaussian predictions, the orientation is also different. The ratio of the volume of error ellipsoids is $1.65$ and $2.67$ for observation times of $1024$ and $10000$ hours respectively, when all the $\mathbf{k}$ modes within the foreground wedge are excluded. With foreground wedge excluded and for $1024$ hours, the 1D marginalized errors are $(\Delta M_{\rm min}/M_{\rm min},\Delta N_{\rm ion}/N_{\rm ion},\Delta R_{\rm mfp}/R_{\rm mfp})=(6.54, 2.71, 7.75) \times 10^{-2}$ which are respectively $2 \%$, $5 \%$ and $23 \%$ larger than the respective Gaussian predictions. The impact of non-Gaussianity increases for longer observations, and it is particularly important for $R_{\rm mfp}$.

The light-cone (LC) effect causes the mean as well as the statistical properties of the redshifted 21-cm signal $T_{\rm b}(\hat{\bf n},\nu)$ to change with frequency $\nu$ (or cosmic time). Consequently, the statistical homogeneity (ergodicity) of the signal along the line of sight (LoS) direction is broken. This is a severe problem particularly during the Epoch of Reionization (EoR) when the mean neutral hydrogen fraction ($\bar{x}_{\rm HI}$) changes rapidly as the universe evolves. This will also pose complications for large bandwidth observations. These effects imply that the 3D power spectrum $P(k)$ fails to quantify the entire second-order statistics of the signal as it assumes the signal to be ergodic and periodic along the LoS. As a proper alternative to $P(k)$, we use the multi-frequency angular power spectrum (MAPS) ${\mathcal C}_{\ell}(\nu_1,\nu_2)$ which does not assume the signal to be ergodic and periodic along the LoS. Here, we study the prospects for measuring the EoR 21-cm MAPS using future observations with the upcoming SKA-Low. Ignoring any contribution from the foregrounds, we find that the EoR 21-cm MAPS can be measured at a confidence level $\ge 5\sigma$ at angular scales $\ell \sim 1300$ for total observation time $t_{\rm obs} \ge 128\,{\rm hrs}$ across $\sim 44\,{\rm MHz}$ observational bandwidth. We also quantitatively address the effects of foregrounds on MAPS detectability forecast by avoiding signal contained within the foreground wedge in $(k_\perp, k_\parallel)$ plane. These results are very relevant for the upcoming large bandwidth EoR experiments as previous predictions were all restricted to individually analyzing the signal over small frequency (or equivalently redshift) intervals.

The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at $1.42\,$GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen (Fialkov & Barkana 2019). We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on $141\,$hours of observations with LOFAR (Mertens et al. 2020) to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 68% C.L.) that the cosmological radio background can be at most 0.7% of the CMB at $1.42\,$GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-0.7$% over the CMB (at $1.42\,$GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$, and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data disfavour scenarios with inefficient X-ray sources and cases where the Universe was ionized by massive halos only.

Minkowski functionals and Shapefinders shed light on the connectedness of large-scale structure by determining its topology and morphology. We use a sophisticated code, SURFGEN2, to measure the Minkowski functionals and Shapefinders of individual clusters by modelling cluster surfaces using the 'Marching Cube 33' triangulation algorithm. In this paper, we study the morphology of simulated neutral hydrogen (HI) density fields using Shapefinders at various stages of reionization from the excursion set approach. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS) to understand the percolation process. Percolation curves demonstrate that the non-Gaussianity in the HI field increases as reionization progresses. The large clusters in both the HI overdense and underdense excursion sets possess similar values of "thickness" ($T$), as well as "breadth" ($B$), but their third Shapefinder - "length" ($L$) - becomes almost proportional to their volume. The large clusters in both HI overdense and underdense segments are overwhelmingly filamentary. The 'cross-section' of a filamentary cluster can be estimated using the product of the first two Shapefinders, $T \times B$. Hence the cross sections of the large clusters at the onset of percolation do not vary much with volume and their sizes only differ in terms of their lengths. This feature appears more vividly in HI overdense regions than in underdense regions and is more pronounced at lower redshifts which correspond to an advanced stage of reionization.

A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 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-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 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.

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR 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 $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 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 $\lesssim 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 $\lesssim 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 Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

Recent simulations show the Epoch of Reionization (EoR) 21-cm signal to be inherently non-Gaussian whereby the error covariance matrix $\mathbf{C}_{ij}$ of the 21-cm power spectrum (PS) contains a trispectrum contribution that would be absent if the signal were Gaussian. Using the binned power spectrum and trispectrum from simulations, here we present a methodology for incorporating these with the baseline distribution and system noise to make error predictions for observations with any radio-interferometric array. Here we consider the upcoming SKA-Low. Non-Gaussianity enhances the errors introducing a positive deviation $\Delta$ relative to the Gaussian predictions. $\Delta$ increases with observation time $t_{\rm obs}$ and saturates as the errors approach the cosmic variance. Considering $t_{\rm obs}=1024$ hours where a $5 \sigma$ detection is possible at all redshifts $7 \le z \le 13$, in the absence of foregrounds we find that the deviations are important at small $k$ where we have $\Delta \sim 40-100 \%$ at $k~\sim 0.04 ~{\rm Mpc}^{-1}$ for some of the redshifts and also at intermediate $k \, (\sim 0.4 ~{\rm Mpc}^{-1})$ where we have $\Delta \sim 200 \%$ at $z=7$. Non-Gaussianity also introduces correlations between the errors in different $k$ bins, and we find both correlations and anticorrelations with the correlation coefficient value spanning $-0.4 \le r_{ij} \le 0.8$. Incorporating the foreground wedge, $\Delta$ continues to be important ($> 50\%$) at $z=7$. We conclude that non-Gaussianity makes a significant contribution to the errors and this is important in the context of the future instruments that aim to achieve high-sensitivity measurements of the EoR 21-cm PS.

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.

The light-cone (LC) effect imprints the cosmological evolution of the redshifted 21-cm signal $T_{\rm b} ({\hat{\bf{n}}}, \nu)$ along the frequency axis which is the line of sight (LoS) direction of an observer. The effect is particularly pronounced during the Epoch of Reionization (EoR) when the mean hydrogen neutral fraction ${\bar{x}_{\rm HI}}(\nu)$ falls rapidly as the universe evolves. The multi-frequency angular power spectrum (MAPS) ${{\mathcal C}_{\ell}}(\nu_1,\nu_2)$ quantifies the entire second-order statistics of $T_{\rm b} ({\hat{\bf{n}}}, \nu)$ considering both the systematic variation along $\nu$ due to the cosmological evolution and also the statistically homogeneous and isotropic fluctuations along all the three spatial directions encoded in ${\hat{\bf{n}}}$ and $\nu$. Here we propose a simple model where the systematic frequency $(\nu_1,\nu_2)$ dependence of ${{\mathcal C}_{\ell}}(\nu_1,\nu_2)$ arises entirely due to the evolution of ${\bar{x}_{\rm HI}}(\nu)$. This provides a new method to observationally determine the reionization history. Considering a LC simulation of the EoR 21-cm signal, we use the diagonal elements $\nu_1=\nu_2$ of ${{\mathcal C}_{\ell}}(\nu_1,\nu_2)$ to validate our model. We demonstrate that it is possible to recover the reionization history across the entire observational bandwidth provided we have the value ${\bar{x}_{\rm HI}}$ at a single frequency as an external input.

Using Shapefinders, which are ratios of Minkowski functionals, we study the morphology of neutral hydrogen (HI) density fields, simulated using semi-numerical technique (inside-out), at various stages of reionization. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS), originally proposed in Klypin and Shandarin (1993), to study the percolation in both neutral and ionized hydrogen. We find that the largest ionized region is percolating below the neutral fraction $x_{HI} \lesssim 0.728$ (or equivalently $z \lesssim 9$). The study of Shapefinders reveals that the largest ionized region starts to become highly filamentary with non-trivial topology near the percolation transition. During the percolation transition, the first two Shapefinders - 'thickness' ($T$) and 'breadth' ($B$) - of the largest ionized region do not vary much, while the third Shapefinder - 'length' ($L$) - abruptly increases. Consequently, the largest ionized region tends to be highly filamentary and topologically quite complex. The product of the first two Shapefinders, $T\times B$, provides a measure of the 'cross-section' of a filament-like ionized region. We find that, near percolation, the value of $T\times B$ for the largest ionized region remains stable at $\sim 7$ Mpc$^2$ (in comoving scale) while its length increases with time. Interestingly all large ionized regions have similar cross-sections. However their length shows a power-law dependence on their volume, $L\propto V^{0.72}$, at the onset of percolation.

An intriguing possibility for the dark sector of our universe is that the dark matter particle could interact with a dark radiation component. If the non-gravitational interactions of the dark matter and dark radiation species with Standard Model particles are highly suppressed, then astrophysics and cosmology could be our only windows into probing the dynamics of such a dark sector. It is well known that such dark sectors would lead to suppression of small scale structure, which would be constrained by measurements of the Lyman-$\alpha$ forest. In this work we consider the cosmological signatures of such dark sectors on the reionization history of our universe. Working within the "ETHOS" (effective theory of structure formation) framework, we show that if such a dark sector exists in our universe, the suppression of low mass dark matter halos would also reduce the total number of ionizing photons, thus affecting the reionization history of our universe. We place constraints on the interaction strengths within such dark sectors by using the measured value of the optical depth from the Planck satellite, as well as from demanding a successful reionization history. We compare and contrast such scenarios with warm dark matter scenarios which also suppress structure formation on small scales. In a model where dark matter interacts with a sterile neutrino, we find a bound on the ETHOS parameter $a_4\lesssim 1.2\times10^6\textrm{ Mpc}^{-1}$. For warm dark matter models, we constrain the mass $m_{\textrm{WDM}}\gtrsim 0.7\textrm{ keV}$, which is comparable to bounds obtained from Lyman-$\alpha$ measurements. Future 21-cm experiments will measure the global history of reionization and the neutral hydrogen power spectrum, which could either lead to stronger constraints or discovery of secret dark sector interactions.

The epoch of reionization (EoR) 21-cm signal is expected to be highly non-Gaussian in nature and this non-Gaussianity is also expected to evolve with the progressing state of reionization. Therefore the signal will be correlated between different Fourier modes ($k$). The power spectrum will not be able capture this correlation in the signal. We use a higher-order estimator -- the bispectrum -- to quantify this evolving non-Gaussianity. We study the bispectrum using an ensemble of simulated 21-cm signal and with a large variety of $k$ triangles. We observe two competing sources driving the non-Gaussianity in the signal: fluctuations in the neutral fraction ($x_{\rm HI}$) field and fluctuations in the matter density field. We find that the non-Gaussian contribution from these two sources vary, depending on the stage of reionization and on which $k$ modes are being studied. We show that the sign of the bispectrum works as a unique marker to identify which among these two components is driving the non-Gaussianity. We propose that the sign change in the bispectrum, when plotted as a function of triangle configuration $\cos{\theta}$ and at a certain stage of the EoR can be used as a confirmative test for the detection of the 21-cm signal. We also propose a new consolidated way to visualize the signal evolution (with evolving $\overline{x}_{\rm HI}$ or redshift), through the trajectories of the signal in a power spectrum and equilateral bispectrum i.e. $P(k)-B(k, k, k)$ space.

The light-cone (LC) effect causes the Epoch of Reionization (EoR) 21-cm signal $T_{\rm b} (\hat{\bf{n}}, \nu)$ to evolve significantly along the line of sight (LoS) direction $\nu$. In the first part of this paper, we present a method to properly incorporate the LC effect in simulations of the EoR 21-cm signal that include peculiar velocities. Subsequently, we discuss how to quantify the second order statistics of the EoR 21-cm signal in the presence of the LC effect. We demonstrate that the 3D power spectrum $P({\bf{k}})$ fails to quantify the entire information because it assumes the signal to be ergodic and periodic, whereas the LC effect breaks these conditions along the LoS. Considering a LC simulation centered at redshift $8$ where the mean neutral fraction drops from $0.65$ to $0.35$ across the box, we find that $P({\bf{k}})$ misses out $\sim 40 \%$ of the information at the two ends of the $17.41 \, {\rm MHz}$ simulation bandwidth. The multi-frequency angular power spectrum (MAPS) ${\mathcal C}_{\ell}(\nu_1,\nu_2)$ quantifies the statistical properties of $T_{\rm b} (\hat{\bf{n}}, \nu)$ without assuming the signal to be ergodic and periodic along the LoS. We expect this to quantify the entire statistical information of the EoR 21-cm signal. We apply MAPS to our LC simulation and present preliminary results for the EoR 21-cm signal.

In this paper we establish the accuracy and robustness of a fast estimator for the bispectrum - the "FFT bispectrum estimator". The implementation of the estimator presented here offers speed and simplicity benefits over a direct sampling approach. We also generalise the derivation so it may be easily be applied to any order polyspectra, such as the trispectrum, with the cost of only a handful of FFTs. All lower order statistics can also be calculated simultaneously for little extra cost. To test the estimator we make use of a non-linear density field, and for a more strongly non-Gaussian test case we use a toy-model of reionization in which ionized bubbles at a given redshift are all of equal size and are randomly distributed. Our tests find that the FFT estimator remains accurate over a wide range of k, and so should be extremely useful for analysis of 21-cm observations. The speed of the FFT bispectrum estimator makes it suitable for sampling applications, such as Bayesian inference. The algorithm we describe should prove valuable in the analysis of simulations and observations, and whilst we apply it within the field of cosmology, this estimator is useful in any field that deals with non-Gaussian data.

The EoR 21-cm signal is expected to become highly non-Gaussian as reionization progresses. This severely affects the error-covariance of the EoR 21-cm power spectrum which is important for predicting the prospects of a detection with ongoing and future experiments. Most earlier works have assumed that the EoR 21-cm signal is a Gaussian random field where (1) the error variance depends only on the power spectrum and the number of Fourier modes in the particular $k$ bin, and (2) the errors in the different $k$ bins are uncorrelated. Here we use an ensemble of simulated 21-cm maps to analyze the error-covariance at various stages of reionization. We find that even at the very early stages of reionization ($\bar{x}_{\rm HI} \sim 0.9 $) the error variance significantly exceeds the Gaussian predictions at small length-scales ($k > 0.5 \,{\rm Mpc}^{-1}$) while they are consistent at larger scales. The errors in most $k$ bins (both large and small scales), are however found to be correlated. Considering the later stages ($\bar{x}_{\rm HI} = 0.15$), the error variance shows an excess in all $k$ bins within $k \ge 0.1 \, {\rm Mpc}^{-1}$, and it is around $200$ times larger than the Gaussian prediction at $k \sim 1 \, {\rm Mpc}^{-1}$. The errors in the different $k$ bins are all also highly correlated, barring the two smallest $k$ bins which are anti-correlated with the other bins. Our results imply that the predictions for different 21-cm experiments based on the Gaussian assumption underestimate the errors, and it is necessary to incorporate the non-Gaussianity for more realistic predictions.

Studying the cosmic dawn and the epoch of reionization through the redshifted 21 cm line are among the major science goals of the SKA1. Their significance lies in the fact that they are closely related to the very first stars in the universe. Interpreting the upcoming data would require detailed modelling of the relevant physical processes. In this article, we focus on the theoretical models of reionization that have been worked out by various groups working in India with the upcoming SKA in mind. These models include purely analytical and semi-numerical calculations as well as fully numerical radiative transfer simulations. The predictions of the 21 cm signal from these models would be useful in constraining the properties of the early galaxies using the SKA data.

The line of sight direction in the redshifted 21-cm signal coming from the cosmic dawn and the epoch of reionization is quite unique in many ways compared to any other cosmological signal. Different unique effects, such as the evolution history of the signal, non-linear peculiar velocities of the matter etc will imprint their signature along the line of sight axis of the observed signal. One of the major goals of the future SKA-LOW radio interferometer is to observe the cosmic dawn and the epoch of reionization through this 21-cm signal. It is thus important to understand how these various effects affect the signal for it's actual detection and proper interpretation. For more than one and half decades, various groups in India have been actively trying to understand and quantify the different line of sight effects that are present in this signal through analytical models and simulations. In many ways the importance of this sub-field under 21-cm cosmology have been identified, highlighted and pushed forward by the Indian community. In this article we briefly describe their contribution and implication of these effects in the context of the future surveys of the cosmic dawn and the epoch of reionization that will be conducted by the SKA-LOW.

The particle nature of dark matter remains a mystery. In this paper, we consider two dark matter models---Late Forming Dark Matter (LFDM) and Ultra-Light Axion (ULA) models---where the matter power spectra show novel effects on small scales. The high redshift universe offers a powerful probe of their parameters. In particular, we study two cosmological observables: the neutral hydrogen (HI) redshifted 21-cm signal from the epoch of reionization, and the evolution of the collapsed fraction of HI in the redshift range $2 < z < 5$. We model the theoretical predictions of the models using CDM-like N-body simulations with modified initial conditions, and generate reionization fields using an excursion-set model. The N-body approximation is valid on the length and halo mass scales studied. We show that LFDM and ULA models predict an increase in the HI power spectrum from the epoch of reionization by a factor between 2--10 for a range of scales $0.1<k<4 \, \rm Mpc^{-1}$. Assuming a fiducial model where a neutral hydrogen fraction $\bar{x}_{HI}=0.5$ must be achieved by $z=8$, the reionization process allows us to put approximate bounds on the redshift of dark matter formation $z_f > 4 \times 10^5$ (for LFDM) and the axion mass $m_a > 2.6 \times 10^{-23} \, \rm eV$ (for ULA). The comparison of the collapsed mass fraction inferred from damped Lyman-$\alpha$ observations to the theoretical predictions of our models lead to the weaker bounds: $z_f > 2 \times 10^5$ and $m_a > 10^{-23} \, \rm eV$. These bounds are consistent with other constraints in the literature using different observables; we briefly discuss how these bounds compare with possible constraints from the observation of luminosity function of galaxies at high redshifts. In the case of ULAs, these constraints are also consistent with a solution to the cusp-core problem of CDM.

The non-Gaussian nature of the epoch of reionization (EoR) 21-cm signal has a significant impact on the error variance of its power spectrum $P({\bf \textit{k}})$. We have used a large ensemble of semi-numerical simulations and an analytical model to estimate the effect of this non-Gaussianity on the entire error-covariance matrix ${\mathcal{C}}_{ij}$. Our analytical model shows that ${\mathcal{C}}_{ij}$ has contributions from two sources. One is the usual variance for a Gaussian random field which scales inversely of the number of modes that goes into the estimation of $P({\bf \textit{k}})$. The other is the trispectrum of the signal. Using the simulated 21-cm signal ensemble, an ensemble of the randomized signal and ensembles of Gaussian random ensembles we have quantified the effect of the trispectrum on the error variance ${\mathcal{C}}_{ij}$. We find that its relative contribution is comparable to or larger than that of the Gaussian term for the $k$ range $0.3 \leq k \leq 1.0 \,{\rm Mpc}^{-1}$, and can be even $\sim 200$ times larger at $k \sim 5\, {\rm Mpc}^{-1}$. We also establish that the off-diagonal terms of ${\mathcal{C}}_{ij}$ have statistically significant non-zero values which arise purely from the trispectrum. This further signifies that the error in different $k$ modes are not independent. We find a strong correlation between the errors at large $k$ values ($\ge 0.5 \,{\rm Mpc}^{-1}$), and a weak correlation between the smallest and largest $k$ values. There is also a small anti-correlation between the errors in the smallest and intermediate $k$ values. These results are relevant for the $k$ range that will be probed by the current and upcoming EoR 21-cm experiments.

The Epoch of Reionization (EoR) 21-cm signal is expected to become increasingly non-Gaussian as reionization proceeds. We have used semi-numerical simulations to study how this affects the error predictions for the EoR 21-cm power spectrum. We expect $SNR=\sqrt{N_k}$ for a Gaussian random field where $N_k$ is the number of Fourier modes in each $k$ bin. We find that non-Gaussianity is important at high $SNR$ where it imposes an upper limit $[SNR]_l$. For a fixed volume $V$, it is not possible to achieve $SNR > [SNR]_l$ even if $N_k$ is increased. The value of $[SNR]_l$ falls as reionization proceeds, dropping from $\sim 500$ at $\bar{x}_{HI} = 0.8-0.9$ to $\sim 10$ at $\bar{x}_{HI} = 0.15 $ for a $[150.08\, {\rm Mpc}]^3$ simulation. We show that it is possible to interpret $[SNR]_l$ in terms of the trispectrum, and we expect $[SNR]_l \propto \sqrt{V}$ if the volume is increased. For $SNR \ll [SNR]_l$ we find $SNR = \sqrt{N_k}/A $ with $A \sim 0.95 - 1.75$, roughly consistent with the Gaussian prediction. We present a fitting formula for the $SNR$ as a function of $N_k$, with two parameters $A$ and $[SNR]_l$ that have to be determined using simulations. Our results are relevant for predicting the sensitivity of different instruments to measure the EoR 21-cm power spectrum, which till date have been largely based on the Gaussian assumption.