Epoch of Reionization

The removal of the Galactic and extragalactic foregrounds remains a major challenge for those wishing to make a detection of the Epoch of Reionization 21-cm signal. Multiple methods of modelling these foregrounds with varying levels of assumption have been trialled and shown promising recoveries on simulated data. Recently however there has been increased discussion of using the expected shape of the foregrounds in Fourier space to define an EoR window free of foreground contamination. By carrying out analysis within this window only, one can avoid the foregrounds and any statistical bias they might introduce. In this paper we discuss the advantages and disadvantages of both foreground removal and foreground avoidance. We create a series of simulations with noise levels in line with both current and future experiments and compare the recovered statistical cosmological signal from both methods. We also relax the assumption that the foregrounds are smooth by introducing a Gaussian random factor along the line-of-sight and then also spatially. We find that, while foreground avoidance enables a slightly better recovery at high kperp, the most significant part of the cosmological signal is lost. While the EoR window is almost entirely obscured for a non-smooth foreground simulation, GMCA is able to recover the signal almost as well as for a smooth foreground for a line-of-sight variation.

After the Big Bang, the intergalactic medium consisted of a dense plasma of fully ionized baryonic materials and it was opaque. While the temperature of the universe was decreasing over time, the ionized particles formed hydrogen and other basic elements such as helium (He) and lithium (Li); however, they were low in concentration. During this time, the universe became transparent and was plunged into darkness. Then, the epoch of reionization arrived and newly-formed stars and galaxies lightened this dark universe. These primitive formations were the ancestors of today's galaxies, black holes, and, of importance for us, of the Earth. According to the current findings, this was the final major era of change in our today's universe. In addition, because of the formative diversity that occurred in this era, this era can be regarded as the rebirth of the universe. In this paper, we will derive the equations that describe this era and construct a theoretical framework for this topic. After this initial stage, the experimentation and computation methodologies to be employed will be discussed. After giving insight into the theory and experimentation of this era, we will try to ascertain why this topic is important as an area of study and discuss the expected future of this field. Most of the identities we will derive are already discussed in textbooks, but I will refine my perspective and understanding through re-deriving the equations of the theory and explaining the current experiments conducted in this field. This paper appeals to a very large audience, from leading scholars to students who show an interest in this topic. Having an intermediate-level calculus background would improve the experience of reading this paper and it is highly recommended for high-school students who wish to read this paper to acquire such prior knowledge. *

The visibility of Lyα-emitting galaxies during the Epoch of Reionization (EoR) is controlled by both diffuse H I patches in large-scale bubble morphology and small-scale absorbers. To investigate their impact on Lyα photons, we apply a powerful and novel combination of analytic and numerical calculations to three scenarios: (i) the 'bubble' model, where only diffuse H I outside ionized bubbles is present; (ii) the 'web' model, where H I exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there exists a degeneracy between the ionization structure of the inter-galactic medium (IGM) and the global H I fraction inferred from Lyα surveys, as the three models suppress Lyα flux equally with (very) different neutral fractions. We argue that a joint analysis of the Lyα luminosity function and the rest-frame equivalent width distribution/Lyα fraction of Lyα-and UV-selected galaxies can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M U V-dependent redshift evolution of the Lyα fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of Lyα fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. Our analysis further highlights the importance of galaxy-absorber cross-correlation and its systematic measurement from the post-reionized era to the EoR. Comparing our grid of models to observations, the neutral fraction at z ∼ 7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit on the neutral fraction of ∼ 1% in the web models, if we allow for a drop in the photoionization rate by a factor of ∼ 100 from the post-reionized universe.

The 21-cm power spectrum (PS) has been shown to be a powerful discriminant of reionization and cosmic dawn astrophysical parameters. However, the 21-cm tomographic signal is highly non-Gaussian. Therefore there is additional information which is wasted if only the PS is used for parameter recovery. Here we showcase astrophysical parameter recovery directly from 21-cm images, using deep learning with convolutional neural networks (CNN). Using a data base of 2D images taken from 10 000 21-cm light-cones (each generated from different cosmological initial conditions), we show that a CNN is able to recover parameters describing the first galaxies: (i) Tvir , their minimum host halo virial temperatures (or masses) capable of hosting efficient star formation; (ii) ζ ⁠, their typical ionizing efficiencies; (iii) LX/SFR , their typical soft-band X-ray luminosity to star formation rate; and (iv) E0 , the minimum X-ray energy capable of escaping the galaxy into the IGM. For most of their allowed ranges, log Tvir and log LX/SFR are recovered with <1 per cent uncertainty, while ζ and E0 are recovered with ∼10 per cent uncertainty. Our results are roughly comparable to the accuracy obtained from Monte Carlo Markov Chain sampling of the PS with 21CMMC for the two mock observations analysed previously, although we caution that we do not yet include noise and foreground contaminants in this proof-of-concept study.

We present an all-sky map synthesized between 139 MHz and 174 MHz using data from both arrays that reaches down to 80 mJy (4.9 K, for a beam size of 2.15e-5 steradians at 154 MHz), with a 10 mJy (620 mK) thermal noise level that indicates what would be achievable with better foreground subtraction. We calculate angular power spectra ($C_\ell$) in a cold patch and determine them to be dominated by point sources, but with contributions from galactic synchrotron emission at lower radio frequencies and angular wavemodes. Although the cosmic variance of foregrounds dominates errors in these power spectra, we measure a thermal noise level of 310 mK at $\ell=100$ for a 1.46-MHz band centered at 164.5 MHz. This sensitivity level is approximately three orders of magnitude in temperature above the level of the fluctuations in 21cm emission associated with reionization.

The angular power spectrum of the cosmic microwave background temperature anisotropy observed by WMAP has an anomalous dip at l~20 and bump at l~40. One explanation for this structure is the presence of features in the primordial curvature power spectrum, possibly caused by a step in the inflationary potential. The detection of these features is only marginally significant from temperature data alone. However, the inflationary feature hypothesis predicts a specific shape for the E-mode polarization power spectrum with a structure similar to that observed in temperature at l~20-40. Measurement of the CMB polarization on few-degree scales can therefore be used as a consistency check of the hypothesis. The Planck satellite has the statistical sensitivity to confirm or rule out the model that best fits the temperature features with 3 sigma significance, assuming all other parameters are known. With a cosmic variance limited experiment, this significance improves to 8 sigma. For tests of inflationary models that can explain both the dip and bump in temperature, the primary source of uncertainty is confusion with polarization features created by a complex reionization history, which at most reduces the significance to 2.5 sigma for Planck and 5-6 sigma for an ideal experiment. Smoothing of the polarization spectrum by a large tensor component only slightly reduces the ability of polarization to test for inflationary features, as does requiring that polarization is consistent with the observed temperature spectrum given the expected low level of TE correlation on few-degree scales. A future polarization satellite would enable a decisive test of the feature hypothesis and provide complementary information about the shape of a possible step in the inflationary potential. (Abridged.)

We examine the possibility of constraining the cosmological mean mass and dark energy densities by an application of the Alcock-Paczynski test on redshifted 21-cm maps of the epoch of reionization. The 21-cm data will be provided as a function of frequency and angular positions on the sky. The ratio of the frequency to angular distance scales can be determined by inspecting the anisotropy pattern of the correlation function of brightness temperature. We assess the sensitivity of the distance ratio to the cosmological parameters and present a technique for disentangling geometrical distortions from redshift space distortions caused by peculiar motions.

Next generation telescopes such as the James Webb Space Telescope (JWST) and the Wide Field Infrared Survey Telescope (WFIRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the astraeus (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above 100% for M UV ∼ −17.5 (−20) at z = 12 (6) for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin 2); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at z > 9, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields) of Lyman Break Galaxies during the EoR. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas.

The detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract H I power spectra that are applicable to tracking observations using an imaging radio interferometer ("Delay Spectrum with Imaging Arrays"). Our method is based on modeling the H I signal taking into account the impact of wide field effects such as the w-term, which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as the Murchison Widefield Array (MWA), Low Frequency Array (LOFAR), Giant Meterwave Radio Telescope (GMRT), Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays. This is the first application of delay spectrum methodology to such data to extract the H I signal. We analyze 3 hr of MWA tracking data on the EoR1 field. We present both two-dimensional ( ^ k k ,) and one-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyze the MWA EoR data.

We study the constraints on reionization from five years of WMAP data, parametrizing the evolution of the average fraction of ionized hydrogen with principal components that provide a complete basis for describing the effects of reionization on large-scale E-mode polarization. Using Markov Chain Monte Carlo methods, we find that the resulting model-independent estimate of the total optical depth is nearly twice as well determined as the estimate from 3-year WMAP data, in agreement with simpler analyses that assume instantaneous reionization. The mean value of the optical depth from principal components is slightly larger than the instantaneous value; we find tau=0.097+/-0.017 using only large-scale polarization, and tau=0.101+/-0.019 when temperature data is included. Likewise, scale invariant n_s=1 spectra are no longer strongly disfavored by WMAP alone. Higher moments of the ionization history show less improvement in the 5-year data than the optical depth. By plotting the distribution of polarization power for models from the MCMC analysis, we show that extracting most of the remaining information about the shape of the reionization history from the CMB requires better measurements of E-mode polarization on scales of l~10-20. Conversely, the quadrupole and octopole polarization power is already predicted to better than cosmic variance given any allowed ionization history at z<30 so that more precise measurements will test the LambdaCDM paradigm.

The Universe's Epoch of Reionization can be studied using a number of observational probes that provide complementary or corroborating information. Each of these probes suffers from its own systematic and statistical uncertainties. It is therefore useful to consider the mutual information that these data sets contain. In this paper, we present a cross-correlation study between the kinetic Sunyaev–Zel'dovich effect – produced by the scattering of cosmic microwave background (CMB) photons off free electrons produced during the reionization process – and the cosmological 21 cm signal – which reflects the neutral hydrogen content of the Universe, as a function of redshift. The study is carried out using a simulated reionization history in 100 h−1 Mpc scale N-body simulations with radiative transfer. In essence, we find that the two probes anticorrelate. The significance of the anticorrelation signal depends on the extent of the reionization process, wherein extended histories result in a much stronger signal compared to instantaneous cases. Unfortunately, however, once the primary CMB fluctuations are included into our simulation they serve as a source of large correlated noise that renders the cross-correlation signal insignificant, regardless of the reionization scenario.

The cross-correlation between high redshift galaxies and 21 cm emission from the high redshift intergalactic medium (IGM) promises to be an excellent probe of the Epoch of Reionization (EoR). On large scales, the 21 cm and galaxy fields are anti-correlated during most of the reionization epoch. However, on scales smaller than the size of the H II regions around detectable galaxies, the two fields become roughly uncorrelated. Consequently, the 21 cm-galaxy cross power spectrum provides a tracer of bubble growth during reionization, with the signal turning over on progressively larger scales as reionization proceeds. The precise turnover scale depends on the minimum host mass of the detectable galaxies, and the galaxy selection technique. Measuring the turnover scale as a function of galaxy luminosity constrains the characteristic bubble size around galaxies of different luminosities. The cross spectrum becomes positive on small scales if ionizing photons fail to escape from low mass galaxies, and these galaxies are detectable longward of the hydrogen ionization edge, because in this case some identifiable galaxies lie outside of ionized regions. LOFAR can potentially measure the 21 cm-galaxy cross spectrum in conjunction with mild extensions to the existing Subaru survey for $z=6.6$ Lyman-alpha emitters, while the MWA is slightly less sensitive for detecting the cross spectrum. A futuristic galaxy survey covering a sizable fraction of the MWA field of view ($\sim 800$ deg$^2$) can probe the scale dependence of the cross spectrum, constraining the filling factor of H II regions at different redshifts during reionization, and providing other valuable constraints on reionization models.

Redshifted 21cm signal is a promising tool to investigate the state of intergalactic medium (IGM) in the Cosmic Dawn (CD) and Epoch of Reionization(EoR). In our previous work (Shimabukuro et al 2015), we studied the variance and skewness to give a clear interpretation of 21cm power spectrum and found that skewness is a good indicator of the epoch when X-ray heating becomes effective. Thus, the non-Gaussian feature of the spatial distribution of the 21cm signal is expected to be useful to investigate the astrophysical effects in the CD and EoR. In this paper, in order to investigate such a non-Gaussian feature in more detail, we focus on the bispectrum of the 21cm signal. It is expected that the 21cm brightness temperature bispectrum is produced by non-gaussianity due to the various astrophysical effects such as Wouthysen-Field (WF) effect, X-ray heating and reionization. We study the various properties of 21cm bispectrum such as scale dependence, shape dependence and redshift evolut...

We use deep HST/ACS observations to calculate the star formation history (SFH) of the Cetus dwarf spheroidal (dSph) galaxy. Our photometry reaches below the oldest main sequence turn-offs, which allows us to estimate the age and duration of the main episode of star formation in Cetus. This is well approximated by a single episode that peaked roughly 12+/-0.5 Gyr ago and lasted no longer than about 1.9+/-0.5 Gyr (FWHM). Our solution also suggests that essentially no stars formed in Cetus during the past 8 Gyrs. This makes Cetus&#39; SFH comparable to that of the oldest Milky Way dSphs. Given the current isolation of Cetus in the outer fringes of the Local Group, this implies that Cetus is a clear outlier in the morphology-Galactocentric distance relation that holds for the majority of Milky Way dwarf satellites. Our results also show that Cetus continued forming stars through z ~ 1, long after the Universe was reionized, and that there is no clear signature of the epoch of reionizati...

Redshifted 21cm signal is a promising tool to investigate the state of intergalactic medium (IGM) in the Cosmic Dawn (CD) and Epoch of Reionization(EoR). In our previous work (Shimabukuro et al 2015), we studied the variance and skewness to give a clear interpretation of 21cm power spectrum and found that skewness is a good indicator of the epoch when X-ray heating becomes effective. Thus, the non-Gaussian feature of the spatial distribution of the 21cm signal is expected to be useful to investigate the astrophysical effects in the CD and EoR. In this paper, in order to investigate such a non-Gaussian feature in more detail, we focus on the bispectrum of the 21cm signal. It is expected that the 21cm brightness temperature bispectrum is produced by non-gaussianity due to the various astrophysical effects such as Wouthysen-Field (WF) effect, X-ray heating and reionization. We study the various properties of 21cm bispectrum such as scale dependence, shape dependence and redshift evolution. And also we study the contribution from each component of 21cm bispectrum. We find that the contribution from each component has characteristic scale-dependent feature, and it is expected that we could obtain more detailed information on the IGM in the CD and EoR by using the 21cm bispectrum in the future experiments, combined with the power spectrum and skewness.