No Arabic abstract
Tomographic three-dimensional 21 cm images from the epoch of reionization contain a wealth of information about the reionization of the intergalactic medium by astrophysical sources. Conventional power spectrum analysis cannot exploit the full information in the 21 cm data because the 21 cm signal is highly non-Gaussian due to reionization patchiness. We perform a Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). We adopt a trained 3D Convolutional Neural Network (CNN) to compress the 3D image data into informative summaries (DELFI-3D CNN). We show that this method recovers accurate posterior distributions for the reionization parameters. Our approach outperforms earlier analysis based on two-dimensional 21 cm images. In contrast, an MCMC analysis of the 3D lightcone-based 21 cm power spectrum alone and using a standard explicit likelihood approximation results in inaccurate credible parameter regions both in terms of the location and shape of the contours. Our proof-of-concept study implies that the DELFI-3D CNN can effectively exploit more information in the 3D 21 cm images than a 2D CNN or power spectrum analysis. This technique can be readily extended to include realistic effects and is therefore a promising approach for the scientific interpretation of future 21 cm observation data.
We present a simulation-based inference framework using a convolutional neural network to infer dynamical masses of galaxy clusters from their observed 3D projected phase-space distribution, which consists of the projected galaxy positions in the sky and their line-of-sight velocities. By formulating the mass estimation problem within this simulation-based inference framework, we are able to quantify the uncertainties on the inferred masses in a straightforward and robust way. We generate a realistic mock catalogue emulating the Sloan Digital Sky Survey (SDSS) Legacy spectroscopic observations (the main galaxy sample) for redshifts $z lesssim 0.09$ and explicitly illustrate the challenges posed by interloper (non-member) galaxies for cluster mass estimation from actual observations. Our approach constitutes the first optimal machine learning-based exploitation of the information content of the full 3D projected phase-space distribution, including both the virialized and infall cluster regions, for the inference of dynamical cluster masses. We also present, for the first time, the application of a simulation-based inference machinery to obtain dynamical masses of around $800$ galaxy clusters found in the SDSS Legacy Survey, and show that the resulting mass estimates are consistent with mass measurements from the literature.
The 21-cm intensity mapping (IM) of neutral hydrogen (HI) is a promising tool to probe the large-scale structures. Sky maps of 21-cm intensities can be highly contaminated by different foregrounds, such as Galactic synchrotron radiation, free-free emission, extragalactic point sources, and atmospheric noise. We here present a model of foreground components and a method of removal, especially to quantify the potential of Five-hundred-meter Aperture Spherical radio Telescope (FAST) for measuring HI IM. We consider 1-year observational time with the survey area of $20,000,{rm deg}^{2}$ to capture significant variations of the foregrounds across both the sky position and angular scales relative to the HI signal. We first simulate the observational sky and then employ the Principal Component Analysis (PCA) foreground separation technique. We show that by including different foregrounds, thermal and $1/f$ noises, the value of the standard deviation between reconstructed 21-cm IM map and the input pure 21-cm signal is $Delta T = 0.034,{rm mK}$, which is well under control. The eigenmode-based analysis shows that the underlying HI eigenmode is just less than $1$ per cent level of the total sky components. By subtracting the PCA cleaned foreground+noise map from the total map, we show that PCA method can recover HI power spectra for FAST with high accuracy.
Next generation observatories will enable us to study the first billion years of our Universe in unprecedented detail. Foremost among these are 21-cm interferometry with the HERA and the SKA, and high-$z$ galaxy observations with the James Webb Space Telescope (JWST). Taking a basic galaxy model, in which we allow the star formation rates and ionizing escape fractions to have a power-law dependence on halo mass with an exponential turnover below some threshold, we quantify how observations from these instruments can be used to constrain the astrophysics of high-$z$ galaxies. For this purpose, we generate mock JWST LFs, based on two different hydrodynamical cosmological simulations; these have intrinsic luminosity functions (LFs) which turn over at different scales and yet are fully consistent with present-day observations. We also generate mock 21-cm power spectrum observations, using 1000h observations with SKA1 and a moderate foreground model. Using only JWST data, we predict up to a factor of 2-3 improvement (compared with HST) in the fractional uncertainty of the star formation rate to halo mass relation and the scales at which the LFs peak (i.e. turnover). Most parameters regulating the UV galaxy properties can be constrained at the level of $sim 10$% or better, if either (i) we are able to better characterize systematic lensing uncertainties than currently possible; or (ii) the intrinsic LFs peak at magnitudes brighter than $M_{rm UV} lesssim -13$. Otherwise, improvement over HST-based inference is modest. When combining with upcoming 21-cm observations, we are able to significantly mitigate degeneracies, and constrain all of our astrophysical parameters, even for our most pessimistic assumptions about upcoming JWST LFs. The 21-cm observations also result in an order of magnitude improvement in constraints on the EoR history.
We present here 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 Lyman-alpha and CMB data. We model AGN by extrapolating the observed M-sigma relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen 21 cm spin temperature is significantly larger than the CMB temperature. Our 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 mK^2 at scales accessible to experiments (k < 1 h/cMpc). This is lower than the sensitivity claimed to have been already reached by ongoing experiments by only a factor of about two or less. When reionization is dominated by AGN, the 21 cm power spectrum is enhanced on all scales due to the enhanced bias of the clustering of the more massive haloes and the peak in the large scale 21 cm power is strongly enhanced and moved to larger scales due to bigger characteristic bubble sizes. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity, assuming successful foreground subtraction and instrument calibration. Conversely, these could become the first non-trivial reionization scenarios to be ruled out by 21 cm experiments, thereby constraining the contribution of AGN to reionization.
The redshifted 21-cm signal of neutral Hydrogen is a promising probe into the period of evolution of our Universe when the first stars were formed (Cosmic Dawn), to the period where the entire Universe changed its state from being completely neutral to completely ionized (Reionization). The most striking feature of this line of neutral Hydrogen is that it can be observed across an entire frequency range as a sky-averaged continuous signature, or its fluctuations can be measured using an interferometer. However, the 21-cm signal is very faint and is dominated by a much brighter Galactic and extra-galactic foregrounds, making it an observational challenge. We have used different physical models to simulate various realizations of the 21-cm Global signals, including an excess radio background to match the amplitude of the EDGES 21-cm signal. First, we have used an artificial neural network (ANN) to extract the astrophysical parameters from these simulated datasets. Then, mock observations were generated by adding a physically motivated foreground model and an ANN was used to extract the astrophysical parameters from such data. The $R^2$ score of our predictions from the mock-observations is in the range of 0.65-0.89. We have used this ANN to predict the signal parameters giving the EDGES data as the input. We find that the reconstructed signal closely mimics the amplitude of the reported detection. The recovered parameters can be used to infer the physical state of the gas at high redshifts.