Source author record

Judd D. Bowman

Judd D. Bowman 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

38works
6topics
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

38 published item(s)

preprint2022arXiv

Analytic approximations of scattering effects on beam chromaticity in 21-cm global experiments

Scattering from objects near an antenna produce correlated signals from strong compact radio sources in a manner similar to those used by the Sea Interferometer to measure the radio source positions using the fine frequency structure in the total power spectrum of a single antenna. These fringes or ripples due to correlated signal interference are present at a low level in the spectrum of any single antenna and are a major source of systematics in systems used to measure the global redshifted 21-cm signal from the early universe. In the Sea Interferometer a single antenna on a cliff above the sea is used to add the signal from the direct path to the signal from the path reflected from the sea thereby forming an interferometer. This was used for mapping radio sources with a single antenna by Bolton and Slee in the 1950s. In this paper we derive analytic expressions to determine the level of these ripples and compare these results in a few simple cases with electromagnetic modeling software to verify that the analytic calculations are sufficient to obtain the magnitude of the scattering effects on the measurements of the global 21-cm signal. These analytic calculations are needed to evaluate the magnitude of the effects in cases that are either too complex or take too much time to be modeled using software.

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 of Diffuse Radio Surveys at 45 and 150 MHz

We use EDGES measurements to determine scale and zero-level corrections to the diffuse radio surveys by Guzmán et al. at $45$ MHz and Landecker & Wielebinski at $150$ MHz. We find that the Guzmán et al. map requires a scale correction of $1.076 \pm 0.034$ ($2σ$) and a zero-level correction of $-160 \pm 78$ K ($2σ$) to best-fit the EDGES data. For the Landecker & Wielebinski map, the scale correction is $1.112 \pm 0.023$ ($2σ$) and the zero-level correction is $0.7 \pm 6.0$ K ($2σ$). The correction uncertainties are dominated by systematic effects, of which the most significant are uncertainty in the calibration of the EDGES receivers, antenna pointing, and tropospheric and ionospheric effects. We propagate the correction uncertainties to estimate the uncertainties in the corrected maps themselves and find that the $2σ$ uncertainty in the map brightness temperature is in the range $3.2-7.5\%$ for the Guzmán et al. map and $2.1-9.0\%$ for the Landecker & Wielebinski map, with the largest percent uncertainties occurring at high Galactic latitudes. The corrected maps could be used to improve existing diffuse low-frequency radio sky models, which are essential tools in analyses of cosmological $21$ cm observations, as well as to investigate the existence of a radio monopole excess above the cosmic microwave background and known Galactic and extragalactic contributions.

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

CHIPS: The Cosmological HI Power Spectrum Estimator

Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$σ$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $Δ_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

preprint2016arXiv

Delay Spectrum with Phase-Tracking Arrays: Extracting the HI power spectrum from the Epoch of Reionization

The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data.

preprint2016arXiv

First Limits on the 21 cm Power Spectrum during the Epoch of X-ray heating

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window.

preprint2016arXiv

First Season MWA EoR Power Spectrum Results at Redshift 7

The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of $Δ^2 \leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications.

preprint2016arXiv

GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue

Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307,455 radio sources with 20 separate flux density measurements across 72--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5σ$, which itself is typically 50mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low frequency flux density scale of the southern sky to better than 10%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results.

preprint2016arXiv

Impact of Instrument Responses on the Detectability of One-point Statistics from Redshifted 21 cm Observations

We study the impact of instrumental systematics on the variance, skewness, and kurtosis of redshifted 21 cm intensity fluctuation observations from the Epoch of Reionization. We simulate realistic 21 cm observations based on the Murchison Widefield Array (MWA) Phase I reionization experiment, using the array's point spread function (PSF) and antenna beam patterns, full-sky 21 cm models, and the FHD imaging pipeline. We measure the observed redshift evolution of pixel probability density functions (PDF) and one-point statistics from the simulated maps, comparing them to the measurements derived from simpler simulations that represent the instrument PSFs with Gaussian kernels. We find that both methods yield statistics with similar trends with greater than 80% correlation. We perform additional simulations based on the Hydrogen Epoch of Reionization Array (HERA), using Gaussian kernels as the instrument PSFs, and study the effect of frequency binning on the statistics. We find that PSF smoothing and sampling variance from measuring the statistics over limited field of view dilute intrinsic features and add fluctuations to the statistics but reveal new detectable features. Observed kurtosis will increase when a few extremely high or low temperature regions are present in the maps. Frequency binning reduces the thermal uncertainty but can also blur regions along the frequency dimension, resulting in kurtosis peaks that only appear in statistics derived from maps of certain frequency bins. We further find that the kurtosis peaks will reach their maxima when the angular resolution of the PSFs match the size scale of the extreme regions that produce the peaks. The HERA array should be capable of charting the evolution of the observed skewness and kurtosis of the 21 cm fluctuations with high sensitivity while the MWA Phase I will likely be capable of detecting the peak in variance.

preprint2016arXiv

One-Port Direct/Reverse Method for Characterizing VNA Calibration Standards

This paper introduces a one-port method for estimating model parameters of VNA calibration standards. The method involves measuring the standards through an asymmetrical passive network connected in direct mode and then in reverse mode, and using these measurements to compute the S-parameters of the network. The free parameters of the calibration standards are estimated by minimizing a figure of merit based on the expected equality of the S-parameters of the network when used in direct and reverse modes. The capabilities of the method are demonstrated through simulations, and real measurements are used to estimate the actual offset delay of a 50-$\mathbfΩ$ calibration load that is assigned zero delay by the manufacturer. The estimated delay is $38.8$ ps with a $1σ$ uncertainty of $2.1$ ps for this particular load. This result is verified through measurements of a terminated airline. The measurements agree better with theoretical models of the airline when the reference plane is calibrated using the new estimate for the load delay.

preprint2016arXiv

The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology

We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds.

preprint2015arXiv

Confirmation of Wide-Field Signatures in Redshifted 21 cm Power Spectra

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

preprint2015arXiv

Foregrounds in Wide-Field Redshifted 21 cm Power Spectra

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.

preprint2015arXiv

Limits on Foreground Subtraction from Chromatic Beam Effects in Global Redshifted 21 cm Measurements

Foreground subtraction in global redshifted 21 cm measurements is limited by frequency-dependent (chromatic) structure in antenna beam patterns. Chromatic beams couple angular structures in Galactic foreground emission to spectral structures that may not be removed by smooth functional forms. We report results for simulations based on two dipole antennas used by the Experiment to Detect the Global EoR Signature (EDGES). The residual levels in simulated foreground-subtracted spectra are found to differ substantially between the two antennas, suggesting that antenna design must be carefully considered. Residuals are also highly dependent on the right ascension and declination of the antenna pointing, with RMS values differing by as much as a factor of 20 across pointings. For EDGES and other ground-based experiments with zenith pointing antennas, right ascension and declination correspond directly to the local sidereal time and the latitude of the deployment site, hence chromatic beam effects should be taken into account when selecting sites. We introduce the "blade" dipole antenna and show, via simulations, that it has better chromatic performance than the "fourpoint" antenna previously used for EDGES. The blade antenna yields 1-5~mK residuals across the entire sky after a 5-term polynomial is removed from simulated spectra, whereas the fourpoint antenna typically requires a 6-term polynomial for comparable residuals. For both antennas, the signal-to-noise ratio of recovered 21 cm input signals peaks for a 5-term polynomial foreground fit given realistic thermal noise levels.

preprint2014arXiv

Overcoming real-world obstacles in 21 cm power spectrum estimation: A method demonstration and results from early Murchison Widefield Array data

We present techniques for bridging the gap between idealized inverse covariance weighted quadratic estimation of 21 cm power spectra and the real-world challenges presented universally by interferometric observation. By carefully evaluating various estimators and adapting our techniques for large but incomplete data sets, we develop a robust power spectrum estimation framework that preserves the so-called "EoR window" and keeps track of estimator errors and covariances. We apply our method to observations from the 32-tile prototype of the Murchinson Widefield Array to demonstrate the importance of a judicious analysis technique. Lastly, we apply our method to investigate the dependence of the clean EoR window on frequency--especially the frequency dependence of the so-called "wedge" feature--and establish upper limits on the power spectrum from z = 6.2 to z = 11.7. Our lowest limit is Delta(k) < 0.3 Kelvin at 95% confidence at a comoving scale k = 0.046 Mpc^-1 and z = 9.5.

preprint2014arXiv

Radiometric Measurements of Electron Temperature and Opacity of Ionospheric Perturbations

Changes in the sky noise spectrum are used to characterize perturbations in the ionosphere. Observations were made at the same sidereal time on multiple days using a calibrated broadband dipole and radio spectrometer covering 80 to 185 MHz. In this frequency range, an ionospheric opacity perturbation changes both the electron thermal emission from the ionosphere and the absorption of the sky noise background. For the first time, these changes are confirmed to have the expected spectral signature and are used to derive the opacity and electron temperature associated with the perturbations as a function of local time. The observations were acquired at the Murchison Radio-astronomy Observatory in Western Australia from 18 April 2014 to 6 May 2014. They show perturbations that increase at sunrise, continue during the day, and decline after sunset. Magnitudes corresponding to an opacity of about 1 percent at 150 MHz with a typical electron temperature of about 800 K, were measured for the strongest perturbations.

preprint2014arXiv

Serendipitous discovery of a dying Giant Radio Galaxy associated with NGC 1534, using the Murchison Widefield Array

Recent observations with the Murchison Widefield Array at 185~MHz have serendipitously unveiled a heretofore unknown giant and relatively nearby ($z = 0.0178$) radio galaxy associated with NGC\,1534. The diffuse emission presented here is the first indication that NGC\,1534 is one of a rare class of objects (along with NGC\,5128 and NGC\,612) in which a galaxy with a prominent dust lane hosts radio emission on scales of $\sim$700\,kpc. We present details of the radio emission along with a detailed comparison with other radio galaxies with disks. NGC1534 is the lowest surface brightness radio galaxy known with an estimated scaled 1.4-GHz surface brightness of just 0.2\,mJy\,arcmin$^{-2}$. The radio lobes have one of the steepest spectral indices yet observed: $α=-2.1\pm0.1$, and the core to lobe luminosity ratio is $<0.1$\%. We estimate the space density of this low brightness (dying) phase of radio galaxy evolution as $7\times10^{-7}$\,Mpc$^{-3}$ and argue that normal AGN cannot spend more than 6\% of their lifetime in this phase if they all go through the same cycle.

preprint2014arXiv

Study of Redshifted HI from the Epoch of Reionization with Drift scan

The detection of the Epoch of Reionization (EoR) in the redshifted 21-cm line is a challenging task. Here we formulate the detection of the EoR signal using the drift scan strategy. This method potentially has better instrumental stability as compared to the case where a single patch of sky is tracked. We demonstrate that the correlation time between measured visibilities could extend up to 1-2 hr for an interferometer array such as the Murchison Widefield Array (MWA), which has a wide primary beam. We estimate the EoR power based on cross-correlation of visibilities across time and show that the drift scan strategy is capable of the detection of the EoR signal with comparable/better signal-to-noise as compared to the tracking case. We also estimate the visibility correlation for a set of bright point sources and argue that the statistical inhomogeneity of bright point sources might allow their separation from the EoR signal.

preprint2014arXiv

The Murchison Widefield Array Commissioning Survey: A Low-Frequency Catalogue of 14,110 Compact Radio Sources over 6,100 Square Degrees

We present the results of an approximately 6,100 square degree 104--196MHz radio sky survey performed with the Murchison Widefield Array during instrument commissioning between 2012 September and 2012 December: the Murchison Widefield Array Commissioning Survey (MWACS). The data were taken as meridian drift scans with two different 32-antenna sub-arrays that were available during the commissioning period. The survey covers approximately 20.5 h < Right Ascension (RA) < 8.5 h, -58 deg < Declination (Dec) < -14 deg over three frequency bands centred on 119, 150 and 180 MHz, with image resolutions of 6--3 arcmin. The catalogue has 3-arcmin angular resolution and a typical noise level of 40 mJy/beam, with reduced sensitivity near the field boundaries and bright sources. We describe the data reduction strategy, based upon mosaiced snapshots, flux density calibration and source-finding method. We present a catalogue of flux density and spectral index measurements for 14,110 sources, extracted from the mosaic, 1,247 of which are sub-components of complexes of sources.

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 study of fundamental limitations to statistical detection of redshifted HI from the epoch of reionization

In this paper we explore for the first time the relative magnitudes of three fundamental sources of uncertainty, namely, foreground contamination, thermal noise and sample variance in detecting the HI power spectrum from the Epoch of Reionization (EoR). We derive limits on the sensitivity of a Fourier synthesis telescope to detect EoR based on its array configuration and a statistical representation of images made by the instrument. We use the Murchison Widefield Array (MWA) configuration for our studies. Using a unified framework for estimating signal and noise components in the HI power spectrum, we derive an expression for and estimate the contamination from extragalactic point-like sources in three-dimensional k-space. Sensitivity for EoR HI power spectrum detection is estimated for different observing modes with MWA. With 1000 hours of observing on a single field using the 128-tile MWA, EoR detection is feasible (S/N > 1 for $k\lesssim 0.8$ Mpc$^{-1}$). Bandpass shaping and refinements to the EoR window are found to be effective in containing foreground contamination, which makes the instrument tolerant to imaging errors. We find that for a given observing time, observing many independent fields of view does not offer an advantage over a single field observation when thermal noise dominates over other uncertainties in the derived power spectrum.

preprint2013arXiv

Science with the Murchison Widefield Array

Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives.

preprint2012arXiv

A lower limit of dz > 0.06 for the duration of the reionization epoch

Observations of the 21-centimetre line of atomic hydrogen in the early Universe directly probe the history of the reionization of the gas between galaxies. The observations are challenging, though, because of the low expected signal strength (~10 mK), and contamination by strong (>100 K) foreground synchrotron emission in the Milky Way and extragalactic continuum sources2. If reionization happened rapidly, there should be a characteristic signature visible against the smooth foreground in an all-sky spectrum. Here we report an all-sky spectrum between 100 and 200 MHz, corresponding to the redshift range 6 < z < 13 for the 21-centimetre line. The data exclude a rapid reionization timescale of dz < 0.06 at the 95% confidence level.

preprint2012arXiv

Absolute calibration of a wideband antenna and spectrometer for sky noise spectral index measurements

A new method of absolute calibration of sky noise temperature using a three-position switched spectrometer, measurements of antenna and low noise amplifier impedance with a vector network analyzer, and ancillary measurements of the amplifier noise waves is described. The details of the method and its application to accurate wideband measurements of the spectral index of the sky noise are described and compared with other methods.

preprint2012arXiv

Fast Holographic Deconvolution: a new technique for precision radio interferometry

We introduce the Fast Holographic Deconvolution method for analyzing interferometric radio data. Our new method is an extension of A-projection/software-holography/forward modeling analysis techniques and shares their precision deconvolution and widefield polarimetry, while being significantly faster than current implementations that use full direction-dependent antenna gains. Using data from the MWA 32 antenna prototype, we demonstrate the effectiveness and precision of our new algorithm. Fast Holographic Deconvolution may be particularly important for upcoming 21 cm cosmology observations of the Epoch of Reionization and Dark Energy where foreground subtraction is intimately related to the precision of the data reduction.

preprint2012arXiv

Low Frequency Imaging of Fields at High Galactic Latitude with the Murchison Widefield Array 32-Element Prototype

The Murchison Widefield Array (MWA) is a new low-frequency, wide field-of-view radio interferometer under development at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. We have used a 32-element MWA prototype interferometer (MWA-32T) to observe two 50-degree diameter fields in the southern sky in the 110 MHz to 200 MHz band in order to evaluate the performance of the MWA-32T, to develop techniques for epoch of reionization experiments, and to make measurements of astronomical foregrounds. We developed a calibration and imaging pipeline for the MWA-32T, and used it to produce ~15' angular resolution maps of the two fields. We perform a blind source extraction using these confusion-limited images, and detect 655 sources at high significance with an additional 871 lower significance source candidates. We compare these sources with existing low-frequency radio surveys in order to assess the MWA-32T system performance, wide field analysis algorithms, and catalog quality. Our source catalog is found to agree well with existing low-frequency surveys in these regions of the sky and with statistical distributions of point sources derived from Northern Hemisphere surveys; it represents one of the deepest surveys to date of this sky field in the 110 MHz to 200 MHz band.

preprint2011arXiv

An MCMC approach to extracting the global 21-cm signal during the cosmic dawn from sky-averaged radio observations

Efforts are being made to observe the 21-cm signal from the 'cosmic dawn' using sky-averaged observations with individual radio dipoles. In this paper, we develop a model of the observations accounting for the 21-cm signal, foregrounds, and several major instrumental effects. Given this model, we apply Markov Chain Monte Carlo techniques to demonstrate the ability of these instruments to separate the 21-cm signal from foregrounds and quantify their ability to constrain properties of the first galaxies. For concreteness, we investigate observations between 40 and 120 MHz with the proposed DARE mission in lunar orbit, showing its potential for science return.

preprint2011arXiv

First Spectroscopic Imaging Observations of the Sun at Low Radio Frequencies with the Murchison Widefield Array Prototype

We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emission features are evident. Our data represent a significant advance in low radio frequency solar imaging, enabling us to follow the spatial, spectral, and temporal evolution of events simultaneously and in unprecedented detail. The rich variety of features seen here reaffirms the coronal diagnostic capability of low radio frequency emission and provides an early glimpse of the nature of radio observations that will become available as the next generation of low frequency radio interferometers come on-line over the next few years.

preprint2009arXiv

Foreground Contamination in Interferometric Measurements of the Redshifted 21 cm Power Spectrum

Subtraction of astrophysical foreground contamination from "dirty" sky maps produced by simulated measurements of the Murchison Widefield Array (MWA) has been performed by fitting a 3rd-order polynomial along the spectral dimension of each pixel in the data cubes. The simulations are the first to include the unavoidable instrumental effects of the frequency-dependent primary antenna beams and synthesized array beams. They recover the one-dimensional spherically-binned input redshifted 21 cm power spectrum to within approximately 1% over the scales probed most sensitively by the MWA (0.01 < k < 1 Mpc^-1) and demonstrate that realistic instrumental effects will not mask the EoR signal. We find that the weighting function used to produce the dirty sky maps from the gridded visibility measurements is important to the success of the technique. Uniform weighting of the visibility measurements produces the best results, whereas natural weighting significantly worsens the foreground subtraction by coupling structure in the density of the visibility measurements to spectral structure in the dirty sky map data cube. The extremely dense uv-coverage of the MWA was found to be advantageous for this technique and produced very good results on scales corresponding to |u| < 500 wavelengths in the uv-plane without any selective editing of the uv-coverage.

preprint2009arXiv

The Murchison Widefield Array: Design Overview

The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities.

preprint2005arXiv

A Search for Radio Gravitational Lenses, Using the Sloan Digital Sky Survey and the Very Large Array

We report on a novel search for radio gravitational lenses. Using the Very Large Array, we imaged ten candidates with both dual redshifts in Sloan Digital Sky Survey spectra and 1.4 GHz radio flux >2 mJy in the FIRST survey. The VLA maps show that in each case the radio emission is associated with the foreground galaxy rather than being lensed emission from the background galaxy, although at least four of our targets are strong lenses at optical wavelengths. These SDSS dual-redshift systems do not have lensed radio emission at the sensitivity of current radio surveys.