The author - with his collaborators - already in years 1995-96 have shown - purely from the analyses of the observations - that the gamma-ray bursts (GRBs) can be till redshift 20. Since that time several other statistical studies of the spatial distribution of GRBs were provided. Remarkable conclusions concerning the star-formation rate and the validity of the cosmological principle were obtained about the regions of the cosmic dawn. In this contribution these efforts are surveyed.

Cosmic reionization is thought to be initiated by the first generation of stars and blackholes. We review recent progress in theoretical studies of early structure formation. Cosmic structure formation is driven by gravitational instability of primeval density fluctuations left over from Big Bang. At early epochs, there are baryonic streaming motions with significant relative velocity with respect to dark matter. The formation of primordial gas clouds is typically delayed by the streaming motions, but then physical conditions for the so-called direct collapse blackhole formation are realized in proto-galactic halos. We present a promising model in which intermediate mass blackholes are formed as early as z = 30.

The evolution of neutral hydrogen (HI) across redshifts is a powerful probe of cosmology, large scale structure in the universe and the intergalactic medium. Using a data-driven halo model to describe the distribution of HI in the post-reionization universe (z ∼ 5 to 0), we obtain the best-fitting parameters from a rich sample of observational data: low redshift 21-cm emission line studies, intermediate redshift intensity mapping experiments, and higher redshift Damped Lyman Alpha (DLA) observations. Our model describes the abundance and clustering of neutral hydrogen across redshifts 0 - 5, and is useful for investigating different aspects of galaxy evolution and for comparison with hydrodynamical simulations. The framework can be applied for forecasting future observations with neutral hydrogen, and extended to the case of intensity mapping with molecular and other line transitions at intermediate redshifts.

The Large Synoptic Survey Telescope (LSST), a next generation astronomical survey, sited on Cerro Pachon in Chile, will provide an unprecedented amount of imaging data for studies of the faint optical sky. The LSST system includes an 8.4m (6.7m effective) primary mirror and a 3.2 Gigapixel camera with a 9.6 sq. deg. field of view. This system will enable about 10,000 sq. deg. of sky to be covered twice per night, every three to four nights on average, with typical 5-sigma depth for point sources of r = 24.5 (AB). With over 800 observations in the ugrizy bands over a 10-year period, these data will enable coadded images reaching r = 27.5 (about 5 magnitudes deeper than SDSS) as well as studies of faint time-domain astronomy. The measured properties of newly discovered and known astrometric and photometric transients will be publicly reported within 60 sec after closing the shutter. The resulting hundreds of petabytes of imaging data for about 40 billion objects will be used for scientific investigations ranging from the properties of near-Earth asteroids to characterizations of dark matter and dark energy. For example, simulations estimate that LSST will discover about 1,000 quasars at redshifts exceeding 7; this sample will place tight constraints on the cosmic environment at the end of the reionization epoch. In addition to a brief introduction to LSST, I review the value of LSST data in support of epoch of reionization experiments and discuss how international participants can join LSST.

In the last decade, it has become clear that the dust-enshrouded star formation contributes significantly to early galaxy evolution. Detection of dust is therefore essential in determining the properties of galaxies in the high-redshift universe. This requires observations at the (sub-)millimeter wavelengths. Unfortunately, sensitivity and background confusion of single dish observations on the one hand, and mapping efficiency of interferometers on the other hand, pose unique challenges to observers. One promising route to overcome these difficulties is intensity mapping of fluctuations which exploits the confusion-limited regime and measures the collective light emission from all sources, including unresolved faint galaxies. We discuss in this contribution how 2D and 3D intensity mapping can measure the dusty star formation at high redshift, through the Cosmic Infrared Background (2D) and [CII] fine structure transition (3D) anisotropies.

The large scatter in Lyman-α opacity at z > 5.3 has been an ongoing mystery, prompting a flurry of numerical models. A uniform ultra-violet background has been ruled out at those redshifts, but it is unclear whether any proposed models produce sufficient inhomogeneities. In this paper we provide an update on the measurement which first highlighted the issue: Lyman-α effective optical depth along high-z quasar lines of sight. We nearly triple on the previous sample size in such a study thanks to the cooperation of the DES-VHS, SHELLQs, and SDSS collaborations as well as new reductions and spectra. We find that a uniform UVB model is ruled out at 5.1 < z < 5.3, as well as higher redshifts, which is perplexing. We provide the first such measurements at z ∼ 6. None of the numerical models we confronted to this data could reproduce the observed scatter.

Euclid enables the exploration of large sky areas with diffraction limited resolution in the optical and near-infrared, and is sensitive enough to detect targets at cosmological distances. This combination of capabilities gives Euclid a clear advantage over telescope facilities with larger apertures, both on ground and in space. The decision to mount in the NISP instrument one extra grism for the wavelength range 0.92-1.3 μm with a spectral resolution of R ≈260 makes possible a rest-frame UV survey of the early Universe in the redshift range 6.5 < z < 9.7. Euclid’s standard imaging with VIS in the 0.55-0.9 μm band and with NISP in the Y, J, H bands provide complementary photometry for further target identification and characterization. Euclid is a suitable facility to discover and map the spatial distribution of rare high-redshift targets and to collect statistically relevant samples, in particular of high redshift Lyα emitters and QSOs, which can be used as signposts of the cosmic structures. The Euclid surveys are also a starting point for deeper follow up observations of the individual high-z objects. We present the Euclid mission and discuss the detectability of high-z objects to probe the epoch of ionization.

We study the inhomogeneity of the reionization process by comparing the reionization times of z = 0 galaxies as a function of their mass. For this purpose, we combine the results of the CODA-I AMR radiative hydrodynamics simulation of the Reionization with the halo merger trees of a pure dark matter tree-code z = 0 simulation evolved from the same set of initial conditions. We find that galaxies with M(z = 0) > 1011M⊙ are reionized earlier than the whole Universe, with e.g. MW-like haloes reionized between 100 and 300 million years before the diffuse IGM. Lighter galaxies reionized as late as the global volume, probably from external radiation.

We present predictions for the spatial distribution of 21 cm brightness temperature fluctuations from high-dynamic-range simulations for AGN-dominated reionization histories that have been tested against available Lyα and CMB data. We model AGN by extrapolating the observed Mbh–σ relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. AGN-dominated reionization histories increase the variance of the 21 cm emission by a factor of up to ten compared to similar reionization histories dominated by faint galaxies, to values close to 100 mK2 at scales accessible to experiments (k ≲ 1 cMpc−1h). This is lower than the sensitivity reached by ongoing experiments by only a factor of about two or less. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity.

Characterizing the properties and the evolution of the first stars and galaxies is a challenging task for traditional galaxy surveys since they are sensitivity limited and can only detect the brightest light sources. Three-dimensional intensity mapping (IM) of transition lines can be a valuable alternative to study the high redshift Universe given that this technique avoids sensitivity limitation problems by measuring the overall emission of a line, with a low resolution, without resolving its sources. While 21cm line IM surveys probe neutral hydrogen gas and can, therefore, be used to probe the state of the IGM and the evolution of the ionization field during the Epoch of Reionization (EoR). IM surveys of other lines, such as CO, CII, Ly-alpha or H-alpha, can be used to probe the galaxies which emitted most of the ionizing radiation responsible for the EoR. These lines will trace the different ISM gas phases, the excitation state of this gas, its metallicity, etc. This study addresses IM of multiple transition lines and how it can be used to probe the EoR and to constrain the redshift evolution of galaxy properties.

The James Webb Space Telescope (JWST) will observe the rest-frame ultraviolet/optical spectra of galaxies from the epoch of reionization (EoR) in unprecedented detail. While escaping into the intergalactic medium, hydrogen-ionizing (Lyman continuum; LyC) photons from the galaxies will contribute to the bluer end of the UV slope and make nebular emission lines less prominent. We present a method to constrain leakage of the LyC photons using the spectra of high redshift (z ≳ 6) galaxies. We simulate JWST/NIRSpec observations of galaxies at z =6–9 by matching the fluxes of galaxies observed in the Frontier Fields observations of galaxy cluster MACS-J0416. Our method predicts the escape fraction fesc with a mean absolute error Δfesc ≈ 0.14. The method also predicts the redshifts of the galaxies with an error $\left\langle \frac{\Delta z}{(1+z)}\right\rangle \approx 0.0003$.

The EoR foregrounds can be up to three magnitudes greater than the cosmological signal we wish to detect. Multiple methods have been developed in order to extract the cosmological signal, falling roughly into three categories: foreground removal, foreground suppression and foreground avoidance. These main approaches are briefly discussed in this review and consideration taken to the future application of these methods as a multi-layered approach.

21-cm observations of the Cosmic dawn (CD) and Epoch of Reionization (EoR) are one of the high priority science objectives for SKA Low. One of the most difficult aspects of the 21-cm measurement is the presence of foreground emission, due to our Galaxy and extragalactic sources, which is about four orders of magnitude brighter than the cosmological signal. While end-to-end simulations are being produced to investigate in details the foreground subtraction strategy, it is useful to complement this thorough but time-consuming approach with simpler, quicker ways to evaluate performance and identify possible critical steps. In this work, I present a forecast method, based on Bonaldi et al. (2015), Bonaldi & Ricciardi (2011), to understand the level of residual contamination after a component separation step, and its impact on our ability to investigate CD and EoR.

Bayesian component separation techniques have played a central role in the data reduction process of Planck. The most important strength of this approach is its global nature, in which a parametric and physical model is fitted to the data. Such physical modeling allows the user to constrain very general data models, and jointly probe cosmological, astrophysical and instrumental parameters. This approach also supports statistically robust goodness-of-fit tests in terms of data-minus-model residual maps, which are essential for identifying residual systematic effects in the data. The main challenges are high code complexity and computational cost. Whether or not these costs are justified for a given experiment depends on its final uncertainty budget. We therefore predict that the importance of Bayesian component separation techniques is likely to increase with time for intensity mapping experiments, similar to what has happened in the CMB field, as observational techniques mature, and their overall sensitivity improves.

Observations of HI 21cm transition line is a promising probe into the Dark Ages and Epoch-of-Reionization. Detection of this redshifted 21cm signal is one of the key science goal for several upcoming low-frequency radio telescopes like HERA, SKA and DARE. Other global signal experiments include EDGES, LEDA, BIGHORNS, SCI-HI, SARAS. One of the major challenges for the detection of this signal is the accuracy of the foreground source removal. Several novel techniques have been explored already to remove bright foregrounds from both interferometric as well as total power experiments. Here, we present preliminary results from our investigation on application of ANN to detect 21cm global signal amidst bright galactic foreground. Following the formalism of representing the global 21cm signal by ’tanh’ model, this study finds that the global 21cm signal parameters can be accurately determined even in the presence of bright foregrounds represented by 3rd order log-polynomial or higher.

Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emission and of its contaminants due to instrument chromaticity, ionosphere and imperfect calibration. To solve for this complex modeling, we propose a new method based on Gaussian Process Regression (GPR) which is able to cleanly separate the cosmological signal from most of the foregrounds contaminants. We also propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere. Using this method, chromatic effects causing the so-called “wedge” are effectively eliminated (i.e. deconvolved) in the cylindrical (k⊥, k∥) power spectrum.

Intensity mapping (IM) is a new observational technique to survey the large-scale structure of matter using spectral emission lines. IM observations are contaminated by instrumental noise and astrophysical foregrounds. The foregrounds are at least three orders of magnitude larger than the searched signals. In this work, we apply the Generalized Needlet Internal Linear Combination (GNILC) method to subtract radio foregrounds and to recover the cosmological HI and CO signals within the IM context. For the HI IM case, we find that GNILC can reconstruct the HI plus noise power spectra with 7.0% accuracy for z = 0.13 − 0.48 (960 − 1260 MHz) and ℓ ≲ 400, while for the CO IM case, we find that it can reconstruct the CO plus noise power spectra with 6.7% accuracy for z = 2.4 − 3.4 (26 − 34 GHz) and ℓ ≲ 3000.

The cross power spectrum of the 21 cm signal and Lyman-α emitters (LAEs) is a probe of the Epoch of Reionization. Astrophysical foregrounds do not correlate with the LAE distribution, though the foregrounds contribute to the error. To study the impact of foregrounds on the measurement, we assume realistic observation by the Murchison Widefield Array using a catalogue of radio galaxies, a LAE survey by the Subaru Hyper Supreme-Cam and the redshift of LAEs is determined by the Prime Focus Spectrograph. The HI distribution is estimated from a radiative transfer simulation with models based on results of radiation hydrodynamics simulation. Using these models, we found that the error of cross power spectrum is dominated by foreground terms. Furthermore, we estimate the effects of foreground removal, and find 99% of the foreground removal is required to detect the 21 cm-LAE signal at k ∼ 0.4 h Mpc−1.