- Split View
-
Views
-
Cite
Cite
Michael Sitwell, Andrei Mesinger, Yin-Zhe Ma, Kris Sigurdson, The imprint of warm dark matter on the cosmological 21-cm signal, Monthly Notices of the Royal Astronomical Society, Volume 438, Issue 3, 01 March 2014, Pages 2664–2671, https://doi.org/10.1093/mnras/stt2392
- Share Icon Share
Abstract
We investigate the effects of warm dark matter (WDM) on the cosmic 21-cm signal. If dark matter exists as WDM instead of cold dark matter (CDM), its non-negligible velocities can inhibit the formation of low-mass haloes that normally form first in CDM models, therefore delaying star formation. The absence of early sources delays the build-up of UV and X-ray backgrounds that affect the 21-cm radiation signal produced by neutral hydrogen. With use of the 21CMFAST code, we demonstrate that the pre-reionization 21-cm signal can be changed significantly in WDM models with a free-streaming length equivalent to that of a thermal relic with mass mX of up to ∼10–20 keV. In such a WDM cosmology, the 21-cm signal traces the growth of more massive haloes, resulting in a delay of the 21-cm absorption signature and followed by accelerated X-ray heating. CDM models where astrophysical sources have a suppressed photon-production efficiency can delay the 21-cm signal as well, although its subsequent evolution is not as rapid as compared to WDM. This motivates using the gradient of the global 21-cm signal to differentiate between some CDM and WDM models. Finally, we show that the degeneracy between the astrophysics and mX can be broken with the 21-cm power spectrum, as WDM models should have a bias-induced excess of power on large scales. This boost in power should be detectable with current interferometers for models with mX ≲ 3 keV, while next-generation instruments will easily be able to measure this difference for all relevant WDM models.
INTRODUCTION
Hierarchical structure formation within the ΛCDM model has been exceptionally accurate in describing the large-scale Universe within the range ∼ 10 Mpc-1 Gpc, as demonstrated from studies of the cosmic microwave background (CMB) and the clustering of galaxies. However, for over a decade, concerns have been raised over whether the standard assumption of cold dark matter (CDM) provides an adequate fit to data on smaller, sub-Mpc scales. These include predictions from N-body simulations that yield an overabundance of galactic satellites in our galaxy and in the field (Klypin et al. 1999; Moore et al. 1999; Papastergis et al. 2011), as well as in voids (Peebles 2001), and produce overly dense galactic centres with ‘cuspy’ density profiles (de Blok et al. 2001; Donato et al. 2009; Newman et al. 2009) and are inconsistent with observations of the kinetic properties of bright Milky Way satellites (Boylan-Kolchin, Bullock & Kaplinghat 2011, 2012).
One possible explanation lies with baryonic feedback processes (Governato et al. 2007; Pontzen & Governato 2012; Garrison-Kimmel et al. 2013; Sobacchi & Mesinger 2013; Teyssier et al. 2013), although accurately modelling these mechanisms is often challenging and difficulties may persist in matching to observations.
Another possible explanation is to change the properties of dark matter so it is warm (warm dark matter – WDM).1 This may alleviate these small-scale problems due to the higher velocities of the dark matter. In this case, structures are smoothed on scales below the dark matter's free-streaming length. Non-relativistic residual velocities can delay halo collapse and star formation. These effects may reduce the number of sub-haloes and low-mass galaxies that are formed as well as flatten out galactic centres.
The two most popular WDM candidates in the literature motivated by particle physics have been the sterile neutrino (Dodelson & Widrow 1994; Abazajian, Fuller & Patel 2001; Boyarsky et al. 2009) and the gravitino (Bond, Szalay & Turner 1982; Pagels & Primack 1982). While WDM may be produced in a number of different ways, it is most often described as a thermal relic that decouples while relativistic, but is non-relativistic by matter-radiation equality as to preserve structure beyond the Mpc scale. In this case, the WDM would have a particle mass mX of the order of 1 keV. Although for our purposes, the free-streaming scale of the dark matter is a more fundamental quantity, we use the standard convention of discussing the WDM mass of a thermal relic instead. We caution that for other WDM-production mechanisms the correspondence between free-streaming length and mass will be different. We also remark that the results presented in this paper can be applicable to models other than WDM that have similar cutoff scales in their power spectrum (see, e.g. Cyr-Racine & Sigurdson 2013).
As WDM suppresses growth of small structures, which form first in the hierarchical structure formation of CDM, early star formation is delayed in WDM models. Detection of signals emitted from high-redshift objects either directly, such as from gamma-ray bursts (GRBs; Mesinger, Perna & Haiman 2005) or strongly lensed galaxies (Pacucci, Mesinger & Haiman 2013), or indirectly through the redshift of reionization (Barkana, Haiman & Ostriker 2001) can place constraints on mX. Recently, de Souza et al. (2013) using GRB catalogues placed a constraint of mX > 1.6 − 1.8 keV at 95 per cent CL. Requiring WDM models to be able to reproduce both the stellar mass function and Tully–Fisher relation places a lower bound of mX ≥ 0.75 keV (Kang, Macciò & Dutton 2013). The Lyman α forest can probe scales down to ∼ 1 Mpc and can provide strict limits on mX (Narayanan et al. 2000; Viel et al. 2005, 2008; Seljak et al. 2006), with the most recent and stringent constraint of mX > 3.3 keV at 2σ (Viel et al. 2013). Although it has been claimed that the less dense galactic cores formed in WDM models may provide a better fit to the kinematic data of bright Milky Way satellites (Lovell et al. 2012), there is an ongoing debate as to whether WDM with a mass above current lower bounds can create a large enough galactic core as needed to solve the ‘cusp–core’ problem (Villaescusa-Navarro & Dalal 2011; Macciò et al. 2012; though see de Vega, Falvella & Sanchez 2013).
Highly redshifted 21-cm radiation emitted from the hyperfine spin-flip of neutral hydrogen is a promising new tool to probe the high-redshift Universe (Madau, Meiksin & Rees 1997; Zaldarriaga, Furlanetto & Hernquist 2004; Furlanetto, Peng Oh & Briggs 2006; Morales & Wyithe 2010; Mesinger, Ewall-Wice & Hewitt 2013a). If WDM is present in sufficient quantities to significantly delay structure formation, it could potentially leave a trace within the 21-cm radiation signal. Light emitted by the first astrophysical sources can couple the spin temperature of neutral hydrogen to the kinetic temperature of the IGM through the Wouthuysen–Field (WF) mechanism (Wouthuysen 1952; Field 1958) as well as heat and ionize the IGM. Thus, a delay in the appearance of these early sources can alter the 21-cm signal and delay milestones in the signal. In this paper, we will examine the effects of WDM on the pre-reionization 21-cm signal. This era may be especially useful for examining WDM since WDM inhibits the formation of low-mass haloes that form first in CDM models and thus differences between the halo populations in CDM and WDM increase with redshift. As astrophysics is very poorly known at high redshifts (z ≥ 6), we will focus on characterizing degeneracies between the unknown astrophysics and the presence of WDM.
The outline of this paper is as follows. In Section 2, we review the effects of the free-streaming of the WDM on the linear power spectrum and its residual velocities on halo collapse. The basic properties of the 21-cm signal are outlined in Section 3 and its simulation is described in Section 4 with the simulation results discussed in Section 5. Throughout this paper, we assume cosmological parameter values of ΩΛ = 0.73, Ωm = 0.27, Ωb = 0.046, h = 0.7, σ8 = 0.82, ns = 0.96. We quote all quantities in comoving units, unless stated otherwise.
EFFECT OF WDM ON STRUCTURE FORMATION
Free-streaming
The free-streaming of WDM particles smears out perturbations on small scales, as WDM particles stream out of overdense regions and into underdense regions. Perturbations are suppressed on scales below that corresponding to the WDM particle horizon.
Residual velocities
In addition, the residual velocity dispersion of the WDM delays the growth of non-linear perturbations and consequently collapse into virialized haloes. This can be thought of as an ‘effective pressure’. BHO modelled the collapse in WDM by studying the collapse in an analogous system comprised of an adiabatic gas, so its root-mean-square velocity evolves as vrms∝1/a, as the case with WDM, and whose initial temperature is set such that it shares the same vrms with the WDM.
Halo abundances
The mean collapse fraction in CDM and WDM models can be seen in Fig. 1. At high redshifts, small haloes begin to collapse in CDM while no or few such haloes collapse in WDM, resulting in a large relative difference between the collapse fractions in these models. However, this difference becomes smaller with lower redshifts as objects on scales larger than that inhibited by WDM start to collapse in both models. At late times, in the CDM scenario, the mass within haloes of sizes suppressed by WDM only represents a small fraction of the total mass within all collapsed structures, so the relative difference between the mean collapse fraction in CDM and WDM models is small at those times. Therefore, while structure formation is delayed in WDM models, the mean collapse fraction rises more rapidly as compared to CDM.
COSMIC 21-cm SIGNAL
The earliest possible measurable cosmic 21-cm signal would be emitted during the ‘dark ages’ before significant star formation occurs. At these early times, the gas is dense enough so that collisional coupling is strong and TS ≈ TK. Before z ∼ 150, residual free electrons strongly couple the gas kinetic temperature to the CMB through Compton scattering, so TS ≈ TK ≈ Tγ and no 21-cm signal can be observed at this time. After this point, the remaining free electrons are so defuse that the gas is decoupled from the CMB and cools adiabatically as TK∝(1 + z)2. Since the CMB temperature decreases at the slower pace of Tγ ∝ (1 + z), a 21-cm signal in absorption may be observed (at least in principle) at this time (Bharadwaj & Ali 2004; Loeb & Zaldarriaga 2004; Naoz & Barkana 2005; Lewis & Challinor 2007). As the gas continues to cool, the collisional coupling becomes less efficient, driving TS back up to the CMB temperature. As this scenario is relatively unaffected by structure formation, we do not expect the presence of WDM to significantly affect this era of the 21-cm signal and will restrict our attention to later times with redshifts below z ∼ 35.2
It will be important to keep in mind that the kinetic temperature of the gas will be lower than the CMB temperature when WF coupling first becomes effective. As the Lyman α background grows, the increasing strength of the WF coupling will drive TS from a value near the CMB temperature to the lower kinetic temperature of the gas, thus producing another absorption signal. As WDM delays structure formation, the production of significant UV and X-ray backgrounds will be delayed, which in turn modifies the WF coupling, X-ray heating and reionization. We therefore focus our attention to the astrophysical epochs in the 21-cm signal.
SIMULATION OF 21-cm SIGNAL
The 21-cm signal is simulated using the publicly available 21 CMFAST code.3 This is a seminumerical simulation that generates density, velocity, ionization and spin-temperature fields in a 3D box with length size ∼Gpc. In this section, we briefly summarize the code. See Mesinger, Furlanetto & Cen (2011), Mesinger, Ferrara & Spiegel (2013b) and references within for further details.
An initial linear density field is generated as a Gaussian random field described by a power spectrum. The initial linear density field is then evolved using the Zel'Dovich approximation.
Ionization fields are generated by assuming that a region is ionized if it contains more ionizing photons than neutral hydrogen atoms (multiplied by |$1+\bar{n}_{{\rm rec}}$|, where |$\bar{n}_{{\rm rec}}$| is the mean number of recombinations per baryon). The excursion-set formalism is used with the condition that |$\zeta f_{{\rm coll}}({\boldsymbol x}, z, R) \ge 1 - x_{{\rm e}}({\boldsymbol x}, z, R)$| for a cell centred at location |${\boldsymbol x}$| to be fully ionized, where |$f_{{\rm coll}}({\boldsymbol x}, z, R)$| is the collapse fraction smoothed on scale R, ζ is the ionization efficiency and |$1-x_{{\rm e}}({\boldsymbol x}, z, R)$| is the remaining fraction of neutral hydrogen within R. This criterion is evaluated at deceasing scales R and if the cell is not marked as fully ionized as the scale of the pixel length is reached, the cell's ionization fraction is marked as |$\zeta f_{{\rm coll}}({\boldsymbol x}, z, R_{{\rm cell}}) + x_{{\rm e}}({\boldsymbol x}, z)$|. Lastly, we note that the ionization efficiency can be decomposed as |$\zeta =A_{{{\rm He}}}f_*f_{{\rm esc}}N_{{\rm ion}}/(1+\bar{n}_{{\rm rec}})$|, where fesc is the fraction of ionizing photons that escape their host galaxy, Nion is the number of ionizing photons per baryon inside stars and AHe is a correction factor due to the presence of Helium.
SIMULATION RESULTS
As much is unknown about astrophysical properties during high-redshift eras, we will examine possible degeneracies in the 21-cm signal between WDM and astrophysical quantities. As a first step, we will compare the delayed WDM 21-cm signal with that in CDM with a reduced photon-production efficiency. Specifically, we decrease the efficiency uniformly over frequency by decreasing f*, but note that f* is degenerate with other parameters used to calculate photon-production efficiencies.
The box used in our simulation runs was 750 Mpc on a side and comprised of 3003 cells. The 21-cm signal was simulated in the redshift range z = 5.6–35. We set the minimum halo virial temperature that supports star formation to be Tvir = 104 K as to approximate the minimum temperature needed to efficiently cool the halo gas through atomic cooling, neglecting possible feedback processes.4 Our fiducial model uses an f* value of f*fid = 10 per cent and an ionization efficiency ζ = 31.5.
Examples of the mean spin and kinetic temperatures for CDM and WDM models are plotted in Fig. 2. As expected, for WDM, TS stays near Tγ for a longer time and the lowest point in the absorption trough, where the X-ray heating rate first surpasses the adiabatic cooling rate, occurs later. As mentioned in Section 2.3, although the mean collapse fraction is lower in WDM models, it grows more rapidly, which is reflected in the heating of the gas. In addition, Fig. 2 shows curves for CDM with the lower f* value of f*/f*fid = 0.1, which in our model happens to delay star formation such that the minimum value of |$\bar{T}_{{\rm S}}$| occurs roughly at the same time as in the WDM example used. In this case, the X-ray heating rate increases at a much slower rate after the minimum in |$\bar{T}_{{\rm S}}$| as compared to the two other cases shown, since lowering f* reduces the photon-production efficiency in stars of all masses. In both non-fiducial cases shown, |$\bar{T}_{{\rm S}}$| and thus |$\delta \bar{ T}_{{\rm b}}$| reach a lower value in their absorption troughs since the gas undergoes further cooling in the extra time needed for the X-ray heating to become efficient.
The evolution of the mean brightness temperatures for WDM models with mX = 2, 3, 4 keV is shown in Fig. 3.5 It is readily seen that having WDM with a particle mass of a few keV can substantially change the mean 21-cm brightness temperature evolution. While lowering f* within CDM models can delay the strong absorption signal, the resulting absorption trough is much wider than in WDM. For the same delay in the minimum of |$\delta \bar{ T}_{{\rm b}}$|, the delay in reionization is greater for CDM than for WDM. Although reionization may be greatly delayed well past z = 6 in models with low values of f*, our primary focus is on the pre-reionization 21-cm signal. We caution against automatically discarding these models, as the star formation efficiency may diverge from earlier values by reionization.
Examining the gradient of the global signal in Fig. 3(b), we see that suppressing f* in CDM models only shifts the mean signal to lower redshifts. On the other hand, decreasing mX in WDM models increases the gradients of the mean signal. In CDM models, |$\mathrm{\partial} \delta \bar{ T}_{{\rm b}} / \mathrm{\partial} z$| attains values near 33 mK ( − 45 mK) near its maximum (minimum) regardless of its f* value. This can increase significantly in WDM models, for example to ∼ 64 mK ( ∼ − 77 mK) at its maximum (minimum) for WDM with mX = 2 keV.
The effect of WDM on the global 21-cm signal can be tracked through different ‘critical points’ in the signal's evolution. We choose these points to be the redshift zmin at which |$\delta \bar{ T}_{{\rm b}}$| reaches its minimum value, the redshift zh when the kinetic temperature of the gas is heated above the CMB temperature and the redshift of reionization zr taken to be the redshift where the mean ionized fraction is |$\bar{x}_i(z_{{\rm r}})=0.5$|. These points are plotted for both CDM and WDM in Fig. 4(a). The solid curves track the effect of lowering f* on the redshifts of the critical points in CDM models (the values of f* can be read from the upper horizontal axis). The dashed curves show the effect of WDM on these redshifts, where the value of mX for each model can be read from the lower horizontal axis.
We begin to explore possible degeneracies between CDM and WDM cosmologies by finding the value of f* required in CDM that would have a particular critical point occurring at the same redshift as it would in WDM with a particular value of mX. In other words, for a particular event that occurs at redshift ze, we would like to find the curve that satisfies ze(f*|CDM) = ze(mX|WDM). These curves for zmin, zh and zr can be seen in Fig. 4(b). We can see that if one uses the milestone zr to distinguish between CDM and WDM with mX = 2, 3, 4 keV, then f* has to be known within a factor of 3.0, 1.8 and 1.4, respectively. Using zmin instead, f* only has to be known within a factor of 50, 13 and 4.8 for mX = 2, 3, 4 keV, respectively, since the impact of WDM is larger at higher redshifts. Near mX = 15 keV, using zmin to distinguish WDM from CDM requires f* to be known within a factor of 1.1 and drops to 1.01 by mX ∼ 20 keV (although the astrophysical motivations for WDM as mentioned in the introduction loses much of its appeal past a few keV).
As the value of mX is lowered, the curves in Fig. 4(b) diverge from one another, as the more rapid growth of structure in WDM changes the relative timing of the milestones. Therefore, if f* is approximately constant throughout the epochs under consideration, adjusting the value of f* in CDM so that a particular critical point occurs at the same redshift as it does in WDM will misalign other critical points and thus cannot reproduce the whole history of |$\delta \bar{ T}_{{\rm b}}$| in WDM models.
However, we can mimic the WDM mean brightness temperature evolution with CDM if we allow f* to vary in time. To illustrate this, Fig. 5 shows the form of f*(z) needed to reproduce the mean 21-cm signal for WDM with mX = 2, 4 keV. At high redshifts (z ≳ 15, 25 for mX = 2, 4 keV), f* is more than an order of magnitude smaller than its value at the end of reionization to compensate for the delay of structure formation in WDM. When more massive haloes start to collapse (near z = 10, 20 for mX = 2, 4 keV), f* rises quickly by roughly an order of magnitude to mimic the more rapid change of the collapse fraction in WDM and finally levels off during reionization. While this evolution of f* may be possible, it seems contrived without an underlying model of such evolution.
Even in cases where f* evolves in such a way as to mimic the mean brightness temperature in WDM, one can differentiate between WDM and CDM by examining the spectrum of perturbations in the 21-cm signal at certain points in its evolution. Perturbations in the UV and X-ray fields add power to the 21-cm power spectrum |$\Delta _{21}^2$| on large scales. Since the bias of sources in WDM can be greater than that in CDM (Smith & Markovic 2011), more power is added on large scales in WDM than in CDM. This effect is most easily seen at times when inhomogeneities in xα or TK are at their maximum. Fig. 6 shows the evolution of the power spectrum for the modes k = 0.08 and 0.18 Mpc− 1, showing a three-peak structure, where the peaks from high to low redshift are associated with inhomogeneities in xα, TK and |$x_{{\rm H\,\small {I}}}$|, respectively. When inhomogeneities in TK are at their maximum, the power at k = 0.08, 0.18 Mpc− 1 can be boosted in WDM by as much as a factor of 2.4, 2.0 (1.3, 1.1) for mX = 2 keV (4 keV). When inhomogeneities in xα are near their height, the power at k = 0.08 Mpc− 1 can be increased by a factor of 1.5 (1.2) for WDM with mX = 2 keV (4 keV).
Current and next-generation interferometric radio telescopes may be used to detect the boost in power associated with WDM models. The dotted curves in Fig. 6 show forecasts for the 1σ power spectrum thermal noise levels for 2000 h of observation time, computed by Mesinger et al. (2013a), for the Murchison Widefield Array (MWA),6 the Square Kilometre Array (SKA)7 and for the proposed Hydrogen Epoch of Reionization Array (HERA).8 This estimate is quite conservative in that it ignores the contribution of foreground-contaminated modes (Pober et al. 2013). From these forecasts, we can see that the MWA may be able to at least marginally detect the boost in power for the mX = 2 keV model at the reionization and X-ray heating peaks. In addition, these estimates indicate that next-generation instruments will be able to easily measure the excess of power at these scales for mX = 2, 4 keV models over a wide range of redshifts.
The 21-cm power spectrum during a redshift near the time when TK is at its most inhomogeneous state is plotted in Fig. 7 for WDM with mX = 2, 4 keV and their CDM counterparts. One can see that the boost in power in WDM may continue to k values lower than those used in Fig. 6. In particular, the power near k = 0.01 Mpc− 1 in WDM models with mX = 2 keV (4 keV) may be larger by a factor of 3 (1.3) as compared to in CDM models at these times.
Finally, we mention that for simplicity we have chosen to vary only one astrophysical property. By allowing other astrophysical parameters to vary as a function of redshift, most notably Mmin, it might be possible to produce a 21-cm power spectrum degenerate with WDM throughout the redshifts under investigation and we leave this question for future work.
CONCLUSIONS
In WDM models, the abundance of small haloes is suppressed, which can leave a strong imprint at high redshifts. Since structure formation is delayed but more rapid in WDM, a delayed, deeper and more narrow absorption trough in the mean 21-cm signal will be produced in WDM models. These effects can easily be seen in the global 21-cm signal for WDM with free-streaming lengths above current observational bounds for thermal relic masses as high as mX ∼ 10–20 keV (|$R_{{\rm c}}^0\sim 6\hbox{--}13\,{\rm kpc}$|).
Suppressing the photon-production efficiency of astrophysical sources can delay the 21-cm signal as well. As such, to discriminate between WDM and CDM models by measuring the redshift of reionization, the photon-production efficiency must be known within a factor of 3.0, 1.8 and 1.4 for WDM with mX = 2, 3, 4 keV (|$R_{{\rm c}}^0\approx 86, 54, 39\,{\rm kpc}$|), respectively. Since the impact of WDM is larger at higher redshifts, if milestones in the mean 21-cm signal that occur at higher redshift are used to differentiate WDM and CDM models, the precision to which this efficiency must be known decreases. For example, if measuring the redshift of the minimum of the mean 21-cm signal (during the astrophysical epoch of the signal) the efficiency must only be known within a factor of 50, 13 and 4.8 for mX = 2, 3, 4 keV, respectively.
If the star formation remains approximately constant over the range of redshifts under consideration, degeneracy between CDM and WDM models may be broken by examining the gradient of the mean 21-cm signal, which is larger in WDM due to its more rapid pace of structure formation. In addition, the spectrum of perturbations in the 21-cm signal may as well be used to break this degeneracy, as the 21-cm power spectrum in WDM has an excess of power on large scales owing to the stronger biasing of sources in WDM. This is true even if the photon-production efficiency evolves with redshift in such a way as to reproduce with CDM the global 21-cm signal in WDM models. For WDM with mX = 2 keV (4 keV), the power in the 21-cm signal at k = 0.08, 0.18 Mpc− 1 can be increased by a factor as high as 2.4, 2.0 (1.3, 1.1) as compared to that in CDM. Power spectrum measurements made by current interferometric telescopes, such as the MWA, should be able to discriminate between CDM and WDM models with mX ≲ 3 keV, while next-generation telescopes will easily be able to differentiate between CDM and all relevant WDM models.
In this work, we assume that atomically cooled haloes drive the 21-cm signal. If instead smaller, molecularly cooled haloes, whose production is suppressed in WDM, play a significant role in producing the 21-cm signal in CDM, then the effects differentiating WDM from CDM described above would be even more pronounced. On the other hand, if star formation was not efficient in haloes with Tvir = 104 K, the differences between CDM and WDM in the 21-cm signal would be diminished.
This research was supported in part by the National Science and Engineering Research Council of Canada (NSERC). MS is supported in part by an NSERC Canada Graduate Scholarship. The research of KS is supported in part by an NSERC Discovery Grant. KS thanks the Aspen Center for Physics, where part of this work was completed, for their hospitality. YZM is supported by a CITA National Fellowship.
Other possible alterations to the standard CDM model that may resolve these small-scale problems include self-interacting dark matter (Burkert 2000; Spergel & Steinhardt 2000; Davé et al. 2001) and atomic dark matter or other models with acoustic damping of dark matter fluctuations (Kaplan et al. 2010; Cyr-Racine & Sigurdson 2013).
Although the very first stars were likely formed within smaller haloes with Tvir of the order of 103 K that were molecularly cooled, star formation in such haloes can easily be disrupted by feedback processes (Haiman, Abel & Rees 2000; Mesinger, Bryan & Haiman 2009), and we therefore neglect radiation from sources located in such haloes.
We caution the reader that WDM models with mX = 2, 3 keV are disfavoured by recent Lyman α observations (Viel et al. 2013). However, Lyman α forest constraints are still susceptible to astrophysical (thermal and ionization history) and observational (sky and continuum subtraction) degeneracies. Therefore, it is still useful to confirm these constraints using the redshifted 21-cm signal.