A possible way to study the reionization of cosmic hydrogen is by observing the large ionized regions (bubbles) around bright individual sources, e.g., quasars, using the redshifted 21 cm signal. It has already been shown that matched filter-based methods are not only able to detect the weak 21 cm signal from these bubbles but also aid in constraining their properties. In this work, we extend the previous studies to develop a rigorous Bayesian framework to explore the possibility of constraining the parameters that characterize the bubbles. To check the accuracy with which we can recover the bubble parameters, we apply our method on mock observations appropriate for the upcoming SKA1-low. For a region of size $\gtrsim 50$ cMpc around a typical quasar at redshift 7, we find that $\approx 20$ h of integration with SKA1-low will be able to constrain the size and location of the bubbles, as well as the difference in the neutral hydrogen fraction inside and outside the bubble, with $\lesssim 10\%$ precision. The recovery of the parameters are more precise and the SNR of the detected signal is higher when the bubble sizes are larger and their shapes are close to spherical. Our method can be useful in identifying regions in the observed field which contain large ionized regions and hence are interesting for following up with deeper integration times.

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

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

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 resonance scattering of Ly$\alpha$ photons with neutral hydrogen atoms in the intergalactic medium not only couples the spin temperature to the kinetic temperature but also leads to a heating of the gas. We investigate the impact of this heating on the average brightness temperature of the 21-cm signal from the Cosmic Dawn in the context of the claimed detection by the EDGES low-band experiment. We model the evolution of the global signal taking into account the Ly$\alpha$ coupling and heating and a cooling which can be stronger than the Hubble cooling. Using the claimed detection of a strong absorption signal at $z\approx 17$ as a constraint, we find that a strong Ly$\alpha$ background is ruled out. Instead the results favour a weak Ly$\alpha$ background combined with an excess cooling mechanism which is substantially stronger than previously considered. We also show that the cooling mechanism driven by the interaction between millicharged baryons and dark matter particles no longer provides a viable explanation for the EDGES result when Ly$\alpha$ heating is taken into account.

We have investigated the long term average spectral properties of galactic X-ray binaries in the energy range of 3-20 keV, using long term monitoring data from MAXI -Gas Slit Camera (GSC). These long term average spectra are used to construct separately the composite spectra of galactic High Mass X-ray binaries (HMXBs) and Low Mass X-ray binaries (LMXBs). These composite spectra can be described empirically with piece-wise power-law with three components. X-rays from HMXBs are considered as important contributors to heating and ionization of neutral hydrogen in the intergalactic medium during the Epoch of Reionization. Using the above empirical form of the composite HMXB spectra extrapolated to lower energies as an input, we have studied the impact of these sources on the 21-cm signal using the outputs of N-body simulation and 1D radiative transfer. The heating due to the composite spectrum is less patchy compared to power-law spectrum with a spectral index $\alpha = 1.5$, used in previous studies. The amplitude of the heating peak of large scale power spectrum, when plotted as a function of the redshift, is less for the composite spectrum.

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-$\alpha$ 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-$\alpha$ saturation occurs before heating dominates the 21-cm signal.

Potential small-scale discrepancies in the picture of galaxy formation painted by the $\Lambda$CDM paradigm have led to considerations of modified dark matter models. One such dark matter model that has recently attracted much attention is fuzzy dark matter (FDM). In FDM models, the dark matter is envisaged to be an ultra-light scalar field with a particle mass $m_{\rm FDM} \sim 10^{-22} $ eV. This yields astronomically large de Broglie wavelengths which can suppress small-scale structure formation and give rise to the observed kpc-sized density cores in dwarf galaxies. We investigate the evolution of the 21-cm signal during Cosmic Dawn and the Epoch of Reionization (EoR) in $\Lambda$FDM cosmologies using analytical models. The delay in source formation and the absence of small halos in $\Lambda$FDM significantly postpone the Ly$\alpha$ coupling, heating, as well as the reionization of the neutral hydrogen of the intergalactic medium. As a result, the absorption feature in the evolution of the global 21-cm signal has a significantly smaller full width at half maximum ($\Delta z \lesssim 3$), than $\Lambda$CDM ($\Delta z \simeq 6$). This alone rules out $m_{\rm FDM} < 6 \times 10^{-22}$ eV as a result of the $2\sigma$ lower limit $\Delta z \gtrsim 4$ from EDGES High-Band. As a result, $\Lambda$FDM is not a viable solution to the potential small-scale problems facing $\Lambda$CDM. Finally, we show that any detection of the 21-cm signal at redshifts $z > 14$ by interferometers such as the SKA can also exclude $\Lambda$FDM models.

The 21-cm signal from the epoch of reionization is non-Gaussian. Current radio telescopes are focused on detecting the 21-cm power spectrum, but in the future the Square Kilometre Array is anticipated to provide a first measurement of the bispectrum. Previous studies have shown that the position-dependent power spectrum is a simple and efficient way to probe the squeezed-limit bispectrum. In this approach, the survey is divided into subvolumes and the correlation between the local power spectrum and the corresponding mean density of the subvolume is computed. This correlation is equivalent to an integral of the bispectrum in the squeezed limit, but is much simpler to implement than the usual bispectrum estimators. It also has a clear physical interpretation: it describes how the small-scale power spectrum of tracers such as galaxies and the 21-cm signal respond to a large-scale environment. Reionization naturally couples large and small scales as ionizing radiation produced by galactic sources can travel up to tens of Megaparsecs through the intergalactic medium during this process. Here we apply the position-dependent power spectrum approach to fluctuations in the 21-cm background from reionization. We show that this statistic has a distinctive evolution in time that can be understood with a simple analytic model. We also show that the statistic can easily distinguish between simple "inside-out" and "outside-in" models of reionization. The position-dependent power spectrum is thus a promising method to validate the reionization signal and to extract higher-order information on this process.

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.

One of the main challenges in probing the reionization epoch using the redshifted 21 cm line is that the magnitude of the signal is several orders smaller than the astrophysical foregrounds. One of the methods to deal with the problem is to avoid a wedge-shaped region in the Fourier $k_{\perp} - k_{\parallel}$ space which contains the signal from the spectrally smooth foregrounds. However, measuring the spherically averaged power spectrum using only modes outside this wedge (i.e., in the reionization window), leads to a bias. We provide a prescription, based on expanding the power spectrum in terms of the shifted Legendre polynomials, which can be used to compute the angular moments of the power spectrum in the reionization window. The prescription requires computation of the monopole, quadrupole and hexadecapole moments of the power spectrum using the theoretical model under consideration and also the knowledge of the effective extent of the foreground wedge in the $k_{\perp} - k_{\parallel}$ plane. One can then calculate the theoretical power spectrum in the window which can be directly compared with observations. The analysis should have implications for avoiding any bias in the parameter constraints using 21 cm power spectrum data.

Three-dimensional radiative transfer simulations of the epoch of reionization can produce realistic results, but are computationally expensive. On the other hand, simulations relying on one-dimensional radiative transfer solutions are faster but limited in accuracy due to their more approximate nature. Here, we compare the performance of the reionization simulation codes grizzly and C2-ray which use 1D and 3D radiative transfer schemes respectively. The comparison is performed using the same cosmological density fields, halo catalogues and source properties. We find that the ionization maps, as well as the 21-cm signal maps from these two simulations are very similar even for complex scenarios which include thermal feedback on low mass halos. The comparison between the schemes in terms of the statistical quantities such as the power spectrum of the brightness temperature fluctuation agree with each other within 10% error throughout the entire reionization history. grizzly seems to perform slightly better than the semi-numerical approaches considered in Majumdar et al. (2014) which are based on the excursion set principle. We argue that grizzly can be efficiently used for exploring parameter space, establishing observations strategies and estimating parameters from 21-cm observations.

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.

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.

Detection of individual luminous sources during the reionization epoch and cosmic dawn through their signatures in the HI 21-cm signal is one of the direct approaches to probe the epoch. Here, we summarize our previous works on this and present preliminary results on the prospects of detecting such sources using the SKA1-low experiment. We first discuss the expected HI 21-cm signal around luminous sources at different stages of reionization and cosmic dawn. We then introduce two visibility based estimators for detecting such signal: one based on the matched filtering technique and the other relies on simply combing the visibility signal from different baselines and frequency channels. We find that that the SKA1-low should be able to detect ionized bubbles of radius $R_b \gtrsim 10$ Mpc with $\sim 100$ hr of observations at redshift $z \sim 8$ provided that the mean outside neutral Hydrogen fraction $ x_{\rm HI} \gtrsim 0.5$. We also investigate the possibility of detecting HII regions around known bright QSOs such as around ULASJ1120+0641 discovered by Mortlock et al. 2011. We find that a $5 \sigma$ detection is possible with $600$ hr of SKA1-low observations if the QSO age and the outside $ x_{\rm HI} $ are at least $\sim 2 \times 10^7$ Myr and $\sim 0.2$ respectively. Finally, we investigate the possibility of detecting the very first X-ray and Ly-$\alpha$ sources during the cosmic dawn. We consider mini-QSOs like sources which emits in X-ray frequency band. We find that with a total $\sim 1000$ hr of observations, SKA1-low should be able to detect those sources individually with a $\sim 9 \sigma$ significance at redshift $z=15$. We summarize how the SNR changes with various parameters related to the source properties.

Understanding properties of the first sources in the Universe using the redshifted \HI ~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 \HI ~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 \sim 15$ are detectable with $\sim 4 - 9 \, \sigma$ 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 \HI ~21-cm maps of the cosmic dawn that are based on $N$-body simulations and a one-dimensional radiative transfer code. We find that these sources should be detectable in the SKA1-low images at $z = 15$ with an SNR of $\sim 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 $\gtrsim 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.

Currently several low-frequency experiments are being planned to study the nature of the first stars using the redshifted 21-cm signal from the cosmic dawn and epoch of reionization. Using a one-dimensional radiative transfer code, we model the 21-cm signal pattern around the early sources for different source models, i.e., the metal-free Population III (PopIII) stars, primordial galaxies consisting of Population II (PopII) stars, mini-QSOs and high-mass X-ray binaries (HMXBs). We investigate the detectability of these sources by comparing the 21-cm visibility signal with the system noise appropriate for a telescope like the SKA1-low. Upon integrating the visibility around a typical source over all baselines and over a frequency interval of 16 MHz, we find that it will be possible to make a $\sim 9-\sigma$ detection of the isolated sources like PopII galaxies, mini-QSOs and HMXBs at $z \sim 15$ with the SKA1-low in 1000 hours. The exact value of the signal to noise ratio (SNR) will depend on the source properties, in particular on the mass and age of the source and the escape fraction of ionizing photons. The predicted SNR decreases with increasing redshift. We provide simple scaling laws to estimate the SNR for different values of the parameters which characterize the source and the surrounding medium. We also argue that it will be possible to achieve a SNR $\sim 9$ even in the presence of the astrophysical foregrounds by subtracting out the frequency-independent component of the observed signal. These calculations will be useful in planning 21-cm observations to detect the first sources.

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-of-sight 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$\alpha$ 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 \sim 0.05\, \rm 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 21 cm brightness temperature $\delta T_{\rm b}$ fluctuations from reionization promise to provide information on the physical processes during that epoch. We present a formalism for generating the $\delta T_{\rm b}$ distribution using dark matter simulations and an one-dimensional radiative transfer code. Our analysis is able to account for the spin temperature $T_{\rm S}$ fluctuations arising from inhomogeneous X-ray heating and Ly$\alpha$ coupling during cosmic dawn. The $\delta T_{\rm b}$ power spectrum amplitude at large scales ($k \sim 0.1$ Mpc$^{-1}$) is maximum when $\sim 10\%$ of the gas (by volume) is heated above the CMB 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 $\lesssim 0.2$) this is because the signal is dominated by fluctuations in $T_{\rm S}$. For reionization models that are solely driven by stars within high mass ($\gtrsim 10^9\, \rm M_{\odot}$) haloes, the peculiar velocity effects are prominent only at smaller scales ($k \gtrsim 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.