Source author record

Aaron R. Parsons

Aaron R. Parsons 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

30works
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

30 published item(s)

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.

preprint2021arXiv

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.

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

Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

preprint2020arXiv

Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1σ$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $κ_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $κ_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

preprint2020arXiv

Foreground modelling via Gaussian process regression: an application to HERA data

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

preprint2020arXiv

Measuring HERA's primary beam in-situ: methodology and first results

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

preprint2019arXiv

Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

preprint2016arXiv

Constraining cosmology and ionization history with combined 21 cm power spectrum and global signal measurements

Improvements in current instruments and the advent of next-generation instruments will soon push observational 21 cm cosmology into a new era, with high significance measurements of both the power spectrum and the mean ("global") signal of the 21 cm brightness temperature. In this paper we use the recently commenced Hydrogen Epoch of Reionization Array as a worked example to provide forecasts on astrophysical and cosmological parameter constraints. In doing so we improve upon previous forecasts in a number of ways. First, we provide updated forecasts using the latest best-fit cosmological parameters from the Planck satellite, exploring the impact of different Planck datasets on 21 cm experiments. We also show that despite the exquisite constraints that other probes have placed on cosmological parameters, the remaining uncertainties are still large enough to have a non-negligible impact on upcoming 21 cm data analyses. While this complicates high-precision constraints on reionization models, it provides an avenue for 21 cm reionization measurements to constrain cosmology. We additionally forecast HERA's ability to measure the ionization history using a combination of power spectrum measurements and semi-analytic simulations. Finally, we consider ways in which 21 cm global signal and power spectrum measurements can be combined, and propose a method by which power spectrum results can be used to train a compact parameterization of the global signal. This parameterization reduces the number of parameters needed to describe the global signal, increasing the likelihood of a high significance measurement.

preprint2016arXiv

Eliminating the optical depth nuisance from the CMB with 21 cm cosmology

Amongst standard model parameters that are constrained by cosmic microwave background (CMB) observations, the optical depth $τ$ stands out as a nuisance parameter. While $τ$ provides some crude limits on reionization, it also degrades constraints on other cosmological parameters. Here we explore how 21 cm cosmology---as a direct probe of reionization---can be used to independently predict $τ$ in an effort to improve CMB parameter constraints. We develop two complementary schemes for doing so. The first uses 21 cm power spectrum observations in conjunction with semi-analytic simulations to predict $τ$. The other uses global 21 cm measurements to directly constrain low redshift (post-reheating) contributions to $τ$ in a relatively model-independent way. Forecasting the performance of the upcoming Hydrogen Epoch of Reionization Array, we find that significant reductions in the errors on $τ$ can be achieved. These results are particularly effective at breaking the CMB degeneracy between $τ$ and the amplitude of the primordial fluctuation spectrum $A_s$, with errors on $\ln (10^{10} A_s)$ reduced by up to a factor of four. Stage 4 CMB constraints on the neutrino mass sum are also improved, with errors potentially reduced to $12\,\textrm{meV}$ regardless of whether CMB experiments can precisely measure the reionization bump in polarization power spectra. Observations of the 21 cm line are therefore capable of improving not only our understanding of reionization astrophysics, but also of cosmology in general.

preprint2016arXiv

Optimized Beam Sculpting with Generalized Fringe-Rate Filters

We generalize the technique of fringe-rate filtering, whereby visibilities measured by a radio interferometer are re-weighted according to their temporal variation. As the Earth rotates, radio sources traverse through an interferometer's fringe pattern at rates that depend on their position on the sky. Capitalizing on this geometric interpretation of fringe rates, we employ time-domain convolution kernels to enact fringe-rate filters that sculpt the effective primary beam of antennas in an interferometer. As we show, beam sculpting through fringe-rate filtering can be used to optimize measurements for a variety of applications, including mapmaking, minimizing polarization leakage, suppressing instrumental systematics, and enhancing the sensitivity of power-spectrum measurements. We show that fringe-rate filtering arises naturally in minimum variance treatments of many of these problems, enabling optimal visibility-based approaches to analyses of interferometric data that avoid systematics potentially introduced by traditional approaches such as imaging. Our techniques have recently been demonstrated in Ali et al. (2015), where new upper limits were placed on the 21 cm power spectrum from reionization, showcasing the ability of fringe-rate filtering to successfully boost sensitivity and reduce the impact of systematics in deep observations.

preprint2016arXiv

Radio Emission from Red-Giant Hot Jupiters

When planet-hosting stars evolve off the main sequence and go through the red-giant branch, the stars become orders of magnitudes more luminous and, at the same time, lose mass at much higher rates than their main-sequence counterparts. Accordingly, if planetary companions exist around these stars at orbital distances of a few AU, they will be heated up to the level of canonical hot Jupiters and also be subjected to a dense stellar wind. Given that magnetized planets interacting with stellar winds emit radio waves, such "Red-Giant Hot Jupiters" (RGHJs) may also be candidate radio emitters. We estimate the spectral auroral radio intensity of RGHJs based on the empirical relation with the stellar wind as well as a proposed scaling for planetary magnetic fields. RGHJs might be intrinsically as bright as or brighter than canonical hot Jupiters and about 100 times brighter than equivalent objects around main-sequence stars. We examine the capabilities of low-frequency radio observatories to detect this emission and find that the signal from an RGHJ may be detectable at distances up to a few hundred parsecs with the Square Kilometer Array.

preprint2016arXiv

Redundant Array Configurations for 21 cm Cosmology

Realizing the potential of 21 cm tomography to statistically probe the intergalactic medium before and during the Epoch of Reionization requires large telescopes and precise control of systematics. Next-generation telescopes are now being designed and built to meet these challenges, drawing lessons from first-generation experiments that showed the benefits of densely packed, highly redundant arrays--in which the same mode on the sky is sampled by many antenna pairs--for achieving high sensitivity, precise calibration, and robust foreground mitigation. In this work, we focus on the Hydrogen Epoch of Reionization Array (HERA) as an interferometer with a dense, redundant core designed following these lessons to be optimized for 21 cm cosmology. We show how modestly supplementing or modifying a compact design like HERA's can still deliver high sensitivity while enhancing strategies for calibration and foreground mitigation. In particular, we compare the imaging capability of several array configurations, both instantaneously (to address instrumental and ionospheric effects) and with rotation synthesis (for foreground removal). We also examine the effects that configuration has on calibratability using instantaneous redundancy. We find that improved imaging with sub-aperture sampling via "off-grid" antennas and increased angular resolution via far-flung "outrigger" antennas is possible with a redundantly calibratable array configuration.

preprint2016arXiv

Spherical Harmonic Analyses of Intensity Mapping Power Spectra

Intensity mapping is a promising technique for surveying the large scale structure of our Universe from $z=0$ to $z \sim 150$, using the brightness temperature field of spectral lines to directly observe previously unexplored portions of out cosmic timeline. Examples of targeted lines include the $21\,\textrm{cm}$ hyperfine transition of neutral hydrogen, rotational lines of carbon monoxide, and fine structure lines of singly ionized carbon. Recent efforts have focused on detections of the power spectrum of spatial fluctuations, but have been hindered by systematics such as foreground contamination. This has motivated the decomposition of data into Fourier modes perpendicular and parallel to the line-of-sight, which has been shown to be a particularly powerful way to diagnose systematics. However, such a method is well-defined only in the limit of a narrow-field, flat-sky approximation. This limits the sensitivity of intensity mapping experiments, as it means that wide surveys must be separately analyzed as a patchwork of smaller fields. In this paper, we develop a framework for analyzing intensity mapping data in a spherical Fourier-Bessel basis, which incorporates curved sky effects without difficulty. We use our framework to generalize a number of techniques in intensity mapping data analysis from the flat sky to the curved sky. These include visibility-based estimators for the power spectrum, treatments of interloper lines, and the "foreground wedge" signature of spectrally smooth foregrounds.

preprint2015arXiv

Mapmaking for Precision 21 cm Cosmology

In order to study the "Cosmic Dawn" and the Epoch of Reionization with 21 cm tomography, we need to statistically separate the cosmological signal from foregrounds known to be orders of magnitude brighter. Over the last few years, we have learned much about the role our telescopes play in creating a putatively foreground-free region called the "EoR window." In this work, we examine how an interferometer's effects can be taken into account in a way that allows for the rigorous estimation of 21 cm power spectra from interferometric maps while mitigating foreground contamination and thus increasing sensitivity. This requires a precise understanding of the statistical relationship between the maps we make and the underlying true sky. While some of these calculations would be computationally infeasible if performed exactly, we explore several well-controlled approximations that make mapmaking and the calculation of map statistics much faster, especially for compact and highly redundant interferometers designed specifically for 21 cm cosmology. We demonstrate the utility of these methods and the parametrized trade-offs between accuracy and speed using one such telescope, the upcoming Hydrogen Epoch of Reionization Array, as a case study.

preprint2015arXiv

PAPER-64 Constraints On Reionization II: The Temperature Of The z=8.4 Intergalactic Medium

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

preprint2014arXiv

Multi-redshift limits on the 21cm power spectrum from PAPER

The epoch of reionization power spectrum is expected to evolve strongly with redshift, and it is this variation with cosmic history that will allow us to begin to place constraints on the physics of reionization. The primary obstacle to the measurement of the EoR power spectrum is bright foreground emission. We present an analysis of observations from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) telescope which place new limits on the HI power spectrum over the redshift range of $7.5<z<10.5$, extending previously published single redshift results to cover the full range accessible to the instrument. To suppress foregrounds, we use filtering techniques that take advantage of the large instrumental bandwidth to isolate and suppress foreground leakage into the interesting regions of $k$-space. Our 500 hour integration is the longest such yet recorded and demonstrates this method to a dynamic range of $10^4$. Power spectra at different points across the redshift range reveal the variable efficacy of the foreground isolation. Noise limited measurements of $Δ^2$ at $k=$0.2hMpc$^{-1}$ and z$=7.55$ reach as low as (48mK)$^2$ ($1σ$). We demonstrate that the size of the error bars in our power spectrum measurement as generated by a bootstrap method is consistent with the fluctuations due to thermal noise. Relative to this thermal noise, most spectra exhibit an excess of power at a few sigma. The likely sources of this excess include residual foreground leakage, particularly at the highest redshift, and unflagged RFI. We conclude by discussing data reduction improvements that promise to remove much of this excess.

preprint2014arXiv

New Limits on 21cm EoR From PAPER-32 Consistent with an X-Ray Heated IGM at z=7.7

We present new constraints on the 21cm Epoch of Reionization (EoR) power spectrum derived from 3 months of observing with a 32-antenna, dual-polarization deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. In this paper, we demonstrate the efficacy of the delay-spectrum approach to avoiding foregrounds, achieving over 8 orders of magnitude of foreground suppression (in $\textrm{mK}^2$). Combining this approach with a procedure for removing off-diagonal covariances arising from instrumental systematics, we achieve a best 2-sigma upper limit of $(41\,\textrm{mK})^2$ for $k=0.27 h\textrm{Mpc}^{-1}$ at $z=7.7$. This limit falls within an order of magnitude of the brighter predictions of the expected 21cm EoR signal level. Using the upper limits set by these measurements, we generate new constraints on the brightness temperature of 21cm emission in neutral regions for various reionization models. We show that for several ionization scenarios, our measurements are inconsistent with cold reionization. That is, heating of the neutral intergalactic medium (IGM) is necessary to remain consistent with the constraints we report. Hence, we have suggestive evidence that by $z=7.7$, the HI has been warmed from its cold primordial state, probably by X-rays from high-mass X-ray binaries or mini-quasars. The strength of this evidence depends on the ionization state of the IGM, which we are not yet able to constrain. This result is consistent with standard predictions for how reionization might have proceeded.

preprint2014arXiv

The Epoch of Reionization Window: I. Mathematical Formalism

The 21 cm line provides a powerful probe of astrophysics and cosmology at high redshifts, but unlocking the potential of this probe requires the robust mitigation of foreground contaminants that are typically several orders of magnitude brighter than the cosmological signal. Recent simulations and observations have shown that the smooth spectral structure of foregrounds combines with instrument chromaticity to contaminate a "wedge"-shaped region in cylindrical Fourier space. While previous efforts have explored the suppression of foregrounds within this wedge, as well as the avoidance of this highly contaminated region, all such efforts have neglected a rigorous examination of the error statistics associated with the wedge. Using a quadratic estimator formalism applied to the interferometric measurement equation, we provide a framework for such a rigorous analysis (incorporating a fully covariant treatment of errors). Additionally, we find that there are strong error correlations at high spatial wavenumbers that have so far been neglected in sensitivity derivations. These error correlations substantially degrade the sensitivity of arrays relying on contributions from long baselines, compared to what one would estimate assuming uncorrelated errors.

preprint2014arXiv

The Epoch of Reionization Window: II. Statistical Methods for Foreground Wedge Reduction

For there to be a successful measurement of the 21 cm Epoch of Reionization (EoR) power spectrum, it is crucial that strong foreground contaminants be robustly suppressed. These foregrounds come from a variety of sources (such as Galactic synchrotron emission and extragalactic point sources), but almost all share the property of being spectrally smooth, and when viewed through the chromatic response of an interferometer, occupy a signature "wedge" region in cylindrical $k_\perp k_\parallel$ Fourier space. The complement of the foreground wedge is termed the "EoR window", and is expected to be mostly foreground-free, allowing clean measurements of the power spectrum. This paper is a sequel to a previous paper that established a rigorous mathematical framework for describing the foreground wedge and the EoR window. Here, we use our framework to explore statistical methods by which the EoR window can be enlarged, thereby increasing the sensitivity of a power spectrum measurement. We adapt the FKP approximation (commonly used in galaxy surveys) for 21 cm cosmology, and also compare the optimal quadratic estimator to simpler estimators that ignore covariances between different Fourier modes. The optimal quadratic estimator is found to suppress foregrounds by an extra factor of $\sim 10^5$ in power at the peripheries of the EoR window, boosting the detection of the cosmological signal from $12σ$ to $50σ$ at the midpoint of reionization in our fiducial models. If numerical issues can be finessed, decorrelation techniques allow the EoR window to be further enlarged, enabling measurements to be made deep within the foreground wedge. These techniques do not assume that foreground are Gaussian-distributed, and we additionally prove that a final round of foreground subtraction can be performed after decorrelation in a way that is guaranteed to have no cosmological signal loss.

preprint2014arXiv

What Next-Generation 21 cm Power Spectrum Measurements Can Teach Us About the Epoch of Reionization

A number of experiments are currently working towards a measurement of the 21 cm signal from the Epoch of Reionization. Whether or not these experiments deliver a detection of cosmological emission, their limited sensitivity will prevent them from providing detailed information about the astrophysics of reionization. In this work, we consider what types of measurements will be enabled by a next-generation of larger 21 cm EoR telescopes. To calculate the type of constraints that will be possible with such arrays, we use simple models for the instrument, foreground emission, and the reionization history. We focus primarily on an instrument modeled after the $\sim 0.1~\rm{km}^2$ collecting area Hydrogen Epoch of Reionization Array (HERA) concept design, and parameterize the uncertainties with regard to foreground emission by considering different limits to the recently described "wedge" footprint in k-space. Uncertainties in the reionization history are accounted for using a series of simulations which vary the ionizing efficiency and minimum virial temperature of the galaxies responsible for reionization, as well as the mean free path of ionizing photons through the IGM. Given various combinations of models, we consider the significance of the possible power spectrum detections, the ability to trace the power spectrum evolution versus redshift, the detectability of salient power spectrum features, and the achievable level of quantitative constraints on astrophysical parameters. Ultimately, we find that $0.1~\rm{km}^2$ of collecting area is enough to ensure a very high significance ($\gtrsim30σ$) detection of the reionization power spectrum in even the most pessimistic scenarios. This sensitivity should allow for meaningful constraints on the reionization history and astrophysical parameters, especially if foreground subtraction techniques can be improved and successfully implemented.

preprint2013arXiv

A Flux Scale for Southern Hemisphere 21cm EoR Experiments

We present a catalog of spectral measurements covering a 100-200 MHz band for 32 sources, derived from observations with a 64-antenna deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. For transit telescopes such as PAPER, calibration of the primary beam is a difficult endeavor, and errors in this calibration are a major source of error in the determination of source spectra. In order to decrease reliance on accurate beam calibration, we focus on calibrating sources in a narrow declination range from -46d to -40d. Since sources at similar declinations follow nearly identical paths through the primary beam, this restriction greatly reduces errors associated with beam calibration, yielding a dramatic improvement in the accuracy of derived source spectra. Extrapolating from higher frequency catalogs, we derive the flux scale using a Monte-Carlo fit across multiple sources that includes uncertainty from both catalog and measurement errors. Fitting spectral models to catalog data and these new PAPER measurements, we derive new flux models for Pictor A and 31 other sources at nearby declinations. 90% of these confirm and refine a power-law model for flux density. Of note is the new Pictor A flux model, which is accurate to 1.4% and shows, in contrast to previous models, that between 100 MHz and 2 GHz, the spectrum of Pictor A is consistent with a single power law given by a flux at 150 MHz of 382+/-5.4 Jy, and a spectral index of -0.76+/-0.01. This accuracy represents an order of magnitude improvement over previous measurements in this band, and is limited by the uncertainty in the catalog measurements used to estimate the absolute flux scale. The simplicity and improved accuracy of Pictor A's spectrum make it an excellent calibrator for experiments seeking to measure 21cm emission from the Epoch of Reionization.

preprint2013arXiv

Opening the 21cm EoR Window: Measurements of Foreground Isolation with PAPER

We present new observations with the Precision Array for Probing the Epoch of Reionization (PAPER) with the aim of measuring the properties of foreground emission for 21cm Epoch of Reionization experiments at 150 MHz. We focus on the footprint of the foregrounds in cosmological Fourier space to understand which modes of the 21cm power spectrum will most likely be compromised by foreground emission. These observations confirm predictions that foregrounds can be isolated to a "wedge"-like region of 2D (k-perpendicular, k-parallel)-space, creating a window for cosmological studies at higher k-parallel values. We also find that the emission extends past the nominal edge of this wedge due to spectral structure in the foregrounds, with this feature most prominent on the shortest baselines. Finally, we filter the data to retain only this "unsmooth" emission and image specific k-parallel modes of it. The resultant images show an excess of power at the lowest modes, but no emission can be clearly localized to any one region of the sky. This image is highly suggestive that the most problematic foregrounds for 21cm EoR studies will not be easily identifiable bright sources, but rather an aggregate of fainter emission.

preprint2013arXiv

The Baryon Acoustic Oscillation Broadband and Broad-beam Array: Design Overview and Sensitivity Forecasts

This work describes a new instrument optimized for a detection of the neutral hydrogen 21cm power spectrum between redshifts of 0.5-1.5: the Baryon Acoustic Oscillation Broadband and Broad-beam (BAOBAB) Array. BAOBAB will build on the efforts of a first generation of 21cm experiments which are targeting a detection of the signal from the Epoch of Reionization at z ~ 10. At z ~ 1, the emission from neutral hydrogen in self-shielded overdense halos also presents an accessible signal, since the dominant, synchrotron foreground emission is considerably fainter than at redshift 10. The principle science driver for these observations are Baryon Acoustic Oscillations in the matter power spectrum which have the potential to act as a standard ruler and constrain the nature of dark energy. BAOBAB will fully correlate dual-polarization antenna tiles over the 600-900MHz band with a frequency resolution of 300 kHz and a system temperature of 50K. The number of antennas will grow in staged deployments, and reconfigurations of the array will allow for both traditional imaging and high power spectrum sensitivity operations. We present calculations of the power spectrum sensitivity for various array sizes, with a 35-element array measuring the cosmic neutral hydrogen fraction as a function of redshift, and a 132-element system detecting the BAO features in the power spectrum, yielding a 1.8% error on the z ~ 1 distance scale, and, in turn, significant improvements to constraints on the dark energy equation of state over an unprecedented range of redshifts from ~0.5-1.5.

preprint2013arXiv

The Effects of Polarized Foregrounds on 21cm Epoch of Reionization Power Spectrum Measurements

Experiments aimed at detecting highly-redshifted 21 centimeter emission from the Epoch of Reionization (EoR) are plagued by the contamination of foreground emission. A potentially important source of contaminating foregrounds may be Faraday-rotated, polarized emission, which leaks into the estimate of the intrinsically unpolarized EoR signal. While these foregrounds' intrinsic polarization may not be problematic, the spectral structure introduced by the Faraday rotation could be. To better understand and characterize these effects, we present a simulation of the polarized sky between 120 and 180 MHz. We compute a single visibility, and estimate the three-dimensional power spectrum from that visibility using the delay spectrum approach presented in Parsons et al. (2012b) . Using the Donald C. Backer Precision Array to Probe the Epoch of Reionization (PAPER) as an example instrument, we show the expected leakage into the unpolarized power spectrum to be several orders of magnitude above the expected 21cm EoR signal.

preprint2012arXiv

A Per-Baseline, Delay-Spectrum Technique for Accessing the 21cm Cosmic Reionization Signature

A critical challenge in measuring the power spectrum of 21cm emission from cosmic reionization is compensating for the frequency dependence of an interferometer's sampling pattern, which can cause smooth-spectrum foregrounds to appear unsmooth and degrade the separation between foregrounds and the target signal. In this paper, we present an approach to foreground removal that explicitly accounts for this frequency dependence. We apply the delay transformation introduced in Parsons & Backer (2009) to each baseline of an interferometer to concentrate smooth-spectrum foregrounds within the bounds of the maximum geometric delays physically realizable on that baseline. By focusing on delay-modes that correspond to image-domain regions beyond the horizon, we show that it is possible to avoid the bulk of smooth-spectrum foregrounds. We map the point-spread function of delay-modes to k-space, showing that delay-modes that are uncorrupted by foregrounds also represent samples of the three-dimensional power spectrum, and can be used to constrain cosmic reionization. Because it uses only spectral smoothness to differentiate foregrounds from the targeted 21cm signature, this per-baseline analysis approach relies on spectrally- and spatially-smooth instrumental responses for foreground removal. For sufficient levels of instrumental smoothness relative to interfering foregrounds, this technique substantially reduces the level of calibration previously thought necessary to detect 21cm reionization. As a result, this approach places fewer constraints on antenna configuration within an array, facilitating the adoption of configurations optimized for power-spectrum sensitivity. Under these assumptions, we demonstrate the potential for PAPER to detect 21cm reionization at an amplitude of 10 mK^2 near k~0.2h Mpc^-1 with 132 dipoles in 7 months of observing.

preprint2011arXiv

A Technique for Primary Beam Calibration of Drift-Scanning, Wide-Field Antenna Elements

We present a new technique for calibrating the primary beam of a wide-field, drift-scanning antenna element. Drift-scan observing is not compatible with standard beam calibration routines, and the situation is further complicated by difficult-to-parametrize beam shapes and, at low frequencies, the sparsity of accurate source spectra to use as calibrators. We overcome these challenges by building up an interrelated network of source "crossing points" -- locations where the primary beam is sampled by multiple sources. Using the single assumption that a beam has 180 degree rotational symmetry, we can achieve significant beam coverage with only a few tens of sources. The resulting network of crossing points allows us to solve for both a beam model and source flux densities referenced to a single calibrator source, circumventing the need for a large sample of well-characterized calibrators. We illustrate the method with actual and simulated observations from the Precision Array for Probing the Epoch of Reionization (PAPER).

preprint2011arXiv

New 145-MHz Source Measurements by PAPER in the Southern Sky

We present observations from the Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa, observed in May and September 2010. Using two nights of drift scanning with PAPER's 60\arcdeg\ FWHM beam we have made a map covering the entire sky below +10 degrees declination with an effective center frequency of 145 MHz, a 70-MHz bandwidth, and a resolution of 26\arcmin. A 4800 square-degree region of this large map with the lowest Galactic emission reaches an RMS of 0.7 Jy. We establish an absolute flux scale using sources from the 160-MHz Culgoora catalog. Using the 408-MHz Molonglo Reference Catalog (MRC) as a finding survey, we identify counterparts to 480 sources in our maps, and compare our fluxes to the MRC and to 332 sources in the Culgoora catalog. For both catalogs, the ratio of PAPER to catalog flux averages to 1, with a standard deviation of 50%. This measured variation is consistent with comparisons between independent catalogs observed at different bands. The PAPER data represent new 145-MHz flux measurements for a large number of sources in the band expected to encompass cosmic reionization, and represents a significant step toward establishing a model for removing foregrounds to the reionization signal.

preprint2009arXiv

Calibration of Low-Frequency, Wide-Field Radio Interferometers Using Delay/Delay-Rate Filtering

We present a filtering technique that can be applied to individual baselines of wide-bandwidth, wide-field interferometric data to geometrically select regions on the celestial sphere that contain primary calibration sources. The technique relies on the Fourier transformation of wide-band frequency spectra from a given baseline to obtain one-dimensional "delay images", and then the transformation of a time-series of delay images to obtain two-dimensional "delay/delay-rate images." Source selection is possible in these images given appropriate combinations of baseline, bandwidth, integration time and source location. Strong and persistent radio frequency interference (RFI) limits the effectiveness of this source selection owing to the removal of data by RFI excision algorithms. A one-dimensional, complex CLEAN algorithm has been developed to compensate for RFI-excision effects. This approach allows CLEANed, source-isolated data to be used to isolate bandpass and primary beam gain functions. These techniques are applied to data from the Precision Array for Probing the Epoch of Reionization (PAPER) as a demonstration of their value in calibrating a new generation of low-frequency radio interferometers with wide relative bandwidths and large fields-of-view.

preprint2009arXiv

The Precision Array for Probing the Epoch of Reionization: 8 Station Results

We are developing the Precision Array for Probing the Epoch of Reionization (PAPER) to detect 21cm emission from the early Universe, when the first stars and galaxies were forming. We describe the overall experiment strategy and architecture and summarize two PAPER deployments: a 4-antenna array in the low-RFI environment of Western Australia and an 8-antenna array at our prototyping site in Green Bank, WV. From these activities we report on system performance, including primary beam model verification, dependence of system gain on ambient temperature, measurements of receiver and overall system temperatures, and characterization of the RFI environment at each deployment site. We present an all-sky map synthesized between 139 MHz and 174 MHz using data from both arrays that reaches down to 80 mJy (4.9 K, for a beam size of 2.15e-5 steradians at 154 MHz), with a 10 mJy (620 mK) thermal noise level that indicates what would be achievable with better foreground subtraction. We calculate angular power spectra ($C_\ell$) in a cold patch and determine them to be dominated by point sources, but with contributions from galactic synchrotron emission at lower radio frequencies and angular wavemodes. Although the cosmic variance of foregrounds dominates errors in these power spectra, we measure a thermal noise level of 310 mK at $\ell=100$ for a 1.46-MHz band centered at 164.5 MHz. This sensitivity level is approximately three orders of magnitude in temperature above the level of the fluctuations in 21cm emission associated with reionization.