Source author record

Saleem Zaroubi

Saleem Zaroubi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2025arXiv

Velocity fields and turbulence from cosmic filaments to galaxy clusters

Galaxy clusters are currently the endpoint of the hierarchical structure formation; they form via the accretion of dark matter and cosmic gas from their local environment. In particular, filaments contribute grandly by accreting gas from cosmic matter sheets and underdense regions and feeding it to the galaxy clusters. Along the way, the gas in filaments is shocked and heated, which, together with the velocity structure within the filament, induces swirling and, thus, turbulence. In this work, we study a constrained hydrodynamical simulation replica of the Virgo cluster to characterise the velocity field in the two cosmic filaments connected to the cluster with unprecedented high resolution. First, we conduct a qualitative examination of slices extracted from the simulation. We study the temperature, the velocity field, and derived quantities in longitudinal cuts to study the general structure of the filaments and in transverse cuts to study their inner organisation and connection to cosmic matter sheets and underdense regions. Then, we conduct a quantitative study of velocities in Virgo's filaments by computing the 2D energy spectrum from 1 and 5~Mpc square maps extracted from the slices and centred on the core of the filaments. We show that the velocity field goes from mostly compressive far in the filaments to mostly solenoidal in Virgo's core. Moreover, we observe that the total energy spectrum in the filaments gains in amplitude and steepens towards Virgo.

preprint2022arXiv

Astrophysical information from the Rayleigh-Jeans Tail of the CMB

One of the explanations for the recent EDGES-LOW band 21-cm measurements of a strong absorption signal around 80~MHz is the presence of an excess radio background to the Cosmic Microwave Background (CMB). Such excess can be produced by the decay of unstable particles into small mass dark photons which have a non-zero mixing angle with electromagnetism. We use the EDGES-LOW band measurements to derive joint constraints on the properties of the early galaxies and the parameters of such a particle physics model for the excess radio background. A Bayesian analysis shows that a high star formation efficiency and X-ray emission of $4-7 \times 10^{48} ~\rm erg$ per solar mass in stars are required along with a suppression of star formation in halos with virial temperatures $\lesssim 2\times 10^4$ K. The same analysis also suggests a 68 percent credible intervals for the mass of the decaying dark matter particles, it's lifetime, dark photon mass and the mixing angle of the dark and ordinary photon oscillation of $[10^{-3.5}, 10^{-2.4}]$ eV, $[10^{1.1}, 10^{2.7}]\times 13.8 ~\rm Gyr$, $[10^{-12.2}, 10^{-10}]$ eV and $[10^{-7}, 10^{-5.6}]$ respectively. This implies an excess radio background which is $\approx 5.7$ times stronger than the CMB around 80~MHz. This value is a factor $\sim 3$ higher than the previous predictions which used a simplified model for the 21-cm signal.

preprint2022arXiv

Convolutional Neural Network-reconstructed velocity for kinetic SZ detection

We report the detection of the kinetic Sunyaev-Zel'dovich (kSZ) effect in galaxy clusters with a 4.9 sigma significance using the latest 217 GHz Planck map from data release 4. For the detection, we stacked the Planck map at the positions of 30,431 galaxy clusters from the Wen-Han-Liu (WHL) catalog. To align the sign of the kSZ signals, the line-of-sight velocities of galaxy clusters were estimated with a machine-learning approach, in which the relation between the galaxy distribution around a cluster and its line-of-sight velocity was trained through a convolutional neural network. To train our network, we used the simulated galaxies and galaxy clusters in the Magneticum cosmological hydrodynamic simulations. The trained model was applied to the large-scale distribution of the Sloan Digital Sky Survey galaxies to derive the line-of-sight velocities of the WHL galaxy clusters. Assuming a standard beta-model for the intracluster medium, we obtained the gas mass fraction in R500 to be fgas,500 = 0.09 +- 0.02 within the galaxy clusters with the average mass of M500 ~ 1.0 x 10^14 Msun/h.

preprint2021arXiv

Direct detection of the kinetic Sunyaev-Zel'dovich effect in galaxy clusters

We report the direct detection of the kinetic Sunyaev-Zel'dovich (kSZ) effect in galaxy clusters with a 3.5 sigma significance level. The measurement was performed by stacking the Planck map at 217 GHz at the positions of galaxy clusters from the Wen-Han-Liu (WHL) catalog. To avoid the cancelation of positive and negative kSZ signals, we used the large-scale distribution of the Sloan Digital Sky Survey (SDSS) galaxies to estimate the peculiar velocities of the galaxy clusters along the line of sight and incorporated the sign in the velocity-weighted stacking of the kSZ signals. Using this technique, we were able to measure the kSZ signal around galaxy clusters beyond 3R500. Assuming a standard beta-model, we also found that the gas fraction within R500 is fgas,500 = 0.12 +- 0.04 for the clusters with the mass of M500 ~ 1e14 Msun/h. We compared this result to predictions from the Magneticum cosmological hydrodynamic simulations as well as other kSZ and X-ray measurements, most of which show a lower gas fraction than the universal baryon fraction for the same mass of clusters. Our value is statistically consistent with results from the measurements and simulations and also with the universal value within our measurement uncertainty.

preprint2020arXiv

Predictions for the 21cm-galaxy cross-power spectrum observable with SKA and future galaxy surveys

In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE cross-correlations are in fact a powerful probe of the epoch of reionization as they are expected to provide precious information on the progress of reionization and the typical scale of ionized regions at different redshifts. The next generation observations with SKA will have a noise level much lower than those with its precursor radio facilities, introducing a significant improvement in the measurement of the cross-correlations. We find that an SKA-HSC/PFS observation will allow to investigate scales below ~10 Mpc/h and ~60 Mpc/h at z=7.3 and 6.6, respectively. WFIRST will allow to access also higher redshifts, as it is expected to observe spectroscopically ~900 LAEs per square degree and unit redshift in the range 7.5<z<8.5. Because of the reduction of the shot noise compared to HSC and PFS, observations with WFIRST will result in more precise cross-correlations and increased observable scales.

preprint2019arXiv

Observing the redshifted 21 cm signal around a bright QSO at $z\sim 10$

We use hydrodynamics and radiative transfer simulations to study the 21~cm signal around a bright QSO at $z \sim 10$. Due to its powerful UV and X-ray radiation, the QSO quickly increases the extent of the fully ionized bubble produced by the pre-existing stellar type sources, in addition to partially ionize and heat the surrounding gas. As expected, a longer QSO lifetime, $t_{\rm QSO}$, results in a 21~cm signal in emission located at increasingly larger angular radii, $θ$, and covering a wider range of $θ$. Similar features can be obtained with a higher galactic emissivity efficiency, $f_{\rm UV}$, so that determining the origin of a large ionized bubble (i.e. QSO vs stars) is not straightforward. Such degeneracy could be reduced by taking advantage of the finite light traveltime effect, which is expected to affect an HII region produced by a QSO differently from one created by stellar type sources. From an observational point of view, we find that the 21 cm signal around a QSO at various $t_{\rm QSO}$ could be detected by SKA1-low with a high signal-noise ratio (S/N). As a reference, for $t_{\rm QSO} = 10\,\rm Myr$, a S/N $\sim 8$ is expected assuming that no pre-heating of the IGM has taken place due to high-$z$ energetic sources, while it can reach value above 10 in case of pre-heating. Observations of the 21~cm signal from the environment of a high-$z$ bright QSO could then be used to set constraints on its lifetime, as well as to reduce the degeneracy between $f_{\rm UV}$ and $t_{\rm QSO}$.

preprint2019arXiv

The information content of Cosmic Infrared Background anisotropies

We use analytic computations to predict the power spectrum as well as the bispectrum of Cosmic Infrared Background (CIB) anisotropies. Our approach is based on the halo model and takes into account the mean luminosity-mass relation. The model is used to forecast the possibility to simultaneously constrain cosmological, CIB and halo occupation distribution (HOD) parameters in the presence of foregrounds. For the analysis we use wavelengths in eight frequency channels between 200 and 900$\;\mathrm{GHz}$ with survey specifications given by Planck and LiteBird. We explore the sensitivity to the model parameters up to multipoles of $\ell =1000$ using auto- and cross-correlations between the different frequency bands. With this setting, cosmological, HOD and CIB parameters can be constrained to a few percent. Galactic dust is modeled by a power law and the shot noise contribution as a frequency dependent amplitude which are marginalized over. We find that dust residuals in the CIB maps only marginally influence constraints on standard cosmological parameters. Furthermore, the bispectrum yields tighter constraints (by a factor four in $1σ$ errors) on almost all model parameters while the degeneracy directions are very similar to the ones of the power spectrum. The increase in sensitivity is most pronounced for the sum of the neutrino masses. Due to the similarity of degeneracies a combination of both analysis is not needed for most parameters. This, however, might be due to the simplified bias description generally adopted in such halo model approaches.

preprint2017arXiv

Tomographic Intensity Mapping versus Galaxy Surveys: Observing the Universe in H-alpha emission with new generation instruments

The H-alpha line emission is an important probe for a number of fundamental quantities in galaxies, including their number density, star formation rate (SFR) and overall gas content. A new generation of low-resolution intensity mapping probes, e.g. SPHEREx and CDIM, will observe galaxies in H-alpha emission over a large fraction of the sky from the local Universe till a redshift of z ~ 6 to 10, respectively. This will also be the target line for observations by the high-resolution Euclid and WFIRST instruments in the z ~ 0.7 - 2 redshift range. In this paper, we estimate the intensity and power spectra of the H-alpha line in the z ~ 0 - 5 redshift range using observed line luminosity functions (LFs), when possible, and simulations, otherwise. We estimate the significance of our predictions by accounting for the modelling uncertainties (e.g. SFR, extinction, etc.) and observational contamination. We find that Intensity Mapping (IM) surveys can make a statistical detection of the full H-alpha emission between z ~ 0.8 - 5. Moreover, we find that the high-frequency resolution and the sensitivity of the planned CDIM surveys allow for the separation of H-alpha emission from several interloping lines. We explore ways to use the combination of these line intensities to probe galaxy properties. As expected, our study indicates that galaxy surveys will only detect bright galaxies that contribute up to a few percent of the overall H-alpha intensity. However, these surveys will provide important constraints on the high end of the H-alpha LF and put strong constraints on the AGN LF.

preprint2016arXiv

Mapping the low surface brightness Universe in the UV band with Lya emission from IGM filaments

A large fraction of the baryonic matter in the Universe is located in filaments in the intergalactic medium. However, the low surface brightness of these filaments has not yet allowed their direct detection except in very special regions in the circum-galactic medium (CGM). Here we simulate the intensity and spatial fluctuations in Lyman Alpha ${\rm (Lyα)}$ emission from filaments in the intergalactic medium (IGM) and discuss the prospects for the next generation of space based instruments to detect the low surface brightness universe at UV wavelengths. Starting with a high resolution N-body simulation we obtain the dark matter density fluctuations and associate baryons with the dark matter particles assuming that they follow the same spatial distribution. The IGM thermal and ionization state is set by a model of the UV background and by the relevant cooling processes for a hydrogen and helium gas. The ${\rm Lyα}$ emissivity is then estimated, taking into account recombination and collisional excitation processes. We find that the detection of these filaments through their ${\rm Lyα}$ emission is well in the reach of the next generation of UV space based instruments and so it should be achieved in the next decade. The density field is populated with halos and galaxies and their ${\rm Lyα}$ emission is estimated. Galaxies are treated as foregrounds and so we discuss methods to reduce their contamination from observational maps. Finally, we estimate the UV continuum background as a function of the redshift of the ${\rm Lyα}$ emission line and discuss how this continuum can affect observations.

preprint2016arXiv

Systematic biases in low frequency radio interferometric data due to calibration: the LOFAR EoR case

The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregrounds, which could potentially mean a suppression of the 21 cm signal. The second effect is an excess of noise beyond the thermal noise. The excess noise shows fluctuations on small frequency scales, and hence it can not be easily removed by foreground removal or avoidance methods. Our analysis suggests that sidelobes of residual sources due to the chromatic point spread function and ionospheric scintillation can not be the dominant causes of the excess noise. Rather, both the suppression of diffuse foregrounds and the excess noise can occur due to calibration with an incomplete sky model containing predominantly bright discrete sources. We show that calibrating only on bright sources can cause suppression of other signals and introduce an excess noise in the data. The levels of the suppression and excess noise depend on the relative flux of sources which are not included in the model with respect to the flux of modeled sources. We discuss possible solutions such as using only long baselines to calibrate the interferometric gain solutions as well as simultaneous multi-frequency calibration along with their benefits and shortcomings.

preprint2016arXiv

The Effect of Foreground Mitigation Strategy on EoR Window Recovery

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 by instead removing these foregrounds. 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 foreground avoidance and a simplified, frequency independent foreground removal model. We find that while, for current generation experiments, foreground avoidance enables a better recovery at $k_{perp} > 0.6 \mathrm{Mpc}^{-1}$, foreground removal is able to recover significantly more signal at small $k_{los}$ for both current and future experiments. 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 both methods perform well for foreground models with line-of-sight and spatial variations around $0.1\%$ however at levels larger than this foregrounds removal shows a greater signal recovery.

preprint2015arXiv

Cosmology from the EoR/Cosmic Dawn with the SKA

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

preprint2015arXiv

Effects of the sources of reionization on 21-cm redshift-space distortions

The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.

preprint2015arXiv

On the definition of superclusters

To obtain a physically well-motivated definition of superclusters, we proposed in our previous work to select superclusters with an overdensity criterion that selects only those objects that will collapse in the future, including those that are at a turn-around in the present epoch. In this paper we present numerical values for these criteria for a range of standard cosmological models. We express these criteria in terms of a density ratio or, alternatively, as an infall velocity and show that these two criteria give almost identical results. To better illustrate the implications of this definition, we applied our criteria to some prominent structures in the local Universe, the Local supercluster, Shapley supercluster, and the recently reported Laniakea supercluster to understand their future evolution. We find that for the Local and Shapley superclusters, only the central regions will collapse in the future, while Laniakea does not constitute a significant overdensity and will disperse in the future. Finally, we suggest that those superclusters that will survive the accelerating cosmic expansion and collapse in the future be called "superstes-clusters", where "superstes" means survivor in Latin, to distinguish them from traditional superclusters.

preprint2015arXiv

Predictions for the 21cm-galaxy cross-power spectrum observable with LOFAR and Subaru

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.

preprint2015arXiv

Realistic Simulations of the Galactic Polarized Foreground: Consequences for 21-cm Reionization Detection Experiments

Experiments designed to measure the redshifted 21~cm line from the Epoch of Reionization (EoR) are challenged by strong astrophysical foreground contamination, ionospheric distortions, complex instrumental response and other different types of noise (e.g. radio frequency interference). The astrophysical foregrounds are dominated by diffuse synchrotron emission from our Galaxy. Here we present a simulation of the Galactic emission used as a foreground module for the LOFAR- EoR key science project end-to-end simulations. The simulation produces total and polarized intensity over $10^\circ \times 10^\circ$ maps of the Galactic synchrotron and free-free emission, including all observed characteristics of the emission: spatial fluctuations of amplitude and spectral index of the synchrotron emission, together with Faraday rotation effects. The importance of these simulations arise from the fact that the Galactic polarized emission could behave in a manner similar to the EoR signal along the frequency direction. As a consequence, an improper instrumental calibration will give rise to leakages of the polarized to the total signal and mask the desired EoR signal. In this paper we address this for the first time through realistic simulations.

preprint2014arXiv

Constraining the epoch of reionization with the variance statistic: simulations of the LOFAR case

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($Δz$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $Δz = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $Δz = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.

preprint2014arXiv

Probing large scale filaments with HI and $^3$HeII

We explore the observability of the neutral hydrogen (HI) and the singly-ionized isotope helium-3 ($^3$HeII) in the intergalactic medium (IGM) from the Epoch of Reionization down to the local Universe. The hyperfine transition of $^3$HeII, which is not as well known as the HI transition, has energy splitting corresponding to 8 cm. It also has a larger spontaneous decay rate than that of neutral hydrogen, whereas its primordial abundance is much smaller. Although both species are mostly ionized in the IGM, the balance between ionization and recombination in moderately high density regions renders them abundant enough to be observed. We estimate the emission signal of both hyperfine transitions from large scale filamentary structures and discuss the prospects for observing them with current and future radio telescopes. We conclude that HI in filaments is possibly observable even with current telescopes after 100 hours of observation. On the other hand, $^3$HeII is only detectable with future telescopes, such as SKA, after the same amount of time.

preprint2014arXiv

Stars and Reionization: The Cross-Correlation of the 21cm Line and the Near Infrared Background

With improving telescopes, it may now be possible to observe the Epoch of Reionization in multiple ways. We examine two of these observables - the excess light in the near-infrared background that may be due to high redshift stars and ionized HII bubbles, and the 21 cm emission from neutral hydrogen. Because these two forms of emission should result from different, mutually exclusive regions, an anticorrelation should exist between them. We discuss the strengths of using cross-correlations between these observations to learn more about high redshift star formation and reionization history. In particular, we create simulated maps of emission from both the near-infrared background and 21 cm emission. We find that these observations are anticorrelated, with the strongest anticorrelation originating from times when the universe is half ionized. This result is robust and does not depend on the properties of the stars themselves. Rather, it depends on the ionization history. Cross-correlations can provide redshift information, which the near-infrared background cannot provide alone. In addition, cross-correlations can help separate foreground emission from the true high redshift component, making it possible to say with greater certainty that we are indeed witnessing the Epoch of Reionization.

preprint2013arXiv

Clustered Calibration: An Improvement to Radio Interferometric Direction Dependent Self-Calibration

The new generation of radio synthesis arrays, such as LOFAR and SKA, have been designed to surpass existing arrays in terms of sensitivity, angular resolution and frequency coverage. This evolution has led to the development of advanced calibration techniques that ensure the delivery of accurate results at the lowest possible computational cost. However, the performance of such calibration techniques is still limited by the compact, bright sources in the sky, used as calibrators. It is important to have a bright enough source that is well distinguished from the background noise level in order to achieve satisfactory results in calibration. We present "clustered calibration" as a modification to traditional radio interferometric calibration, in order to accommodate faint sources that are almost below the background noise level into the calibration process. The main idea is to employ the information of the bright sources' measured signals as an aid to calibrate fainter sources that are nearby the bright sources. In the case where we do not have bright enough sources, a source cluster could act as a bright source that can be distinguished from background noise. We construct a number of source clusters assuming that the signals of the sources belonging to a single cluster are corrupted by almost the same errors, and each cluster is calibrated as a single source, using the combined coherencies of its sources simultaneously. This upgrades the power of an individual faint source by the effective power of its cluster. We give performance analysis of clustered calibration to show the superiority of this approach compared to the traditional unclustered calibration. We also provide analytical criteria to choose the optimum number of clusters for a given observation in an efficient manner.

preprint2013arXiv

Probing reionization with LOFAR using 21-cm redshift space distortions

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

preprint2013arXiv

Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

preprint2013arXiv

The End of an Era - The Population III to Population II Transition and the Near Infrared Background

There are only a few ways to constrain the Era of Reionization and the properties of high redshift (z>6) stars through observations. Here, we discuss one of these observables - the spectrum of the Near Infrared Background - and how it is potentially affected by the transition from Population III to Population II stars. The stronger Lyman-alpha emission expected from massive Population III stars could result in a 'bump' in the spectrum of the Near Infrared Background (referred to in this work as the Lyman-alpha bump). The strength and shape of this bump can reveal properties of Population III stars. The Lyman-alpha bump is predicted to be higher if Population III stars are more massive and present at lower redshifts. The shape of the bump is governed by the star formation rate and the time it takes Population III stars to transition to Population II stars. If Population III stars are indeed massive, a bump is predicted as long as Population III stars exist at z < 15, even if their star formation rate is as low as 10^-7 M_sun yr^-1 Mpc^-3. This means that there may be some observational signature in the Near Infrared Background of small pockets of metal-free gas forming Population III stars at z ~ 6, even if they are quite rare.

preprint2013arXiv

The Scale of the Problem : Recovering Images of Reionization with GMCA

The accurate and precise removal of 21-cm foregrounds from Epoch of Reionization redshifted 21-cm emission data is essential if we are to gain insight into an unexplored cosmological era. We apply a non-parametric technique, Generalized Morphological Component Analysis or GMCA, to simulated LOFAR-EoR data and show that it has the ability to clean the foregrounds with high accuracy. We recover the 21-cm 1D, 2D and 3D power spectra with high accuracy across an impressive range of frequencies and scales. We show that GMCA preserves the 21-cm phase information, especially when the smallest spatial scale data is discarded. While it has been shown that LOFAR-EoR image recovery is theoretically possible using image smoothing, we add that wavelet decomposition is an efficient way of recovering 21-cm signal maps to the same or greater order of accuracy with more flexibility. By comparing the GMCA output residual maps (equal to the noise, 21-cm signal and any foreground fitting errors) with the 21-cm maps at one frequency and discarding the smaller wavelet scale information, we find a correlation coefficient of 0.689, compared to 0.588 for the equivalently smoothed image. Considering only the central 50% of the maps, these coefficients improve to 0.905 and 0.605 respectively and we conclude that wavelet decomposition is a significantly more powerful method to denoise reconstructed 21-cm maps than smoothing.

preprint2012arXiv

Foreground Removal using FastICA: A Showcase of LOFAR-EoR

We introduce a new implementation of the FastICA algorithm on simulated LOFAR EoR data with the aim of accurately removing the foregrounds and extracting the 21-cm reionization signal. We find that the method successfully removes the foregrounds with an average fitting error of 0.5 per cent and that the 2D and 3D power spectra are recovered across the frequency range. We find that for scales above several PSF scales the 21-cm variance is successfully recovered though there is evidence of noise leakage into the reconstructed foreground components. We find that this blind independent component analysis technique provides encouraging results without the danger of prior foreground assumptions.

preprint2012arXiv

The Epoch of Reionization

The Universe's dark ages end with the formation of the first generation of galaxies. These objects start emitting ultraviolet radiation that carves out ionized regions around them. After a sufficient number of ionizing sources have formed, the ionized fraction of the gas in the Universe rapidly increases until hydrogen becomes fully ionized. This period, during which the cosmic gas went from neutral to ionized, is known as the Universe's Epoch of Reionization . The Epoch of Reionization is related to many fundamental questions in cosmology, such as properties of the first galaxies, physics of (mini-)quasars, formation of very metal-poor stars and a slew of other important research topics in astrophysics. Hence uncovering it will have far reaching implications on the study of structure formation in the early Universe. This chapter reviews the current observational evidence for the occurrence of this epoch, its key theoretical aspects and main characteristics, and finally the various observational probes that promise to uncover it. A special emphasis is put on the redshifted 21 cm probe, the various experiments that are currently being either built or designed, and what we can learn from them about the Epoch of Reionization.

preprint2011arXiv

Clustered Radio Interferometric Calibration

This paper introduces an amendment to radio interferometric calibration of sources below the noise level. The main idea is to employ the information of the stronger sources' measured signals as a plug-in criterion to solve for the weaker ones. For this purpose, we construct a number of source clusters, with centroids mainly near the strongest sources, assuming that the signals of the sources belonging to a single cluster are corrupted by almost the same errors. Due to this characteristic of clusters, each cluster is calibrated as a single source, using all the coherencies of its sources simultaneously. The obtained solutions for every cluster are assigned to all the cluster's sources. An illustrative example reveals the superiority of this calibration compared to the un-clustered calibration.

preprint2011arXiv

Modification of the halo mass function by kurtosis associated with primordial non-Gaussianity

We study the halo mass function in the presence of the kurtosis type of primordial non-Gaussianity. The kurtosis corresponds to the trispectrum as defined in Fourier space. The primordial trispectrum is commonly characterized by two parameters, $τ_{\rm NL}$ and $g_{\rm NL}$. As applications of the derived non-Gaussian mass function, we consider the effect on the abundance of void structure, the effect on early star formation and on formation of the most massive object at high redshift. We show that by comparing the effects of primordial non-Gaussianity on cluster abundance with that on void abundance, we can distinguish between the skewness and the kurtosis types of primordial non-Gaussianity. As for early star formation, we show that the kurtosis type of primordial non-Gaussianity seems not to affect the reionization history of the Universe on average. However, at high redshifts (up to $z\simeq 20$) such non-Gaussianity does somewhat affect the early stages of reionization.

preprint2011arXiv

Second order cross-correlation between kSZ and 21 cm fluctuations from the EoR

The measurement of the brightness temperature fluctuations of neutral hydrogen 21 cm lines from the Epoch of Reionisation (EoR) is expected to be a powerful tool for revealing the reionisation process. We study the 21 cm cross-correlation with Cosmic Microwave Background (CMB) temperature anisotropies, focusing on the effect of the patchy reionisation. We calculate, up to second order, the angular power spectrum of the cross-correlation between 21 cm fluctuations and the CMB kinetic Sunyaev-Zel'dovich effect (kSZ) from the EoR, using an analytical reionisation model. We show that the kSZ and the 21 cm fluctuations are anti-correlated on the scale corresponding to the typical size of an ionised bubble at the observed redshift of the 21 cm fluctuations. The amplitude of the angular power spectrum of the cross-correlation depends on the fluctuations of the ionised fraction. Especially, in a highly inhomogeneous reionisation model, the amplitude reaches the order of $100 μK^2$ at $\ell \sim 3000$. We also show that second order terms may help in distinguishing between reionisation histories.

preprint2010arXiv

Detectability of the 21 cm-CMB cross-correlation from the EoR

The 21-cm line fluctuations and the cosmic microwave background (CMB) are powerful probes of the epoch of reionisation of the universe. We study the potential of the cross-correlation between 21-cm line fluctuations and CMB anisotropy to obtain further constraints on the reionisation history. We compute analytically the 21-cm cross-correlation with the CMB temperature anisotropy and polarisation, and we calculate the signal-to-noise (SN) ratio for its detection with Planck together with LOFAR, MWA and SKA. We find, on the one hand, that the 21-cm cross-correlation signal with CMB polarisation from the instant reionisation can be detected with an SN ratio of $\sim 1$ for LOFAR and $\sim 10$ for SKA. On the other hand, we confirm that the detection of the 21-cm cross-correlation with CMB polarisation is practically infeasible.

preprint2010arXiv

On the spin-temperature evolution during the epoch of reionization

Simulations estimating the differential brightness temperature of the redshifted 21-cm from the epoch of reionization (EoR) often assume that the spin temperature is decoupled from the background CMB temperature and is much larger than it. Although a valid assumption towards the latter stages of the reionization process, it does not necessarily hold at the earlier epochs. Violation of this assumption will lead to fluctuations in differential brightness temperature that are neither driven by density fluctuations nor by HII regions. Therefore, it is vital to calculate the spin temperature self-consistently by treating the Lyman-alpha and collisional coupling of spin temperature to the kinetic temperature. In this paper we develop an extension to the BEARS algorithm, originally developed to model reionization history, to include these coupling effects. Here we simulate the effect in ionization and heating for three models in which the reionization is driven by stars, miniqsos or a mixture of both.We also perform a number of statistical tests to quantify the imprint of the self-consistent inclusion of the spin temperature decoupling from the CMB. We find that the evolution of the spin temperature has an impact on the measured signal specially at redshifts higher than 10 and such evolution should be taken into account when one attempts to interpret the observational data.

preprint2010arXiv

Power spectrum extraction for redshifted 21-cm epoch of reionization experiments: the LOFAR case

One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in particular from the subtraction of foreground contamination. We model the extraction process using realistic simulations of the cosmological signal, the foregrounds and noise, and so estimate the sensitivity of the LOFAR EoR experiment to the redshifted 21-cm power spectrum. Detection of emission from the EoR should be possible within 360 hours of observation with a single station beam. Integrating for longer, and synthesizing multiple station beams within the primary (tile) beam, then enables us to extract progressively more accurate estimates of the power at a greater range of scales and redshifts. We discuss different observational strategies which compromise between depth of observation, sky coverage and frequency coverage. A plan in which lower frequencies receive a larger fraction of the time appears to be promising. We also study the nature of the bias which foreground fitting errors induce on the inferred power spectrum, and discuss how to reduce and correct for this bias. The angular and line-of-sight power spectra have different merits in this respect, and we suggest considering them separately in the analysis of LOFAR data.

preprint2010arXiv

Probing the Epoch of Reionization with Low Frequency Arrays

The Epoch of Reionization (EoR) is the epoch in which hydrogen in the Universe reionize after the "Dark Ages". This is the second of two major phase transitions that hydrogen in the Universe underwent, the first phase being the recombination era in which hydrogen became neutral at redshift about 1100. The EoR, occurs around z of 10 and is probably caused by the first radiation emitting astrophysical sources, hence it is crucial to our understanding of when and how the Universe "decided" to start forming astrophysical objects and how that influenced subsequent structure formation in the Universe. As such, the EoR is related to many fundamental questions in cosmology, galaxy formation, quasars and very metal poor stars; all are foremost research issues in modern astrophysics. The redshifted 21 cm hyperfine line is widely considered as the most promising probe for studying the EoR in detail. In the near future a number of low frequency radio telescopes (LOFAR, MWA, GMRT and SKA) will be able to observe the 21 cm radiation arriving from the high redshift Universe. In this paper I present our current picture of the ionization process, review the 21 cm line physics and discuss the challenges that the current generation experiments are expected to face. Finally, I discuss the potential of SKA in exploring the EoR and the Universe's Dark Ages.

preprint2009arXiv

Cross-correlation Study between the Cosmological 21-cm Signal and the kinetic Sunyaev-Zel'dovich effect

The Universe's Epoch of Reionization (EoR) 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 (kSZ) -- produced by the scattering of CMB photons off free electrons produced during the reionization process -- and the cosmological 21cm 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 100Mpc/h scale N-body simulations with radiative transfer. In essence we find that the two probes anti-correlate. The significance of the anti-correlation 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.

preprint2003arXiv

Redshift estimation of clusters by wavelet decomposition of their Sunyaev-Zel'dovich morphology

A method for estimating redshifts of galaxy clusters based solely on resolved Sunyaev-Zel'dovich (SZ) images is proposed. Given a high resolution SZ cluster image (with FWHM of approx. 1 arcmin), the method indirectly measures its structure related parameters (amplitude, size, etc.) by fitting a model function to the higher order wavelet momenents of the cluster's SZ morphology. The applicability and accuracy of the wavelet method is assessed by applying it to maps of a set of clusters extracted from hydrodynamical simulations of cosmic structure formation. The parameters, derived by a fit to the spectrum of wavelet moments as a function of scale, are found to show a dependence on redshift z that is of the type x(z) = x_1 exp(-z/x_2) + x_3, where the monotony of this functional behaviour and the non-degeneracy of those parameters allow inversion and estimation of the redshift z. The average attainable accuracy in the z-estimation relative to 1+z is approx. 4-5% out to z = 1.2, which is comparable to photometric redshifts. For single-frequency SZ interferometers, where the ambient fluctuating CMB is the main noise source, the accuracy of the method drops slightly to <Delta z/(1+z)> = 6-7%.

preprint1994arXiv

Wiener Reconstruction of Density, Velocity, and Potential Fields from All-Sky Galaxy Redshift Surveys

We present a new method for recovering the cosmological density, velocity, and potential fields from all-sky redshift catalogues. The method is based on an expansion of the fields in orthogonal radial (Bessel) and angular (spherical harmonic) functions. In this coordinate system, peculiar velocities introduce a coupling of the radial harmonics describing the density field in redshift space but leave the angular modes unaffected. In the harmonic transform space, this radial coupling is described by a distortion matrix which can be computed analytically within the context of linear theory; the redshift space harmonics can then be converted to their real space values by inversion of this matrix. Statistical noise is mitigated by regularizing the matrix inversion with a Wiener filter. The method yields a minimum variance estimate of the density field in real space. In this coordinate system, the minimum variance harmonics of the peculiar velocity and potential fields are related to those of the density field by simple linear transformations. Tests of the method with simulations of a CDM universe and comparison with previously proposed methods demonstrate it to be a very promising new reconstruction method for the local density and velocity field. A first application to the 1.2 Jy IRAS redshift survey is presented.