Source author record

Mario G. Santos

Mario G. Santos 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

40works
4topics
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

40 published item(s)

preprint2022arXiv

Antenna characterization for the HIRAX experiment

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) aims to improve constraints on the dark energy equation of state through measurements of large-scale structure at high redshift ($0.8<z<2.5$), while serving as a state-of-the-art fast radio burst detector. Bright galactic foregrounds contaminate the 400--800~MHz HIRAX frequency band, so meeting the science goals will require precise instrument characterization. In this paper we describe characterization of the HIRAX antenna, focusing on measurements of the antenna beam and antenna noise temperature. Beam measurements of the current HIRAX antenna design were performed in an anechoic chamber and compared to simulations. We report measurement techniques and results, which find a broad and symmetric antenna beam for $ν<$650MHz, and elevated cross-polarization levels and beam asymmetries for $ν>$700MHz. Noise temperature measurements of the HIRAX feeds were performed in a custom apparatus built at Yale. In this system, identical loads, one cryogenic and the other at room temperature, are used to take a differential (Y-factor) measurement from which the noise of the system is inferred. Several measurement sets have been conducted using the system, involving CHIME feeds as well as four of the HIRAX active feeds. These measurements give the first noise temperature measurements of the HIRAX feed, revealing a $\sim$60K noise temperature (relative to 30K target) with 40K peak- to-peak frequency-dependent features, and provide the first demonstration of feed repeatability. Both findings inform current and future feed designs.

preprint2022arXiv

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets.

preprint2022arXiv

Measurements of the diffuse Galactic synchrotron spectral index and curvature from MeerKLASS pilot data

21cm intensity mapping experiments are bringing an influx of high spectral resolution observational data in the $\sim100$ MHz $- 1$ GHz regime. We use pilot $971-1075$ MHz data from MeerKAT in single-dish mode, recently used to test the calibration and data reduction scheme of the upcoming MeerKLASS survey, to probe the spectral index of diffuse synchrotron emission below 1 GHz within $145^{\circ} < α< 180^{\circ}$, $-1^{\circ} < δ< 8^{\circ}$. Through comparisons with data from the OVRO Long Wavelength Array and the Maipu and MU surveys, we find an average spectral index of $-2.75 < β< -2.71$ between 45 and 1055 MHz. By fitting for spectral curvature with a spectral index of the form $β+ c \, {\rm{ln}}(ν/ 73~{\rm MHz})$, we measure $β= -2.55 \pm 0.13$ and $c = -0.12 \pm 0.05$ within our target field. Our results are in good agreement (within $1σ$) with existing measurements from experiments such as ARCADE2 and EDGES. These results show the calibration accuracy of current data and demonstrate that MeerKLASS will also be capable of achieving a secondary science goal of probing the interstellar medium.

preprint2022arXiv

MIGHTEE-HI: Evolution of HI scaling relations of star-forming galaxies at $z<0.5$

We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with ${\rm sSFR}\equiv M_*/{\rm SFR}$), obtaining $\gtrsim 5σ$ detections in most cases, the strongest HI-stacking detections to date in this redshift range. With these detections, we are able to measure scaling relations in the probed redshift interval, finding evidence for a moderate evolution from the median redshift of our sample $z_{\rm med}\sim 0.37$ to $z\sim 0$. In particular, low-$M_*$ galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 9$) experience a strong HI depletion ($\sim 0.5$ dex in $\log_{10}(M_{\rm HI}/{\rm M}_\odot)$), while massive galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 11$) keep their HI mass nearly unchanged. When looking at the star formation activity, highly star-forming galaxies evolve significantly in $M_{\rm HI}$ ($f_{\rm HI}$, where $f_{\rm HI}\equiv M_{\rm}/M_*$) at fixed SFR (sSFR), while at the lowest probed SFR (sSFR) the scaling relations show no evolution. These findings suggest a scenario in which low-$M_*$ galaxies have experienced a strong HI depletion during the last $\sim4$ Gyr, while massive galaxies have undergone a significant HI replenishment through some accretion mechanism, possibly minor mergers. Interestingly, our results are in good agreement with the predictions of the SIMBA simulation. We conclude that this work sets novel important observational constraints on galaxy scaling relations.

preprint2022arXiv

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.

preprint2021arXiv

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.

preprint2020arXiv

A Full Treatment of Peculiar Velocities on the Reionization Light Cone

Accurate simulations of the 21-cm signal from the Epoch of Reionization require the generation of maps at specific frequencies given the values of astrophysical and cosmological parameters. The peculiar velocities of the hydrogen atoms producing the 21-cm radiation result in a shift in the observed frequency of the 21-cm radiation and changes the amplitude of the signal itself. This is not an effect we can remove but instead needs to be accurately modelled to ensure we infer the correct physical parameters from an observation. We describe the full calculation of the distortion of the 21-cm signal, and propose a new code that integrates the 21-cm intensity along the line of sight for each individual light cone pixel to fully capture the intensity contributions from multiple redshifts. This algorithm naturally deals with the typical divergences found in standard approaches, allowing for large optical depths and 21-cm absorption events at high redshifts. We find the new method results in up to a 60% decrease in power on the largest scales on the sky, and an increase of over 80% on the smallest scales on the sky. We find that the new implementation of the light cone results in a longer tail of bright temperatures in the brightness temperature distribution, as a result of the successful circumventing of a previous cap that had to be implemented to avoid a divergence in the brightness temperature. We conclude that this full treatment of the evolution of the light cone pixel can be an important effect.

preprint2020arXiv

ARTIST: Fast radiative transfer for large-scale simulations of the epoch of reionisation

We introduce the "Asymmetric Radiative Transfer In Shells Technique" (ARTIST), a new method for photon propagation on large scales that explicitly conserves photons, propagates photons at the speed of light, approximately accounts for photon directionality, and closely reproduces results of more detailed radiative transfer (RT) codes. Crucially, it is computationally fast enough to evolve the large cosmological volumes required to predict the 21cm power spectrum on scales that will be probed by future experiments targeting the Epoch of Reionisation (EoR). Most semi-numerical models aimed at predicting the EoR 21cm signal make use of an excursion set formalism (ESF) approach, which achieves computational viability by compromising on photon conservation, constraining ionised regions to be spherical by construction, and not accounting for light-travel time. By implementing our RT method within the semi-numerical code SimFast21, we show that ARTIST predicts a significantly different evolution for the EoR ionisation field compared to the code's native ESF. In particular, ARTIST predicts a more gradual evolution of the volume-averaged ionisation fraction, and up to an order-of-magnitude difference in the ionisation power, depending on the physical parameters assumed. Its application to large-scale EoR simulations will therefore allow more physically-motivated constraints to be obtained for key EoR parameters, such as the escape fraction.

preprint2020arXiv

Constraining the reionization history using deep learning from 21cm tomography with the Square Kilometre Array

Upcoming 21cm surveys with the SKA1-LOW telescope will enable imaging of the neutral hydrogen distribution on cosmological scales in the early Universe. These surveys are expected to generate huge imaging datasets that will encode more information than the power spectrum. This provides an alternative unique way to constrain the reionization history, which might break the degeneracy in the power spectral analysis. Using Convolutional Neural Networks (CNN), we create a fast estimator of the neutral fraction from the 21cm maps that are produced by our large semi-numerical simulation. Our estimator is able to efficiently recover the neutral fraction ($x_{\rm HI}$) at several redshifts with a high accuracy of 99\% as quantified by the coefficient of determination $R^{2}$. Adding the instrumental effects from the SKA design slightly increases the loss function, but nevertheless we are still able to recover the neutral fraction with a similar high accuracy of 98\%, which is only 1 per cent less. While a weak dependence on redshift is observed, the accuracy increases rapidly with decreasing neutral fraction. This is due to the fact that the instrumental noise increases towards high redshift where the Universe is highly neutral. Our results show the promise of directly using 21cm-tomography to constrain the reionization history in a model independent way, complementing similar efforts, such as those of the optical depth measurements from the Cosmic Microwave Background (CMB) observations by {\it Planck}.

preprint2020arXiv

Prospects for cosmic magnification measurements using HI intensity mapping

We investigate the prospects of measuring the cosmic magnification effect by cross-correlating neutral hydrogen intensity mapping (HI IM) maps with background optical galaxies. We forecast the signal-to-noise ratio for HI IM data from SKA1-MID and HIRAX, combined with LSST photometric galaxy samples. We find that, thanks to their different resolutions, SKA1-MID and HIRAX are highly complementary in such an analysis. We predict that SKA1-MID can achieve a detection with a signal-to-noise ratio of $\sim 15$ on a multipole range of $\ell \lesssim 200$, while HIRAX can reach a signal-to-noise ratio of $\sim 50$ on $200 < \ell < 2000$. We conclude that measurements of the cosmic magnification signal will be possible on a wide redshift range with foreground HI intensity maps up to $z \lesssim 2$, while optimal results are obtained when $0.6 \lesssim z \lesssim 1.3$. Finally, we perform a signal to noise analysis that shows how these measurements can constrain the HI parameters across a wide redshift range.

preprint2020arXiv

The optically-selected 1.4-GHz quasar luminosity function below 1 mJy

We present the radio luminosity function (RLF) of optically-selected quasars below 1~mJy, constructed by applying a Bayesian-fitting stacking technique to objects well below the nominal radio flux-density limit. We test the technique using simulated data, confirming that we can reconstruct the RLF over three orders of magnitude below the typical $5σ$ detection threshold. We apply our method to 1.4-GHz flux-densities from the Faint Images of the Radio Sky at Twenty-cm survey (FIRST), extracted at the positions of optical quasars from the Sloan Digital Sky Survey (SDSS) over seven redshift bins up to $z=2.15$ {and measure the RLF down to two orders of magnitude below the FIRST detection threshold}. In the lowest redshift bin ($0.2<z<0.45$), we find that our measured RLF agrees well with deeper data from the literature. The RLF for the radio-loud quasars flattens below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 25.5$ and becomes steeper again below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 24.8$, where radio-quiet quasars start to emerge. The radio luminosity where radio-quiet quasars emerge coincides with the luminosity where star-forming galaxies are expected to start to dominate the radio source counts. This implies that there could be a significant contribution from star formation in the host galaxies, but additional data is required to investigate this further. The higher-redshift bins show a similar behaviour as for the lowest-$z$ bin, implying that the same physical process may be responsible.

preprint2019arXiv

Measuring the HI mass function below the detection threshold

We present a Bayesian Stacking technique to directly measure the HI mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several sky areas given an assumed HIMF described by a Schechter function and simulate the HI emission lines with different levels of background noise to test the technique. We use Multinest to constrain the parameters of the HIMF in a broad redshift bin, demonstrating that the HIMF can be accurately reconstructed, using the simulated spectral cube far below the HI mass limit determined by the $5σ$ flux-density limit, i.e. down to $M_{\rm HI} = 10^{7.5}$ M$_{\odot}$ over the redshift range $0 < z < 0.55$ for this particular simulation, with a noise level similar to that expected for the MIGHTEE survey. We also find that the constraints on the parameters of the Schechter function, $ϕ_{\star}$, $M_\star$ and $α$ can be reliably fit, becoming tighter as the background noise decreases as expected, although the constraints on the redshift evolution are not significantly affected. All the parameters become better constrained as the survey area increases. In summary, we provide an optimal method for estimating the HI mass at cosmological distances that allows us to constrain the HI mass function below the detection threshold in forthcoming HI surveys. This study is a first step towards the measurement of the HIMF at high ($z>0.1$) redshifts.

preprint2018arXiv

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

preprint2016arXiv

Simulating the 21-cm signal from reionisation including non-linear ionisations and inhomogeneous recombinations

We explore the impact of incorporating physically motivated ionisation and recombination rates on the history and topology of cosmic reionisation, by incorporating inputs from small-volume hydrodynamic simulations into a semi-numerical code, SimFast21, that evolves reionisation on large scales. We employ radiative hydrodynamic simulations to parameterize the ionisation rate Rion and recombination rate Rrec as functions of halo mass, overdensity and redshift. We find that Rion is super-linearly dependent on halo mass (Rion ~ Mh^1.41), in contrast to previous assumptions. We implement these scalings into SimFast21 to identify the ionized regions. We tune our models to be consistent with recent observations of the optical depth, ionizing emissivity, and neutral fraction by the end of reionisation. We require an average photon escape fraction fesc=0.04 within ~ 0.5 cMpc cells, independent of halo mass or redshift, to simultaneously match these data. We present predictions for the 21cm power spectrum, and show that it is converged with respect to simulation volume. We find that introducing superlinearly mass-dependent ionisations increases the duration of reionisation and boosts the small-scale 21cm power by ~ 2-3 at intermediate phases of reionisation. Introducing inhomogeneous recombinations reduces ionised bubble sizes and suppresses large-scale 21cm power by ~ 2-3. Moreover, gas clumping on sub-cell scales has a minimal effect on the 21cm power, indicating that robust predictions do not depend on the behaviour of kpc-scale structures. The superlinear ionisations significantly increase the median halo mass scale for ionising photon output to >10^10 Mo, giving greater hope for detecting most of ionising sources with next-generation facilities. These results highlight the importance of more accurately treating ionising sources and recombinations for modeling reionisation and its 21cm signal.

preprint2015arXiv

Cosmology with a SKA HI intensity mapping survey

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

preprint2015arXiv

Cross-correlating 21cm intensity maps with Lyman Break Galaxies in the post-reionization era

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

preprint2015arXiv

Einstein's legacy in galaxy surveys

Non-Gaussianity in the primordial fluctuations that seeded structure formation produces a signal in the galaxy power spectrum on very large scales. This signal contains vital information about the primordial Universe, but it is very challenging to extract, because of cosmic variance and large-scale systematics - especially after the Planck experiment has already ruled out a large amplitude for the signal. Cosmic variance and experimental systematics can be alleviated by the multi-tracer method. Here we address another systematic - introduced by not using the correct relativistic analysis of the power spectrum on very large scales. In order to reduce the errors on fNL, we need to include measurements on the largest possible scales. Failure to include the relativistic effects on these scales can introduce significant bias in the best-fit value of fNL from future galaxy surveys.

preprint2015arXiv

Epoch of Reionization modelling and simulations for SKA

In this chapter we provide an overview of the current status of the simulations and modelling of the Cosmic Dawn and Epoch of Reionization. We discuss the modelling requirements as dictated by the characteristic scales of the problem and the SKA instrumental properties and the planned survey parameters. Current simulations include most of the relevant physical processes. They can follow the full nonlinear dynamics and are now reaching the required scale and dynamic range, although small-scale physics still needs to be included at sub-grid level. However, despite a significant progress in developing novel numerical methods for efficient utilization of current hardware they remain quite computationally expensive. In response, a number of alternative approaches, particularly semi-analytical/semi-numerical methods, have been developed. While necessarily more approximate, if appropriately constructed and calibrated on simulations they could be used to quickly explore the vast parameter space available. Further work is still required on including some physical processes in both simulations and semi-analytical modelling. This hybrid approach of fast, approximate modelling calibrated on numerical simulations can then be used to construct large libraries of reionization models for reliable interpretation of the observational data.

preprint2015arXiv

HI galaxy simulations for the SKA: number counts and bias

This chapter describes the assumed specifications and sensitivities for HI galaxy surveys with SKA1 and SKA2. It addresses the expected galaxy number densities based on available simulations as well as the clustering bias over the underlying dark matter. It is shown that a SKA1 HI galaxy survey should be able to find around $5\times 10^6$ galaxies over 5,000 deg$^2$ (up to $z\sim 0.8$), while SKA2 should find $\sim 10^9$ galaxies over 30,000 deg$^2$ (up to $z\sim 2.5$). The numbers presented here have been used throughout the cosmology chapters for forecasting.

preprint2015arXiv

Late-time cosmology with 21cm intensity mapping experiments

We present a framework for forecasting cosmological constraints from future neutral hydrogen intensity mapping experiments at low to intermediate redshifts. In the process, we establish a simple way of comparing such surveys with optical galaxy redshift surveys. We explore a wide range of experimental configurations and assess how well a number of cosmological observables (the expansion rate, growth rate, and angular diameter distance) and parameters (the densities of dark energy and dark matter, spatial curvature, the dark energy equation of state, etc.) will be measured by an extensive roster of upcoming experiments. A number of potential contaminants and systematic effects are also studied in detail. The overall picture is encouraging -- if autocorrelation calibration can be controlled to a sufficient level, Phase I of the SKA should be able to constrain the dark energy equation of state about as well as a DETF Stage IV galaxy redshift survey like Euclid, in roughly the same timeframe.

preprint2015arXiv

Measuring baryon acoustic oscillations with future SKA surveys

The imprint of baryon acoustic oscillations (BAO) in large-scale structure can be used as a standard ruler for mapping out the cosmic expansion history, and hence for testing cosmological models. In this article we briefly describe the scientific background to the BAO technique, and forecast the potential of the Phase 1 and 2 SKA telescopes to perform BAO surveys using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys. We find that a 25,000 sq. deg. intensity mapping survey on a Phase 1 array will preferentially constrain the radial BAO, providing a highly competitive 2% constraint on the expansion rate at z ~ 2. A 30,000 sq. deg. galaxy redshift survey on SKA2 will outperform all other planned experiments for z < 1.4.

preprint2015arXiv

Probing primordial non-Gaussianity via iSW measurements with SKA continuum surveys

The Planck CMB experiment has delivered the best constraints so far on primordial non-Gaussianity, ruling out early-Universe models of inflation that generate large non-Gaussianity. Although small improvements in the CMB constraints are expected, the next frontier of precision will come from future large-scale surveys of the galaxy distribution. The advantage of such surveys is that they can measure many more modes than the CMB -- in particular, forthcoming radio surveys with the SKA will cover huge volumes. Radio continuum surveys deliver the largest volumes, but with the disadvantage of no redshift information. In order to mitigate this, we use two additional observables. First, the integrated Sachs-Wolfe effect -- the cross-correlation of the radio number counts with the CMB temperature anisotropies -- helps to reduce systematics on the large scales that are sensitive to non-Gaussianity. Second, optical data allows for cross-identification in order to gain some redshift information. We show that, while the single redshift bin case can provide a sigma(fNL) ~ 20, and is therefore not competitive with current and future constraints on non-Gaussianity, a tomographic analysis could improve the constraints by an order of magnitude, even with only two redshift bins. A huge improvement is provided by the addition of high-redshift sources, so having cross-ID for high-z galaxies and an even higher-z radio tail is key to enabling very precise measurements of fNL. Our results show that SKA continuum surveys could provide constraints competitive with CMB and forthcoming optical surveys, potentially allowing a measurement of sigma(fNL) ~ 1 to be made. Moreover, these measurements would act as a useful check of results obtained with other probes at other redshift ranges with other methods.

preprint2015arXiv

Ultra large-scale cosmology in next-generation experiments with single tracers

Future surveys of large-scale structure will be able to measure perturbations on the scale of the cosmological horizon, and so could potentially probe a number of novel relativistic effects that are negligibly small on sub-horizon scales. These effects leave distinctive signatures in the power spectra of clustering observables and, if measurable, would open a new window on relativistic cosmology. We quantify the size and detectability of the effects for the most relevant future large-scale structure experiments: spectroscopic and photometric galaxy redshift surveys, intensity mapping surveys of neutral hydrogen, and radio continuum surveys. Our forecasts show that next-generation experiments, reaching out to redshifts $z\simeq 4$, will not be able to detect previously-undetected general-relativistic effects by using individual tracers of the density field, although the contribution of weak lensing magnification on large scales should be clearly detectable. We also perform a rigorous joint forecast for the detection of primordial non-Gaussianity through the excess power it produces in the clustering of biased tracers on large scales, finding that uncertainties of $σ(f_{\rm NL})\sim 1-2$ should be achievable. We study the level of degeneracy of these large-scale effects with several tracer-dependent nuisance parameters, quantifying the minimal priors on the latter that are needed for an optimal measurement of the former. Finally, we discuss the systematic effects that must be mitigated to achieve this level of sensitivity, and some alternative approaches that should help to improve the constraints. The computational tools developed to carry out this study, which requires the full-sky computation of the theoretical angular power spectra for $\mathcal{O}(100)$ redshift bins, as well as realistic models of the luminosity function, are publicly available.

preprint2014arXiv

Blind foreground subtraction for intensity mapping experiments

We make use of a large set of fast simulations of an intensity mapping experiment with characteristics similar to those expected of the Square Kilometre Array (SKA) in order to study the viability and limits of blind foreground subtraction techniques. In particular, we consider different approaches: polynomial fitting, principal component analysis (PCA) and independent component analysis (ICA). We review the motivations and algorithms for the three methods, and show that they can all be described, using the same mathematical framework, as different approaches to the blind source separation problem. We study the efficiency of foreground subtraction both in the angular and radial (frequency) directions, as well as the dependence of this efficiency on different instrumental and modelling parameters. For well-behaved foregrounds and instrumental effects we find that foreground subtraction can be successful to a reasonable level on most scales of interest. We also quantify the effect that the cleaning has on the recovered signal and power spectra. Interestingly, we find that the three methods yield quantitatively similar results, with PCA and ICA being almost equivalent.

preprint2014arXiv

Fast simulations for intensity mapping experiments

We present a code to generate mock observations of 21 cm intensity mapping experiments. The emphasis of the code is on reducing the computational cost of running a full-blown simulation, trading computational time for accuracy. The code can be used to generate independent realizations of the cosmological signal and foregrounds, which are necessary, for instance, in order to obtain realistic forecasts for future intensity mapping experiments. The code is able to reproduce the correct angular and radial clustering pattern for the cosmological signal, including redshift-space distortions, lightcone evolution and bias. Furthermore, it is possible to simulate a variety of foregrounds, including the potentially problematic polarized synchrotron emission.

preprint2014arXiv

Foreground contamination in Ly-alpha intensity mapping during the epoch of reionization

The intensity mapping of Ly-alpha emission during the epoch of reionization (EoR) will be contaminated by foreground emission lines from lower redshifts. We calculate the mean intensity and power spectrum of Ly-alpha emission at z~7, and estimate the uncertainties according to the relevant astrophysical processes. We find that the low-redshift emission lines from 6563 A H-alpha, 5007 A [OIII] and 3727 A [OII] will be the strong contaminants on the observed Ly-alpha power spectrum. We make use of both the star formation rate (SFR) and luminosity functions (LF) to estimate the mean intensity and power spectra of the three foreground lines at z~0.5 for H-alpha, z~0.9 for [OIII] and z~1.6 for [OII], as they will contaminate the Ly-alpha emission at z~7. The [OII] line is found to be the strongest. We analyze the masking of the bright survey pixels with a foreground line above some line intensity threshold as a way to reduce the contamination in the intensity mapping survey. We find that the foreground contamination can be neglected if we remove the pixels with fluxes above 1.4x10^-20 W/m^2.

preprint2014arXiv

Radio Galaxy populations and the multi-tracer technique: pushing the limits on primordial non-Gaussianity

We explore the use of different radio galaxy populations as tracers of different mass halos and therefore, with different bias properties, to constrain primordial non-Gaussianity of the local type. We perform a Fisher matrix analysis based on the predicted auto and cross angular power spectra of these populations, using simulated redshift distributions as a function of detection flux and the evolution of the bias for the different galaxy types (Star forming galaxies, Starburst galaxies, Radio-Quiet Quasars, FRI and FRII AGN galaxies). We show that such a multi-tracer analysis greatly improves the information on non-Gaussianity by drastically reducing the cosmic variance contribution to the overall error budget. By using this method applied to future surveys, we predict a constraint of sigma_fnl=3.6 on the local non-Gaussian parameter for a galaxy detection flux limit of 10 \muJy and sigma_fnl=2.2 for 1 \muJy. We show that this significantly improves on the constraints obtained when using the whole undifferentiated populations (sigma_fnl=48 for 10 \muJy and sigma_fnl=12 for 1 \muJy). We conclude that continuum radio surveys alone have the potential to constrain primordial non-Gaussianity to an accuracy at least a factor of two better than the present constraints obtained with Planck data on the CMB bispectrum, opening a window to obtain sigma_fnl~1 with the Square Kilometer Array.

preprint2013arXiv

Beyond Stacking: A Maximum-Likelihood Method to Constrain Radio Source Counts Below the Detection Threshold

We present a statistical method based on a maximum likelihood approach to constrain the number counts of extragalactic sources below the nominal flux-density limit of continuum imaging surveys. We extract flux densities from a radio map using positional information from an auxiliary catalogue and show that we can model the number counts of this undetected population down to flux density levels well below the detection threshold of the radio survey. We demonstrate the capabilities that our method will have with future generation wide-area radio surveys by performing simulations over various sky areas. We show that it is possible to accurately constrain the number counts of the simulated distribution down to one tenth of the flux noise rms with just a sky area of 100 deg$^2$. We then test the application of our method using data from the Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST). We extract flux densities from the FIRST map, sensitive to 150 $μ$Jy/beam (1 $σ$), using the positional information from a catalogue in the same field, also acquired at the same frequency, sensitive to 12 $μ$Jy/beam (1 $σ$). Implementing our method, with known source positions, we are able to recover the right differential number counts of the noise-dominated FIRST map fluxes down to a flux density level which is one-tenth the FIRST detection threshold.

preprint2013arXiv

Cosmology on Ultralarge Scales with Intensity Mapping of the Neutral Hydrogen 21 cm Emission: Limits on Primordial Non-Gaussianity

The large-scale structure of the Universe supplies crucial information about the physical processes at play at early times. Unresolved maps of the intensity of 21 cm emission from neutral hydrogen HI at redshifts z~1-5 are the best hope of accessing the ultralarge-scale information, directly related to the early Universe. A purpose-built HI intensity experiment may be used to detect the large scale effects of primordial non-Gaussianity, placing stringent bounds on different models of inflation. We argue that it may be possible to place tight constraints on the non-Gaussianity parameter f_NL, with an error close to ~1.

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.

preprint2012arXiv

Intensity Mapping of Lyman-alpha Emission During the Epoch of Reionization

We calculate the absolute intensity and anisotropies of the Lyman-alpha radiation field present during the epoch of reionization. We consider emission from both galaxies and the intergalactic medium (IGM) and take into account the main contributions to the production of Lyman-alpha photons: recombinations, collisions, continuum emission from galaxies and scattering of Lyman-n photons in the IGM. We find that the emission from individual galaxies dominates over the IGM with a total Lyman-alpha intensity (times frequency) of about (1.43-3.57)x10^{-8} erg s^{-1} cm^{-2} sr^{-1} at a redshift of 7. This intensity level is low so it is unlikely that the Lyman-αbackground during reionization can be established by an experiment aiming at an absolute background light measurement. Instead we consider Lyman-αintensity mapping with the aim of measuring the anisotropy power spectrum which has rms fluctuations at the level of 1 x 10^{-16} [erg s^[-1} cm^{-2} sr^{-1}]^2 at a few Mpc scales. These anisotropies could be measured with a spectrometer at near-IR wavelengths from 0.9 to 1.4 μm with fields in the order of 0.5 to 1 sq. degrees. We recommend that existing ground-based programs using narrow band filters also pursue intensity fluctuations to study statistics on the spatial distribution of faint Lyman-αemitters. We also discuss the cross-correlation signal with 21 cm experiments that probe HI in the IGM during reionization. A dedicated sub-orbital or space-based Lyman-αintensity mapping experiment could provide a viable complimentary approach to probe reionization, when compared to 21 cm experiments, and is likely within experimental reach.

preprint2012arXiv

The Near-IR Background Intensity and Anisotropies During The Epoch of Reionization

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The absolute intensity and the anisotropies of the near-IR background provide an observational probe of the first-light galaxies and their spatial distribution. We estimate the extragalactic background light intensity during reionization by accounting for the stellar and nebular emission from first-light galaxies. We require the UV photon density from these galaxies to generate a reionization history that is consistent with the optical depth to electron scattering from cosmic microwave background measurements. We also require the bright-end luminosity function of galaxies in our models to reproduce the measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 0.1 to 0.3 nW m$^{-2}$ sr$^{-1}$ at the peak of its spectrum at $\sim$ 1.1 $μ$m. We also discuss the anisotropy power spectrum of the near-IR background using a halo model to describe the galaxy distribution. We compare our predictions for the anisotropy power spectrum to existing measurements from deep near-IR imaging data from {\it Spitzer}/IRAC, {\it Hubble}/NICMOS, and {\it AKARI}. The predicted rms fluctuations at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing measurements. While strong arguments have been made that the measured fluctuations do not have an origin involving faint low-redshift galaxies, we find that the existing measurements are also incompatible with an origin during the era of reionization. The measured near-IR background anisotropies remain unexplained and could be associated with an unidentified non-astrophysical origin.

preprint2011arXiv

Cosmological Measurements with Forthcoming Radio Continuum Surveys

We present forecasts for constraints on cosmological models which can be obtained by forthcoming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the planned surveys to predict measurements obtained with the source auto-correlation, the cross-correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the cross-correlation of radio sources with foreground objects due to cosmic magnification, and a joint analysis together with the CMB power spectrum and supernovae. We show that near future radio surveys will bring complementary measurements to other experiments, probing different cosmological volumes, and having different systematics. Our results show that the unprecedented sky coverage of these surveys combined should provide the most significant measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the ISW effect will significantly tighten constraints on modified gravity parameters, while the best measurements of dark energy models will come from galaxy auto-correlation function analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will be possible to measure the dark energy parameters with an uncertainty of \{$σ(w_0) = 0.05$, $σ(w_a) = 0.12$\} and the modified gravity parameters \{$σ(η_0) = 0.10$, $σ(μ_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors. Finally, we show that radio surveys would detect a primordial non-Gaussianity of $f_{\rm NL}$ = 8 at 1-$σ$ and we briefly discuss other promising probes.

preprint2011arXiv

Intensity Mapping of the [CII] Fine Structure Line during the Epoch of Reionization

The atomic CII fine-structure line is one of the brightest lines in a typical star-forming galaxy spectrum with a luminosity ~ 0.1% to 1% of the bolometric luminosity. It is potentially a reliable tracer of the dense gas distribution at high redshifts and could provide an additional probe to the era of reionization. By taking into account of the spontaneous, stimulated and collisional emission of the CII line, we calculate the spin temperature and the mean intensity as a function of the redshift. When averaged over a cosmologically large volume, we find that the CII emission from ionized carbon in individual galaxies is larger than the signal generated by carbon in the intergalactic medium (IGM). Assuming that the CII luminosity is proportional to the carbon mass in dark matter halos, we also compute the power spectrum of the CII line intensity at various redshifts. In order to avoid the contamination from CO rotational lines at low redshift when targeting a CII survey at high redshifts, we propose the cross-correlation of CII and 21-cm line emission from high redshifts. To explore the detectability of the CII signal from reionization, we also evaluate the expected errors on the CII power spectrum and CII-21 cm cross power spectrum based on the design of the future milimeter surveys. We note that the CII-21 cm cross power spectrum contains interesting features that captures physics during reionization, including the ionized bubble sizes and the mean ionization fraction, which are challenging to measure from 21-cm data alone. We propose an instrumental concept for the reionization CII experiment targeting the frequency range of $\sim$ 200 to 300 GHz with 1, 3 and 10 meter apertures and a bolometric spectrometer array with 64 independent spectral pixels with about 20,000 bolometers.

preprint2011arXiv

Primordial non-Gaussianity from the 21 cm Power Spectrum during the Epoch of Reionization

Primordial non-Gaussianity is a crucial test of inflationary cosmology. We consider the impact of non-Gaussianity on the ionization power spectrum from 21 cm emission during the epoch of reionization. We focus on the power spectrum on large scales at redshifts of 7 to 8 and explore the expected constraint on the local non-Gaussianity parameter f_NL for current and next-generation 21 cm experiments. We show that experiments such as SKA and MWA could measure f_NL values of order 10. This can be improved by an order of magnitude with a fast-Fourier transform telescope like Omniscope.

preprint2011arXiv

Probing Reionization with Intensity Mapping of Molecular and Fine Structure Lines

We propose observations of the molecular gas distribution during the era of reionization. At $z \sim $ 6 to 8, $^{12}$CO$(J=1-0)$ line intensity results in a mean brightness temperature of about 0.5 $μ$K with rms fluctuations of 0.1 $μ$K at 1 to 10 Mpc spatial scales, corresponding to 30 arcminute angular scales. This intensity fluctuations can be mapped with an interferometer, similar to existing and planned 21-cm background experiments, but operating at $\sim$ 12 to 17 GHz. We discuss the feasibility of detecting the cross-correlation between HI and CO molecular gas since such a cross-correlation has the advantage that it will be independent of systematics and most foregrounds in each of the 21-cm and CO$(1-0)$ line experiments. Additional instruments tuned to higher-order transitions of the CO molecule or an instrument operating with high spectral resolution at mm-wavelengths targeting 158 $μ$m CII could further improve the reionization studies with molecular gas. The combined 21-cm and CO line observations has the potential to establish the relative distribution of gas in the inter-galactic medium and molecular gas that are clumped in individual first-light galaxies that are closely connected to the formation of massive stars in these galaxies.

preprint2011arXiv

The OH line contamination of 21 cm intensity fluctuation measurements for z=1~4

The large-scale structure of the Universe can be mapped with unresolved intensity fluctuations of the 21 cm line. The power spectrum of the intensity fluctuations has been proposed as a probe of the baryon acoustic oscillations at low to moderate redshifts with interferometric experiments now under consideration. We discuss the contamination to the low-redshift 21 cm intensity power spectrum generated by the 18 cm OH line since the intensity fluctuations of the OH line generated at a slightly higher redshift contribute to the intensity fluctuations observed in an experiment. We assume the OH megamaser luminosity is correlated with the star formation rate, and use the simulation to estimate the OH signal and the spatial anisotropies. We also use a semi-analytic simulation to predict the 21 cm power spectrum. At z=1 to 3, we find that the OH contamination could reach 0.1% to 1% of the 21 cm rms fluctuations at the scale of the first peak of the baryon acoustic oscillation. When z>3 the OH signal declines quickly, so that the contamination on the 21-cm becomes negligible at high redshifts.

preprint2010arXiv

A Constraint On the Integrated Mass Power Spectrum out to z = 1100 from Lensing of the Cosmic Microwave Background

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.