Abstract
This paper considers the impact of large scale biasing of the IGM on reionization. The two simplest but extreme scenarios for IGM biasing are: an unbiased IGM which has a constant density and an IGM with density equal to the collapsed matter density. In this work, the relationship between the IGM density and the collapsed matter density is defined through an IGM bias parameter. The two extreme scenarios of homogeneous and perfectly biased IGM are produced for two extreme values of this bias parameter. It is found that, for the same level of reionization (i.e., for same global neutral hydrogen fraction). one could get very different 21 cm brightness temperature distributions for different values of this bias parameter. These distributions could give an order of magnitude more or less power as compared to the uniform case. It is also found that there exists a critical value for the IGM bias parameter for which there could be a near washout of the structure in the 21 cm brightness temperature distribution (i.e., zero power or a nearly uniform 21 cm brightness temperature distribution). To address the problem, a new method of generating 21 cm brightness temperature maps is used. The method uses the results of n-body simulations and then employs ray tracing to obtain the 21 cm brightness temperature maps. Towards the end, a prescription for the IGM bias parameter is given. This is derived within the framework of the Press-Schechter theory.
1 Introduction
Epoch of Reionization (EoR) is considered as one of the most difficult challenges in modern cosmology (see e.g., Shapiro and Giroux, 1987; ; , ; ; Zaroubi, 2013; ). Analysis of the EoR is demanding on observational (see e.g., Tozzi et al., 2000; ; ; Parsons et al., 2010; Hazelton et al., 2013; Thyagarajan et al., 2015; ; Thyagarajan et al., 2016; Mertens et al., 2018; ; ) as well as theoretical (see e.g., Madau et al., 1997; ; ; ; Mesinger et al., 2011; ; Qin et al., 2018; La Plante and Ntampaka, 2019; Park et al., 2019; Liu and Shaw, 2020) fronts. The possibility of observing EoR in the radio window using the 21 cm radiation was suggested a long time ago (Wouthuysen, 1952; ; Sunyaev and Zeldovich, 1975). The 21 cm line emitted by the neutral hydrogen in the IGM gets redshifted to much larger wavelengths and a continuous distribution of neutral hydrogen gives rise to a continuum of emission in the radio window (for reviews, see ; ; Choudhury, 2009; Pritchard and Loeb, 2012).
This paper explores the impact of large scale biasing of the IGM on the 21 cm power spectrum from EoR. Many earlier works have considered a clumpy IGM during the EoR (see e.g., Kohler et al., 2007; Pawlik et al., 2009; McQuinn et al., 2011; Raičević and Theuns, 2011; ; ; Jeeson-Daniel et al., 2014; So et al., 2014; Sobacchi and Mesinger, 2014; Trombetti and Burigana, 2014; Mao et al., 2020). Such an IGM is granular or clustered on small scales and is considered on subgrid level. Its eventual impact on the EoR power spectrum is considered by some authors (see e.g., Mao et al., 2020). This paper addresses large scale clumping or large scale biasing of the IGM. By large scale biasing one means how the distribution of the IGM is in the whole simulation box. Many numerical approaches of simulating EoR assume IGM to be having uniform density far from halo boundaries (see e.g., Mellema et al., 2006; Thomas and Zaroubi, 2008; Thomas et al., 2009; ). While the early versions of semi-analytic methods determine the neutral hydrogen fraction, , based on excursion sets arguments and then assume that the hydrogen density follows the dark matter density (see e.g., ; Mesinger et al., 2011). Now, the hydrogen in the IGM could have a more complicated dependence on the local density of the massive haloes. So the regions of the universe that have more or less massive halos also have more or less density of the background IGM. The simplest way to consider this is through following expression:
In this work, α is termed as the IGM bias parameter and it notifies the biasing of the IGM with respect to the collapsed massive haloes1. For , the IGM is uniform which means that the uncollapsed gas is distributed uniformly in the Inter-Galactic regions. For , the IGM is biased in the way that its density is proportional to the collapsed matter density: . These are the two extremes the models suggested by Eq. 1 take. Note that as , introduction of alpha does not change the value of mean or global amount of neutral hydrogen. It only changes its distribution. It is seen in this work that the two extreme distributions, given by and , could have very different power spectra. As 21 cm power spectrum is one the prime observables of the EoR, this work suggests that the distribution of the IGM with respect to the collapsed matter should be studied in more depth in order to gain a proper understanding of the epoch.
The paper is organized as follows: The next section discusses the method of generating the 21 cm maps. The method is ideally suited to address the problem of impact of a biased IGM on EoR. The later section discusses the main results of the paper. A prescription for the IGM bias parameter within the Press-Schechter framework is then derived. This is followed by the conclusion of the work in the last section.
2 Method for Generating 21 cm Maps
The method used in this work (
Raut, 2019) borrows ideas of semi-analytic methods (
;
;
Mesinger et al., 2011) and extensively numerical methods (
Razoumov and Scott, 1999;
;
Nakamoto et al., 2001). The parameters that are borrowed from semi-analytic approach are:
(minimum mass of collapsed haloes),
(number of ionizing photons deposited in the IGM per collapsed baryon in stars) and
(number representing the ratio of ionization rate to recombination rate). This method has similarities with the hybrid approaches of
Trac and Cen (2007);
Lee et al. (2008);
Zheng et al. (2011), but it is computationally less demanding as it emphasizes upon the scales of interests (
). Upcoming radio telescopes like HERA (
)
2and SKA1-Low (
Koopmans et al., 2015)
3will be more sensitive to these modes (
Kulkarni et al., 2016). The fiducial values of reionization parameters used in this work are:
,
and
for
(denotes the mean or global value of neutral hydrogen fraction.) The details of the method used for generating 21 cm maps is
as follows:
• The first step is to do N-body simulations and compute overdensity for each grid cell. I have used Illustris-1 simulations () for this step. The cosmological parameters used are wmap-9 (): , , and . The box-size is and gridsize is . This paper uses a method that is a small modification to the one used in the earlier work. Instead of using excursions sets (Press and Schechter, 1974; ; ; Lacey and Cole, 1993)to find out the collapsed fraction , one uses the total halo distribution itself to find the extent of collapsed structure. The provided by above n-body simulations is about and so this approach is justified.
Second step assigns ionizing photons to each grid cell based on following prescription ().
• Here, is the escape fraction of the ionizing photons, is the star formation efficiency and is the number of ionizing photons produced per baryon in stars. is the number of ionizing photons produced in a grid cell which has collapse fraction .
• The next step is to split all the photons, , into a large number () of photon-rays. Each photon ray carries a weight indicating the amount of photons it carries. Initially most of the photon rays have equal weights. Each photon ray is assigned a random position and direction inside the grid cells. Then each photon ray is propagated in small steps, i.e. 0.05 of grid cell size. Note that the number of photon-rays originating in a grid cell is proportional to count of that cell.
• The next step is to implement reionization. At each step of propagation of photon-rays, one computes the total weights of the photon-rays arriving in the grid cells. This total is compared with the amount of neutral hydrogen present in the grid cells. The photon-ray weights and neutral hydrogen counts are decremented at every step accounting for reionization. It is assumed at this step that each photon-ray can ionize amount of neutral hydrogen equal to the photon-ray weight (i.e., 100% ionization efficiency or a very large cross section for ionization reaction).
• The next step implements recombination. After the photon-ray weights and neutral hydrogen content of each grid cell is revised, the number of neutral hydrogen atoms that were ionized in that particular step, , is computed. This number is computed for each grid cell separately. Before implementation of the next step of photon-ray propagation, number of neutral hydrogen atoms are added back to each grid cell to account for recombinations. Here and is the ratio of recombination rate to ionization rate.
• The final step is to terminate the photon-ray propagation and assign brightness temperature to each grid cell. The first part is implemented when the total of photon-ray weights falls below 0.1% of its initial value. For the second part, we use the well-known expression (Wouthuysen, 1952; ; Madau et al., 1997; ),
3 The Biased IGM Analysis
This section discusses the effect of biased IGM on 21 cm brightness temperature distribution. The biasing of the IGM is modelled through Eq. 1 in the next two subsections.
3.1 Modification of the 21 cm Brightness Temperature Distribution With the IGM Bias Parameter
The 21 cm maps are generated for three values of α: 0, 0.5, and 1. As suggested earlier, corresponds to an IGM that is distributed uniformly in the whole simulation box, while, corresponds to an IGM that is distributed proportional to the extent of the collapsed structures. This way of assigning IGM to each grid cell could be over simplified and actual distribution of the IGM could be more complicatedly biased with respect to the collapsed matter (see the next section). Here, value of one is used for the recombination parameter, (Sobacchi and Mesinger, 2014) and value of is used for the minimum halo mass (). The values were chosen in such a way so as to produce global neutral hydrogen fraction of 0.4, 0.6, and 0.75 for redshifts 7, 8, and 9 respectively. This set of values of neutral hydrogen fraction is consistent with the CMBR and LyAlpha data (Kulkarni et al., 2016; Raut and Choudhury, 2018). To get the values for each redshift, we solve the following equation ( is the Helium fraction):
In this work, the global neutral hydrogen fraction, denoted by , is considered as a parameter instead of . Once and values are fixed, above equation uniquely determines the value for a given .
The reionization data cubes are generated as discussed in the previous section. The 21 cm brightness temperature slices (midway through the simulation box) are shown in Figure 1. This is done for three values of the IGM bias parameter () and for three values of redshift ().
FIGURE 1
As seen from Figure 1, the topology of the 21 cm brightness temperature distribution is changed as the IGM bias parameter is changed. For , it is outside-in with the 21 cm bright regions are surrounding the 21 cm dark regions. While for , it is inside-out with the 21 cm bright regions are embedded inside the 21 cm dark regions. This is expected as gas in the voids will remain neutral and hence 21 cm bright for the uniform IGM or case (The hydrogen near collapsed structure will be ionized while the hydrogen away from collapsed structures will remain neutral.) On the other hand, the voids will be devoid of gas and hence 21 cm dark for case. As expected, for the intermediate value of , the distribution is also intermediate and is close to being uniform.
These conclusions are fortified if one analyzes the 21 cm brightness temperature probability distributions. As seen from Figure 2, the probability distributions are opposite for and case. For , there are more bright regions as compared to the dark ones and for , there are more dark regions are compared to the bright ones. This is because, voids, which occupy a larger volume, are brighter for case and darker for case. For the intermediate case of , the plots indicate that the distributions are more or less uniform as they have a single peak and most of the areas in the probability distribution curves are covered by a small range of values.
FIGURE 2
One can also study the impact of the variation of the IGM bias parameter on the power spectrum. To obtain the power spectrum, the standard procedure is followed. First, the 21 cm brightness temperature distribution () is Fourier transformed and dimensionless power spectrum is obtained using following equations:The spherically averaged dimensionless power spectrum is the obtained by averaging over all possible angles
Here is the cosine of the angle that wavevector makes with the LOS or parallel direction and ϕ is the azimuthal angle.
As seen from Figure 3, for , the power spectrum is nearly same for the two extreme scenarios of and . One should take note that although the power spectrum is almost the same, the topologies are opposite. For the other two redshifts, though, there is significant enhancement of power for the case of . This could be because of combination of larger neutral hydrogen fraction and more isolation of collapsed structures as compared to the case. These two, combined, are producing more contrast in distribution for the case as compared to the case.
FIGURE 3
There is also a notable fall in the power spectrum for the case for redshifts and . As mentioned before, this is because the intermediate value of the IGM bias parameter drives the distribution toward uniformity. This means that distributions of the remnant neutral hydrogen and the collapsed structures are almost perfectly inversely related to each other. Thus the product that appears in the expression is almost constant. This is giving a nearly constant brightness temperature distribution and very small power spectrum. This behaviour is generic and seems to appear for every redshift for some α (see the discussion in the next subsection).
3.2 The Critical Value of the IGM Bias Parameter and the Uniform Distribution
The severe degradation of the 21 cm brightness temperature power spectrum is one of the main results of this paper and hence some further analysis is justified. Figure 4 shows the dependence of the power for three typical modes for three redshifts as a function of α. As mentioned before, these three modes are representative of the scales that would be measured by upcoming and future experiments. As seen from the figure, behavior of all the modes is similar: the plots are U-shaped. The power goes to almost zero for each mode for some value of α. Notably, this value of α seems to depend on the redshift alone and not on the scales. This means that all the scales hit minima at the same value of α. This observation indicates that there is indeed a near washout of fluctuations in the entire simulation box.
FIGURE 4
A critical value of the IGM bias parameter that causes a near uniform brightness temperature distribution is possibly a characteristic of the reionization process. To check this all three reionization parameters are varied from their fiducial values. This is done for the case. As mentioned earlier, the fiducial parameter values are , and . Figure 5 shows the variation of the three k-modes as a function of the IGM bias parameter for six different combinations of reionization parameters. The central region of the plots indicate the changed parameter while the other two parameters are kept unchanged from their fiducial values. For all the cases, one can see that there is either a critical value of the IGM bias parameter or a range of values of the IGM bias parameter for which there is a low signal. As the signal is nearly zero, this work warrants a careful study of distribution of IGM with respect to the collapsed structures. It should also be noted that the topology of the 21 cm brightness temperature distribution would be opposite on two sides of the critical value of the IGM bias parameter.
FIGURE 5
4 Prescription for IGM Distribution Based on Press-Schechter Theory
The IGM bias parameter being a constant throughout the box is ideally not possible but it is the next step to the other two extreme scenarios of it being equal to one or zero. These two extreme scenarios form basis of the EoR simulations done by many earlier works (see e.g.,
Mellema et al., 2006;
Thomas and Zaroubi, 2008;
;
Thomas et al., 2009;
Mesinger et al., 2011;
). A possible next step to a constant value for the IGM bias parameter is an IGM bias parameter that is a function of the overdensity parameter. This could be predicted within the framework of the Press-Schecheter theory. The following assumptions are employed to obtain the IGM distribution within this framework:
• The results of the Press-Schechter theory are exact,
• Baryons follow dark matter,
• All the baryons that are not part of galaxies (i.e., massive collapsed haloes) are considered as being part of the IGM.4
Let,It is shown in the Supplementary Appendix that for a grid cell with overdensity δ and collapse fraction is given by:
Here is average or global value of the collapse fraction. Figure 6 shows the dependence of the IGM bias parameter obtained using above equation.
FIGURE 6
This prescription should be able to give a good zeroth level baryon distribution. Environmental effects, baryonic physics and other astrophysical effects would give rise to modulations to this distribution. Using this prescription for the IGM density, one can now obtain the 21 cm brightness temperature distribution and it is shown in Figure 7.
FIGURE 7
The under dense grid cells are more in number than the over dense grid cells and they dominate the overall behaviour. For , the average value of the IGM bias parameter is about 0.22 and is close to the critical value. As a result, the distribution seen is also nearly uniform. One can find the power spectrum of the 21 cm brightness temperature distribution for the Press-Schechter prescription and compare it to the two extreme scenarios of the uniform and the perfectly biased IGM. The results are shown in Figure 8. As seen from the figure, there is degradation of power for all three cases compared to the other two extreme scenarios. It can also be seen that for there is a severe loss of power for the large scales as the average value of the IGM bias parameter is close to the critical value.
FIGURE 8
5 Discussion and Conclusion
In this work we discussed how an IGM that is having varied bias with respect to the massive halo population could impact the 21 cm brightness temperature distribution. It is found that the distribution as well as the power spectrum are significantly affected with the difference in the IGM bias parameter α (Eq. 1), while, the global neutral hydrogen fraction remains constant. It is also found that for some critical value of the IGM bias parameter, the distribution becomes almost uniform causing a near washout of the 21 cm power spectrum. This critical value of the IGM bias parameter is different for different redshifts. Thus, it is essential to estimate the large scale distribution of the IGM. In this work, we also saw that this distribution could be estimated within the framework of the Press-Schechter theory. It was also observed that for , the framework yielded a nearly uniform distribution. The degradation in the structure of 21 cm brightness temperature distribution occurs because there are more ionizing photons in the regions where there is more neutral hydrogen. This situation can also occur for standard scenarios (i.e. or ) if one has an that is density dependent. In this work has been assumed to be globally constant. But, consider case of for which neutral hydrogen follows the collapsed matter. If one has larger for over dense regions and smaller for under dense regions then one can arrive at a situation in which the extent of neutral hydrogen in each grid cell is about the same. And so the 21 cm brightness temperature distribution is nearly uniform.
Many early large scale structure studies assume that the baryons follow dark matter or . It is found in many works that, based on theory (see e.g. ; Springel et al., 2001; ; Zentner et al., 2005; ; ; ; Moster et al., 2018) and observations (see e.g. ; ; ; McGaugh et al., 2010; ; ), the above assumption need not be valid and more work is needed to get a conclusion on the exact dependence. As suggested by these works, the problem of finding exact dependence of the hydrogen content on the massive collapsed structures is a difficult one. In the realistic scenarios, where baryons do not strictly follow the dark matter, there could be more structure in the 21 cm brightness temperature distribution. This would possibly give more power on small scales. Baryonic physics could alter the bias parameter but the average value over the entire box is likely to remain close to the one obtained using the assumption of baryons following dark matter.
There are possibly two future directions to the problem posed in this work: One could assume that the baryon fraction parameter, , to be dependent on halo mass and find the its correlation to the halo mass based on observations (see e.g., Padmanabhan et al., 2016) or simulations (see e.g., ; ; Pillepich et al., 2018; ). One could also take a more direct approach and use the high resolution simulations themselves to use the hydrogen content of various cosmic structures (Martizzi et al., 2019). A very high resolution simulation could directly yield the hydrogen density for the whole simulation box.
The author would like to thank the Illustris collaboration for their N-body simulations. The author would also like to thank anonymous referees whose comments helped in improving the quality of the manuscript.
Both, and , can be obtained as function of true overdensity δ and hence the IGM bias parameter α can be obtained as a function of overdensity of collapsed massive haloes .
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
The author confirms being the sole contributor of this work and has approved it for publication.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2021.633007/full#supplementary-material.
Footnotes
1.^Note that the bias parameter is defined with respect to collapsed matter density and not with respect to the underlying dark-matter density
3.^http://astronomers.skatelescope.org/
4.^The baryons in uncollapsed haloes would give rise to a clumpy IGM. This will modify the recombination parameter and would make it position dependent. This complication is not considered in this work.
References
1
AbelT.NormanM. L.MadauP. (1999). Photon-conserving radiative transfer around point sources in multidimensional numerical cosmology. Astrophys. J.523, 66–71. 10.1086/307739
2
BardeenJ. M.BondJ. R.KaiserN.SzalayA. S. (1986). The statistics of peaks of Gaussian random fields. Astrophys. J.304, 15. 10.1086/164143
3
BarkanaR.LoebA. (2005). Detecting the earliest galaxies through two new sources of 21 centimeter fluctuations. Astrophys. J.626, 1–11. 10.1086/429954
4
BarkanaR.LoebA. (2001). In the beginning: the first sources of light and the reionization of the universe. Phys. Rep.349, 125–238. 10.1016/S0370-1573(01)00019-9
5
BarryN.WilenskyM.TrottC. M.PindorB.BeardsleyA. P.HazeltonB. J.et al (2019). Improving the epoch of reionization power spectrum results from murchison widefield array season 1 observations. Astrophys. J.884, 16. 10.3847/1538-4357/ab40a8
6
BensonA. J.SugiyamaN.NusserA.LaceyC. G. (2006). The epoch of reionization. Mon. Not. R. Astron. Soc.369, 1055–1080. 10.1111/j.1365-2966.2006.10426.x
7
BernardiG.de BruynA. G.BrentjensM. A.CiardiB.HarkerG.JelićV.et al (2009). Foregrounds for observations of the cosmological 21 cm line. I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field. Astron. Astrophys.500, 965–979. 10.1051/0004-6361/200911627
8
BondJ. R.ColeS.EfstathiouG.KaiserN. (1991). Excursion set mass functions for hierarchical Gaussian fluctuations. Astrophys. J.379, 440. 10.1086/170520
9
CenR.OstrikerJ. P. (1999). Where are the baryons?Astrophys. J.514, 1–6. 10.1086/306949
10
ChanM. H. (2019). A universal constant for dark matter-baryon interplay. Sci. Rep.9, 3570. 10.1038/s41598-019-39717-x
11
ChoudhuryT. R. (2009). Analytical models of the intergalactic medium and reionization. Curr. Sci.97, 841.
12
ChoudhuryT. R.HaehneltM. G.ReganJ. (2009). Inside-out or outside-in: the topology of reionization in the photon-starved regime suggested by Lyman-alpha forest data. Mon. Not. R. Astron. Soc.394, 960–977. 10.1111/j.1365-2966.2008.14383.x
13
CiardiB.FerraraA.GovernatoF.JenkinsA. (2000). Inhomogeneous reionization of the intergalactic medium regulated by radiative and stellar feedbacks. Mon. Not. R. Astron. Soc.314, 611–629. 10.1046/j.1365-8711.2000.03365.x
14
CoorayA.ChangT. C.UnwinS.ZemcovM.CoffeyA.MorrisseyP.et al (2019). Cosmic dawn intensity mapper. Bull. Am. Astron. Soc. 51, 23.
15
CrainR. A.EkeV. R.FrenkC. S.JenkinsA.McCarthyI. G.NavarroJ. F.et al (2007). The baryon fraction of 𝜆-CDM haloes. Mon. Not. R. Astron. Soc.377, 41–49. 10.1111/j.1365-2966.2007.11598.x
16
CrainR. A.SchayeJ.BowerR. G.FurlongM.SchallerM.TheunsT.et al (2015). The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations. Mon. Not. Roy. Astron. Soc.450, 1937–1961. 10.1093/mnras/stv725
17
DaiX.BregmanJ. N.KochanekC. S.RasiaE. (2010). On the baryon fractions in clusters and groups of galaxies. Astrophys. J.719, 119–125. 10.1088/0004-637X/719/1/119
18
DeBoerD. R.ParsonsA. R.AguirreJ. E.AlexanderP.AliZ. S.BeardsleyA. P.et al (2017). Hydrogen epoch of reionization array (HERA). Pub. Astron. Soc. Pac.129, 045001. 10.1088/1538-3873/129/974/045001
19
DiemandJ.KuhlenM.MadauP. (2007). Formation and evolution of galaxy dark matter halos and their substructure. Astrophys. J.667, 859–877. 10.1086/520573
20
DiemerB.StevensA. R. H.LagosC. d. P.CaletteA. R.TacchellaS.HernquistL.et al (2019). Atomic and molecular gas in IllustrisTNG galaxies at low redshift. Mon. Not. R. Astron. Soc.487, 1529–1550. 10.1093/mnras/stz1323
21
DillonJ. S.ParsonsA. R. (2016). Redundant array configurations for 21 cm cosmology. Astrophys. J.826, 181. 10.3847/0004-637X/826/2/181
22
EckertK. D.KannappanS. J.StarkD. V.MoffettA. J.BerlindA. A.NorrisM. A. (2016). Resolve and Eco: the halo mass-dependent shape of galaxy stellar and baryonic mass functions. Astrophys. J.824, 124. 10.3847/0004-637X/824/2/124
23
EmbersonJ. D.ThomasR. M.AlvarezM. A. (2013). The opacity of the intergalactic medium during reionization: resolving small-scale structure. Astrophys. J.763, 146. 10.1088/0004-637X/763/2/146
24
EttoriS. (2003). Are we missing baryons in galaxy clusters?Mon. Not. R. Astron. Soc.344, L13–L16. 10.1046/j.1365-8711.2003.06810.x
25
FanX.CarilliC. L.KeatingB. (2006). Observational constraints on cosmic reionization. Ann. Rev. Astron. Astrophys.44, 415–462. 10.1146/annurev.astro.44.051905.092514
26
FieldG. B. (1959a). An attempt to observe neutral hydrogen between the galaxies. Astrophys. J.129, 525. 10.1086/146652
27
FieldG. B. (1959b). The spin temperature of intergalactic neutral hydrogen. Astrophys. J.129, 536. 10.1086/146653
28
FinlatorK.OhS. P.ÖzelF.DavéR. (2012). Gas clumping in self-consistent reionization models. Mon. Not. R. Astron. Soc.427, 2464–2479. 10.1111/j.1365-2966.2012.22114.x
29
FukugitaM.HoganC. J.PeeblesP. J. E. (1998). The cosmic baryon budget. Astrophys. J.503, 518–530. 10.1086/306025
30
FurlanettoS. R.LidzA.LoebA.McQuinnM.PritchardJ. R.ShapiroP. R.et al (2009). “Cosmology from the highly-redshifted 21 cm line,” in Astro2010: the Astronomy and astrophysics decadal survey (The National Academic Press), 2010, 82.
31
FurlanettoS. R.OhS. P.BriggsF. H. (2006). Cosmology at low frequencies: the 21 cm transition and the high-redshift Universe. Phys. Rep.433, 181–301. 10.1016/j.physrep.2006.08.002
32
FurlanettoS. R.ZaldarriagaM.HernquistL. (2004). The growth of H II regions during reionization. Astrophys. J.613, 1–15. 10.1086/423025
33
GenelS.VogelsbergerM.SpringelV.SijackiD.NelsonD.SnyderG.et al (2014). Introducing the Illustris project: the evolution of galaxy populations across cosmic time. Mon. Not. R. Astron. Soc.445, 175–200. 10.1093/mnras/stu1654
34
GharaR.ChoudhuryT. R.DattaK. K. (2015). 21 cm signal from cosmic dawn: imprints of spin temperature fluctuations and peculiar velocities. Mon. Not. R. Astron. Soc.447, 1806–1825. 10.1093/mnras/stu2512
35
GiocoliC.PieriL.TormenG. (2008). Analytical approach to subhalo population in dark matter haloes. Mon. Not. R. Astron. Soc.387, 689–697. 10.1111/j.1365-2966.2008.13283.x
36
GiriS. K.MellemaG.DixonK. L.IlievI. T. (2018). Bubble size statistics during reionization from 21-cm tomography. Mon. Not. Roy. Astron. Soc.473, 2949–2964. 10.1093/mnras/stx2539
37
GnedinN. Y.OstrikerJ. P. (1997). Reionization of the universe and the early production of metals. Astrophys. J.486, 581–598. 10.1086/304548
38
GreigB.MesingerA.KoopmansL. V. E. (2020). Reionization and cosmic dawn astrophysics from the square kilometre array: impact of observing strategies. Mon. Not. R. Astron. Soc.491, 1398–1407. 10.1093/mnras/stz3138
39
HarkerG.ZaroubiS.BernardiG.BrentjensM. A.de BruynA. G.CiardiB.et al (2010). Power spectrum extraction for redshifted 21-cm Epoch of Reionization experiments: the LOFAR case. Mon. Not. R. Astron. Soc.405, 2492–2504. 10.1111/j.1365-2966.2010.16628.x
40
HeP.FengL. L.FangL. Z. (2005). Distributions of the baryon fraction on large scales in the universe. Astrophys. J.623, 601–611. 10.1086/428708
41
HinshawG.LarsonD.KomatsuE.SpergelD. N.BennettC. L.DunkleyJ.et al (2013). Nine-year wilkinson microwave anisotropy probe (WMAP) observations: cosmological parameter results. Astrophys. J.208, 19. 10.1088/0067-0049/208/2/19
42
Jeeson-DanielA.CiardiB.GrazianiL. (2014). Clumping factors of H II, He II and He III. Mon. Not. R. Astron. Soc.443, 2722–2732. 10.1093/mnras/stu1365
43
HazeltonB. J.MoralesM. F.SullivanI. S. (2013). The fundamental multi-baseline mode-mixing foreground in 21 cm epoch of reionization observations. Astrophys. J.770, 15. 10.1088/0004-637X/770/2/156
44
KohlerK.GnedinN. Y.HamiltonA. J. S. (2007). Large-scale simulations of reionization. Astrophys. J.657, 15–29. 10.1086/509907
45
KoopmansL.PritchardJ.MellemaG.AguirreJ.AhnK.BarkanaR.et al (2015). “The cosmic dawn and epoch of reionisation with SKA,” in Advancing Astrophysics with the Square Kilometre Array (AASKA14), Giardini Naxos, Italy, June 9–13, 2014.
46
KulkarniG.ChoudhuryT. R.PuchweinE.HaehneltM. G. (2016). Models of the cosmological 21 cm signal from the epoch of reionization calibrated with Ly α and CMB data. Mon. Not. R. Astron. Soc.463, 2583–2599. 10.1093/mnras/stw2168
47
La PlanteP.NtampakaM. (2019). Machine learning applied to the reionization history of the universe in the 21 cm signal. Astrophys. J.880, 110. 10.3847/1538-4357/ab2983
48
LaceyC.ColeS. (1993). Merger rates in hierarchical models of galaxy formation. Mon. Not. R. Astron. Soc.262, 627–649. 10.1093/mnras/262.3.627
49
LeeK. G.CenR.GottI.RichardJ.TracH. (2008). The topology of cosmological reionization. Astrophys. J.675, 8–15. 10.1086/525520
50
LiuA.ShawJ. R. (2020). Data analysis for precision 21 cm cosmology. Pub. Astron. Soc. Pac.132, 062001. 10.1088/1538-3873/ab5bfd
51
MadauP.MeiksinA.ReesM. J. (1997). 21 centimeter tomography of the intergalactic medium at high redshift. Astrophys. J.475, 429–444. 10.1086/303549
52
MaoY.KodaJ.ShapiroP. R.IlievI. T.MellemaG.ParkH.et al (2020). The impact of inhomogeneous subgrid clumping on cosmic reionization. Mon. Not. R. Astron. Soc.491, 1600–1621. 10.1093/mnras/stz2986
53
MartizziD.VogelsbergerM.ArtaleM. C.HaiderM.TorreyP.MarinacciF.et al (2019). Baryons in the Cosmic Web of IllustrisTNG—I: gas in knots, filaments, sheets, and voids. Mon. Not. R. Astron. Soc.486, 3766–3787. 10.1093/mnras/stz1106
54
McGaughS. S.SchombertJ. M.de BlokW. J. G.ZagurskyM. J. (2010). The baryon content of cosmic structures. Astrophys. J. Lett.708, L14–L17. 10.1088/2041-8205/708/1/L14
55
McQuinnM.OhS. P.Faucher-GiguèreC.-A. (2011). On lyman-limit Systems and the evolution of the intergalactic ionizing background. Astrophys. J.743, 82. 10.1088/0004-637X/743/1/82
56
MellemaG.IlievI. T.AlvarezM. A.ShapiroP. R. (2006). C2-ray: a new method for photon-conserving transport of ionizing radiation. New Astron.11, 374–395. 10.1016/j.newast.2005.09.004
57
MertensF.GhoshA.KoopmansL. V. E. (2018). “Robust foregrounds removal for 21-cm experiments,” in Peering towards cosmic dawn. Editors JelićV.van der HulstT. (Groningen, The Netherlands: Cambridge University Press), 333, 284–287. 10.1017/S1743921318000546
58
MesingerA.FurlanettoS.CenR. (2011). 21cmFast: a fast, seminumerical simulation of the high-redshift 21-cm signal. Mon. Not. R. Astron. Soc.411, 955–972. 10.1111/j.1365-2966.2010.17731.x
59
MosterB. P.NaabT.WhiteS. D. M. (2018). Emerge–an empirical model for the formation of galaxies since z 10. Mon. Not. R. Astron. Soc.477, 1822–1852. 10.1093/mnras/sty655
60
NakamotoT.UmemuraM.SusaH. (2001). The effects of radiative transfer on the reionization of an inhomogeneous universe. Mon. Not. R. Astron. Soc.321, 593–604. 10.1046/j.1365-8711.2001.04008.x
61
PadmanabhanH.ChoudhuryT. R.RefregierA. (2016). Modelling the cosmic neutral hydrogen from DLAs and 21-cm observations. Mon. Not. R. Astron. Soc.458, 781–788. 10.1093/mnras/stw353
62
ParkJ.MesingerA.GreigB.GilletN. (2019). Inferring the astrophysics of reionization and cosmic dawn from galaxy luminosity functions and the 21-cm signal. Mon. Not. R. Astron. Soc.484, 933–949. 10.1093/mnras/stz032
63
ParsonsA. R.BackerD. C.FosterG. S.WrightM. C. H.BradleyR. F.GugliucciN. E.et al (2010). The precision array for probing the epoch of Re-ionization: eight station results. Astron. J.139, 1468–1480. 10.1088/0004-6256/139/4/1468
64
PawlikA. H.SchayeJ.van ScherpenzeelE. (2009). Keeping the Universe ionized: photoheating and the clumping factor of the high-redshift intergalactic medium. Mon. Not. R. Astron. Soc.394, 1812–1824. 10.1111/j.1365-2966.2009.14486.x
65
PillepichA.NelsonD.HernquistL.SpringelV.PakmorR.TorreyP.et al (2018). First results from the IllustrisTNG simulations: the stellar mass content of groups and clusters of galaxies. Mon. Not. R. Astron. Soc.475, 648–675. 10.1093/mnras/stx3112
66
PressW. H.SchechterP. (1974). Formation of galaxies and clusters of galaxies by self-similar gravitational condensation. Astrophys. J.187, 425–438. 10.1086/152650
67
PritchardJ. R.LoebA. (2012). 21 cm cosmology in the 21st century. Rep. Prog. Phys.75, 086901. 10.1088/0034-4885/75/8/086901
68
QinY.DuffyA. R.MutchS. J.PooleG. B.GeilP. M.MesingerA.et al (2018). Dark-ages Reionization and Galaxy Formation Simulation—XIV. Gas accretion, cooling, and star formation in dwarf galaxies at high redshift. Mon. Not. R. Astron. Soc.477, 1318–1335. 10.1093/mnras/sty767
69
RaičevićM.TheunsT. (2011). Modelling recombinations during cosmological reionization. Mon. Not. R. Astron. Soc.412, L16–L19. 10.1111/j.1745-3933.2010.00993.x
70
RautD. (2019). A hybrid approach toward simulating reionization: coupling ray tracing with excursion sets. Astrophys. J.887, 81. 10.3847/1538-4357/ab5054
71
RautD.ChoudhuryT. R. (2018). Unbiased constraints on reionization model parameters in presence of the foreground wedge. arXiv e-prints arXiv:1806.08687.
72
RazoumovA. O.ScottD. (1999). Three-dimensional numerical cosmological radiative transfer in an inhomogeneous medium. Mon. Not. R. Astron. Soc.309, 287–298. 10.1046/j.1365-8711.1999.02775.x
73
ShapiroP. R.GirouxM. L. (1987). Cosmological H II regions and the photoionization of the intergalactic medium. Astrophys. J.321, L107–L112. 10.1086/185015
74
SoG. C.NormanM. L.ReynoldsD. R.WiseJ. H. (2014). Fully coupled simulation of cosmic reionization. II. Recombinations, clumping factors, and the photon budget for reionization. Astrophys. J.789, 149. 10.1088/0004-637X/789/2/149
75
SobacchiE.MesingerA. (2014). Inhomogeneous recombinations during cosmic reionization. Mon. Not. R. Astron. Soc.440, 1662–1673. 10.1093/mnras/stu377
76
SpringelV.WhiteS. D. M.TormenG.KauffmannG. (2001). Populating a cluster of galaxies - I. Results at [formmu2]z=0. Mon. Not. R. Astron. Soc.328, 726–750. 10.1046/j.1365-8711.2001.04912.x
77
SunyaevR. A.ZeldovichI. B. (1975). On the possibility of radioastronomical investigation of the birth of galaxies. Mon. Not. R. Astron. Soc.171, 375–379. 10.1093/mnras/171.2.375
78
ThomasR. M.ZaroubiS.CiardiB.PawlikA. H.LabropoulosP.JelićV.et al (2009). Fast large-scale reionization simulations. Mon. Not. R. Astron. Soc.393, 32–48. 10.1111/j.1365-2966.2008.14206.x
79
ThomasR. M.ZaroubiS. (2008). Time-evolution of ionization and heating around first stars and miniqsos. Mon. Not. R. Astron. Soc.384, 1080–1096. 10.1111/j.1365-2966.2007.12767.x
80
ThyagarajanN.JacobsD. C.BowmanJ. D.BarryN.BeardsleyA. P.BernardiG.et al (2015). Foregrounds in wide-field redshifted 21 cm power spectra. Astrophys. J.804, 14–15. 10.1088/0004-637X/804/1/14
81
ThyagarajanN.ParsonsA. R.DeBoerD. R.BowmanJ. D.Ewall-WiceA. M.NebenA. R.et al (2016). Effects of antenna beam chromaticity on redshifted 21 cm power spectrum and implications for hydrogen epoch of reionization array. Astrophys. J.825, 9. 10.3847/0004-637X/825/1/9
82
TozziP.MadauP.MeiksinA.ReesM. J. (2000). Radio signatures of H I at high redshift: mapping the end of the “dark ages”. Astrophys. J.528, 597–606. 10.1086/308196
83
TracH.CenR. (2007). Radiative transfer simulations of cosmic reionization. I. Methodology and initial results. Astrophys. J.671, 1–13. 10.1086/522566
84
TrombettiT.BuriganaC. (2014). Semi-analytical description of clumping factor and cosmic microwave background free-free distortions from reionization. Mon. Not. R. Astron. Soc.437, 2507–2520. 10.1093/mnras/stt2063
85
WouthuysenS. A. (1952). On the excitation mechanism of the 21-cm (radio-frequency) interstellar hydrogen emission line. Astron. J.57, 31–32. 10.1086/106661
86
ZaroubiS. (2013). The epoch of reionization. Astrophy. Space Sci. Lib.396, 45. 10.1007/978-3-642-32362-1˙2
87
ZentnerA. R.BerlindA. A.BullockJ. S.KravtsovA. V.WechslerR. H. (2005). The physics of galaxy clustering. I. A model for subhalo populations. Astrophys. J.624, 505–525. 10.1086/428898
88
ZhengZ.CenR.TracH.Miralda-EscudéJ. (2011). Radiative transfer modeling of lyα emitters. II. New effects on galaxy clustering. Astrophys. J.726, 38. 10.1088/0004-637X/726/1/38
Summary
Keywords
cosmology—theory, large scale structure, epoch of reionization, cosmology—galaxies, large scale structure formation
Citation
Raut D (2021) Homogeneous vs Biased IGM: Impact on Reionization. Front. Astron. Space Sci. 8:633007. doi: 10.3389/fspas.2021.633007
Received
24 November 2020
Accepted
20 January 2021
Published
29 March 2021
Volume
8 - 2021
Edited by
Carlo Burigana, Istituto di Radioastronomia di Bologna (INAF), Italy
Reviewed by
Emanuela Dimastrogiovanni, University of New South Wales, Australia
Kazuharu Bamba, Fukushima University, Japan
Updates
Copyright
© 2021 Raut.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dinesh Raut, dinesh.v.raut@gmail.com
This article was submitted to Cosmology, a section of the journal Frontiers in Astronomy and Space Sciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.