We report on the design, deployment, and first results from a scintillation detector deployed at the Murchison Radio-astronomy Observatory (MRO). The detector is a prototype for a larger array -- the Square Kilometre Array Particle Array (SKAPA) -- planned to allow the radio-detection of cosmic rays with the Murchison Widefield Array and the low-frequency component of the Square Kilometre Array. The prototype design has been driven by stringent limits on radio emissions at the MRO, and to ensure survivability in a desert environment. Using data taken from Nov.\ 2018 to Feb.\ 2019, we characterize the detector response while accounting for the effects of temperature fluctuations, and calibrate the sensitivity of the prototype detector to through-going muons. This verifies the feasibility of cosmic ray detection at the MRO. We then estimate the required parameters of a planned array of eight such detectors to be used to trigger radio observations by the Murchison Widefield Array.

The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $\Delta^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

The Murchison Widefield Array (MWA) is an electronically steered low frequency ($<300$\,MHz) radio interferometer, with a `slew' time less than 8seconds. Low frequency ($\sim 100$\,MHz) radio telescopes are ideally suited for rapid-response follow-up of transients due to their large field of view, the inverted spectrum of coherent emission, and the fact that the dispersion delay between a 1GHz and 100MHz pulse is on the order of $1-10$\,min for dispersion measures of $100-2000$\,pc/cm$^3$. The MWA has previously been used to provide fast follow up for transient events including gamma-ray bursts, fast radio bursts, and gravitational waves, using systems that respond to gamma-ray coordinates network (GCN) packet-based notifications. We describe a system for automatically triggering MWA observations of such events, based on VOEvent triggers, which is more flexible, capable, and accurate than previous systems. The system can respond to external multi-messenger triggers, which makes it well-suited to searching for prompt coherent radio emission from gamma-ray bursts, the study of fast radio bursts and gravitational waves, single pulse studies of pulsars, and rapid follow-up of high-energy superflares from flare stars. The new triggering system has the capability to trigger observations in both the regular correlator mode (limited to $\geq 0.5$\,s integrations) or using the Voltage Capture System (VCS, $0.1$\,ms integration) of the MWA, and represents a new mode of operation for the MWA. The upgraded standard correlator triggering capability has been in use since MWA observing semester 2018B (July-Dec 2018), and the VCS and buffered mode triggers will become available for observing in a future semester.

One of the major challenges for pulsar timing array (PTA) experiments is the mitigation of the effects of the turbulent interstellar medium (ISM) from timing data. These can potentially lead to measurable delays and/or distortions in the pulse profiles and scale strongly with the inverse of the radio frequency. Low-frequency observations are therefore highly appealing for characterizing them. However, in order to achieve the necessary time resolution to resolve profile features of short-period millisecond pulsars, phase-coherent de-dispersion is essential, especially at frequencies below $300$ MHz. We present the lowest-frequency ($80$-$220$ MHz), coherently de-dispersed detections of one of the most promising pulsars for current and future PTAs, PSR J2241$-$5236, using our new beam-former software for the MWA's voltage capture system (VCS), which reconstructs the time series at a much higher time resolution of $\sim 1 \mu$s by re-synthesizing the recorded voltage data at $10$-kHz/$100$-$\mu$s native resolutions. Our data reveal a dual-precursor type feature in the pulse profile that is either faint or absent in high-frequency observations from Parkes. The resultant high-fidelity detections have enabled dispersion measure (DM) determinations with very high precision, of the order of $(2$-$6)\times10^{-6}$ $\rm pc\,cm^{-3}$, owing to the microsecond level timing achievable for this pulsar at the MWA's low frequencies. This underscores the usefulness of low-frequency observations for probing the ISM toward PTA pulsars and informing optimal observing strategies for PTA experiments.

We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.

The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molonglo Observatory Synthesis Telescope (UTMOST) at 835 MHz, as part of an on-going observing campaign. During a 30-minute simultaneous observation, we detected the pulsar in its bright emission state for approximately 15 minutes, where 86 and 283 pulses were detected above a signal-to-noise threshold of 6 with the MWA and UTMOST, respectively. Of the detected pulses, 51 had counterparts at both frequencies and exhibited steep spectral indices for both the bright main pulse component and the precursor component. We find that the bright state pulse energy distribution is best parameterised by a log-normal distribution at both frequencies, contrary to previous results which suggested a power law distribution. Further low-frequency observations are required in order to explore in detail aspects such as pulse-to-pulse variability, intensity modulations and to better constrain the signal propagation effects due to the interstellar medium and intermittency characteristics at these frequencies. The spectral index, extended profile emission covering a large fraction of pulse longitude, and the broadband intermittency of PSR J1107-5907 suggests that future low-frequency pulsar searches, for instance those planned with SKA-Low, will be in an excellent position to find and investigate new pulsars of this type.

We present the results of a coordinated campaign conducted with the Murchison Widefield Array (MWA) to shadow Fast Radio Bursts (FRBs) detected by the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz, which resulted in simultaneous MWA observations of seven ASKAP FRBs. We de-dispersed the $24$ $\times$ $1.28$ MHz MWA images across the $170-200$ MHz band taken at 0.5 second time resolution at the known dispersion measures (DMs) and arrival times of the bursts and searched both within the ASKAP error regions (typically $\sim$ $10$ arcmin $\times$ $10$ arcmin), and beyond ($4$ deg $\times$ $4$ deg). We identified no candidates exceeding a $5\sigma$ threshold at these DMs in the dynamic spectra. These limits are inconsistent with the mean fluence scaling of $\alpha=-1.8 \pm 0.3$ (${\cal F}_\nu \propto \nu^\alpha$, where $\nu$ is the observing frequency) that is reported for ASKAP events, most notably for the three high fluence (${\cal F}_{1.4\,{\rm GHz}} \gtrsim 100$ Jy ms) FRBs 171020, 180110 and 180324. Our limits show that pulse broadening alone cannot explain our non-detections, and that there must be a spectral turnover at frequencies above 200 MHz. We discuss and constrain parameters of three remaining plausible spectral break mechanisms: free-free absorption, intrinsic spectral turn-over of the radiative processes, and magnification of signals at ASKAP frequencies by caustics or scintillation. If free-free absorption were the cause of the spectral turnover, we constrain the thickness of the absorbing medium in terms of the electron temperature, $T$, to $< 0.03$ $(T/10^4 K)^{-1.35}$ pc for FRB 171020.

We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These new tiles enhance the surface brightness sensitivity of the MWA and will improve the ability of the MWA to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ~3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u,v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of MWA continuum images. The upgrade retains all of the features that have underpinned the MWA's success (large field-of-view, snapshot image quality, pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies.

We provide the first in situ measurements of antenna element (tile) beam shapes of the Murchison Widefield Array (MWA), a low radio-frequency interferometer and an SKA precursor. Most current MWA processing pipelines use an assumed beam shape, errors in which can cause absolute and relative flux density errors, as well as polarisation 'leakage'. This makes understanding the primary beam of paramount importance, especially for sensitive experiments such as a measurement of the 21 cm line from the epoch of reionisation (EoR). The calibration requirements for measuring the EoR 21 cm line are so extreme that tile to tile beam variations may affect our ability to make a detection. Measuring the primary beam shape from visibilities alone is challenging, as multiple instrumental, atmospheric, and astrophysical factors contribute to uncertainties in the data. Building on the methods of Neben et al. (2015), we tap directly into the receiving elements of the MWA before any digitisation or correlation of the signal. Using ORBCOMM satellite passes we are able to produce all-sky maps for 4 separate tiles in the XX polarisation. We find good agreement with the cutting-edge 'fully' embedded element (FEE) model of Sokolowski et al. (2017), and observe that the MWA beamformers consistently recreate beam shapes to within ~1dB in the reliable areas of our beam maps. We also clearly observe the effects of a missing dipole from a tile in one of our beam maps, and show that the FEE model is able to reproduce this modified beam shape. We end by motivating and outlining additional onsite experiments to further constrain the primary beam behaviour.

Studying the gravitational-wave sky with pulsar timing arrays (PTAs) is a key science goal for the Square Kilometre Array (SKA) and its pathfinder telescopes. With current PTAs reaching sub-microsecond timing precision, making accurate measurements of interstellar propagation effects and mitigating them effectively has become increasingly important to realise PTA goals. As these effects are much stronger at longer wavelengths, low-frequency observations are most appealing for characterizing the interstellar medium (ISM) along the sight lines toward PTA pulsars. The Murchison Widefield Array (MWA) and the Engineering Development Array (EDA), which utilizes MWA technologies, present promising opportunities for undertaking such studies, particularly for PTA pulsars located in the southern sky. Such pulsars are also the prime targets for PTA efforts planned with the South African MeerKAT, and eventually with the SKA. In this paper we report on observations of two bright southern millisecond pulsars PSRs J0437-4715 and J2145-0750 made with these facilities; MWA observations sampling multiple frequencies across the 80-250 MHz frequency range, while the EDA providing direct-sampled baseband data to yield a large instantaneous usable bandwidth of $\sim$200 MHz. Using these observations, we investigate various aspects relating to pulsar emission and ISM properties, such as spectral evolution of the mean pulse shape, scintillation as a function of frequency, chromaticity in interstellar dispersion, and flux density spectra at low frequencies. Systematic and regular monitoring observations will help ascertain the role of low-frequency measurements in PTA experiments, while simultaneously providing a detailed characterization of the ISM toward the pulsars, which will be useful in devising optimal observing strategies for future PTA experiments.

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.

Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_\nu \propto \nu^\alpha$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteor emission and other transient events using the Murchison Widefield Array (MWA) at 72-103 MHz. In our 322-hour survey, down to a $5\sigma$ detection threshold of 3.5 Jy/beam, no transient candidates were identified as intrinsic emission from meteors. We derived an upper limit of -3.7 (95% confidence limit) on the spectral index in our frequency range. We also report detections of other transient events, like reflected FM broadcast signals from small satellites, conclusively demonstrating the ability of the MWA to detect and track space debris on scales as small as 0.1 m in low Earth orbits.

We examine data from the Murchison Widefield Array (MWA) in the frequency range 72 -- 102 MHz for a field-of-view that serendipitously contained the interstellar object 'Oumuamua on 2017 November 28. Observations took place with time resolution of 0.5 s and frequency resolution of 10 kHz. %This observation was undertaken for another purpose but due to the MWA's extremely large field-of-view, 'Oumuamua was serendipitously observed simultaneously. Based on the interesting but highly unlikely suggestion that 'Oumuamua is an interstellar spacecraft, due to some unusual orbital and morphological characteristics, we examine our data for signals that might indicate the presence of intelligent life associated with 'Oumuamua. We searched our radio data for: 1) impulsive narrow-band signals; 2) persistent narrow-band signals; and 3) impulsive broadband signals. We found no such signals with non-terrestrial origins and make estimates of the upper limits on Equivalent Isotropic Radiated Power (EIRP) for these three cases of approximately 7 kW, 840 W, and 100 kW, respectively. These transmitter powers are well within the capabilities of human technologies, and are therefore plausible for alien civilizations. While the chances of positive detection in any given Search for Extraterrestrial Intelligence (SETI) experiment are vanishingly small, the characteristics of new generation telescopes such as the MWA (and in the future, the Square Kilometre Array) make certain classes of SETI experiment easy, or even a trivial by-product of astrophysical observations. This means that the future costs of SETI experiments are very low, allowing large target lists to partially balance the low probability of a positive detection.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We describe the design and performance of the Engineering Development Array (EDA), which is a low frequency radio telescope comprising 256 dual-polarisation dipole antennas working as a phased-array. The EDA was conceived of, developed, and deployed in just 18 months via re-use of Square Kilometre Array (SKA) precursor technology and expertise, specifically from the Murchison Widefield Array (MWA) radio telescope. Using drift scans and a model for the sky brightness temperature at low frequencies, we have derived the EDA's receiver temperature as a function of frequency. The EDA is shown to be sky-noise limited over most of the frequency range measured between 60 and 240 MHz. By using the EDA in interferometric mode with the MWA, we used calibrated visibilities to measure the absolute sensitivity of the array. The measured array sensitivity matches very well with a model based on the array layout and measured receiver temperature. The results demonstrate the practicality and feasibility of using MWA-style precursor technology for SKA-scale stations. The modular architecture of the EDA allows upgrades to the array to be rolled out in a staged approach. Future improvements to the EDA include replacing the second stage beamformer with a fully digital system, and to transition to using RF-over-fibre for the signal output from first stage beamformers.

Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts.

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 $\Delta^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.

We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ~54' resolution with a sensitivity of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz.

Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction.

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.

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.

Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle \alpha_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population.

In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum.

The redshifted 21cm line of neutral hydrogen (Hi), potentially observable at low radio frequencies (~50-200 MHz), is a promising probe of the physical conditions of the inter-galactic medium during Cosmic Dawn and the Epoch of Reionisation (EoR). The sky-averaged Hi signal is expected to be extremely weak (~100 mK) in comparison to the Galactic foreground emission (~$10^4$ K). Moreover, the sky-averaged spectra measured by ground-based instruments are affected by chromatic propagation effects (of the order of tens of Kelvins) originating in the ionosphere. We analyze data collected with the upgraded BIGHORNS system deployed at the Murchison Radio-astronomy Observatory to assess the significance of ionospheric effects (absorption, emission and refraction) on the detection of the global EoR signal. We measure some properties of the ionosphere, such as the electron temperature ($T_e \approx$470 K at nighttime), magnitude, and variability of optical depth ($\tau_{100 MHz} \approx$0.01 and $\delta \tau \approx$0.005 at nighttime). According to the results of a statistical test applied on a large data sample, very long integrations lead to increased signal to noise even in the presence of ionospheric variability. This is further supported by the structure of the power spectrum of the sky temperature fluctuations, which has flicker noise characteristics at frequencies $\gtrsim 10^{-5}$ Hz, but becomes flat below $\approx 10^{-5}$ Hz. We conclude that the stochastic error introduced by the chromatic ionospheric effects tends to zero in an average. Therefore, the ionospheric effects and fluctuations are not fundamental impediments preventing ground-based instruments from integrating down to the precision required by global EoR experiments.

Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We directly examine 51 hours of data observed over 10 nights under quiet geomagnetic conditions (global storm index Kp < 2), analysing the behaviour of short-timescale angular position and peak flux density variations of around ten thousand unresolved sources. We find that while much of the variation in angular position can be attributed to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching under typical conditions. By examining bulk trends in amplitude variability, we place upper limits on the modulation index associated with ionospheric scintillation of 1-3% for the various nights. For sources fainter than ~1 Jy, this variation is below the image noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is not a significant impediment to the goals of time-domain science with the MWA at 154 MHz.

We compare first order (refractive) ionospheric effects seen by the Murchison Widefield Array (MWA) with the ionosphere as inferred from Global Positioning System (GPS) data. The first order ionosphere manifests itself as a bulk position shift of the observed sources across an MWA field of view. These effects can be computed from global ionosphere maps provided by GPS analysis centres, namely the Center for Orbit Determination in Europe (CODE), using data from globally distributed GPS receivers. However, for the more accurate local ionosphere estimates required for precision radio astronomy applications, data from local GPS networks needs to be incorporated into ionospheric modelling. For GPS observations, the ionospheric parameters are biased by GPS receiver instrument delays, among other effects, also known as receiver Differential Code Biases (DCBs). The receiver DCBs need to be estimated for any non-CODE GPS station used for ionosphere modelling, a requirement for establishing dense GPS networks in arbitrary locations in the vicinity of the MWA. In this work, single GPS station-based ionospheric modelling is performed at a time resolution of 10 minutes. Also the receiver DCBs are estimated for selected Geoscience Australia (GA) GPS receivers, located at Murchison Radio Observatory (MRO1), Yarragadee (YAR3), Mount Magnet (MTMA) and Wiluna (WILU). The ionospheric gradients estimated from GPS are compared with the ionospheric gradients inferred from radio source position shifts observed with the MWA. The ionospheric gradients at all the GPS stations show a correlation with the gradients observed with the MWA. The ionosphere estimates obtained using GPS measurements show promise in terms of providing calibration information for the MWA.

Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celestial sources are proportional to spatial gradients in the electron column density transverse to the line of sight. These can be used to probe plasma structures and waves in the ionosphere. The regional (10-100 km) scales probed by the MWA, determined by the size of its FoV and the spatial density of radio sources (typically thousands in a single FoV), complement the global (100-1000 km) scales of GPS studies and local (0.01-1 km) scales of radar scattering measurements. Our data exhibit a range of complex structures and waves. Some fluctuations have the characteristics of travelling ionospheric disturbances (TIDs), while others take the form of narrow, slowly-drifting bands aligned along the Earth's magnetic field.

We have carried out multiwavelength observations of the near-by ($z=0.046$) rich, merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the eastern and western peripheries of the cluster, were detected at all the frequencies. The relics, with a linear extent of $\sim$ 1 Mpc each, are separated by $\sim$ 2 Mpc. Combining the current observations with those in the literature, we have obtained the spectra of these relics over the frequency range 80 -- 1400 MHz. The spectra follow power laws, with $\alpha$ = $-1.17\pm0.06$ and $-1.37\pm0.08$ for the west and east relics, respectively ($S \propto \nu^{\alpha}$). Assuming the break frequency to be near the lower end of the spectrum we estimate the age of the relics to be $\sim$ 0.4 Gyr. No diffuse radio emission from the central regions of the cluster (halo) was detected. The upper limit on the radio power of any possible halo that might be present in the cluster is a factor of 35 lower than that expected from the radio power and X-ray luminosity correlation for cluster halos. From this we conclude that the cluster halo is very extended ($>$ 500 kpc) and/or most of the radio emission from the halo has decayed. The current limit on the halo radio power is a factor of ten lower than the existing upper limits with possible implications for models of halo formation.

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40 kHz frequency resolution and 0.5 s time resolution; the second year using 10 kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing and image weighting scheme. At 154 MHz the image resolution is approximately 2.5 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA's primary beam.

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.

An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements.

The science cases for incorporating high time resolution capabilities into modern radio telescopes are as numerous as they are compelling. Science targets range from exotic sources such as pulsars, to our Sun, to recently detected possible extragalactic bursts of radio emission, the so-called fast radio bursts (FRBs). Originally conceived purely as an imaging telescope, the initial design of the Murchison Widefield Array (MWA) did not include the ability to access high time and frequency resolution voltage data. However, the flexibility of the MWA's software correlator allowed an off-the-shelf solution for adding this capability. This paper describes the system that records the 100 micro-second and 10 kHz resolution voltage data from the MWA. Example science applications, where this capability is critical, are presented, as well as accompanying commissioning results from this mode to demonstrate verification.

The Murchison Widefield Array (MWA) is a Square Kilometre Array (SKA) Precursor. The telescope is located at the Murchison Radio--astronomy Observatory (MRO) in Western Australia (WA). The MWA consists of 4096 dipoles arranged into 128 dual polarisation aperture arrays forming a connected element interferometer that cross-correlates signals from all 256 inputs. A hybrid approach to the correlation task is employed, with some processing stages being performed by bespoke hardware, based on Field Programmable Gate Arrays (FPGAs), and others by Graphics Processing Units (GPUs) housed in general purpose rack mounted servers. The correlation capability required is approximately 8 TFLOPS (Tera FLoating point Operations Per Second). The MWA has commenced operations and the correlator is generating 8.3 TB/day of correlation products, that are subsequently transferred 700 km from the MRO to Perth (WA) in real-time for storage and offline processing. In this paper we outline the correlator design, signal path, and processing elements and present the data format for the internal and external interfaces.

We present an in-situ antenna characterization method and results for a "low-frequency" radio astronomy engineering prototype array, characterized over the 75-300 MHz frequency range. The presence of multiple cosmic radio sources, particularly the dominant Galactic noise, makes in-situ characterization at these frequencies challenging; however, it will be shown that high quality measurement is possible via radio interferometry techniques. This method is well-known in the radio astronomy community but seems less so in antenna measurement and wireless communications communities, although the measurement challenges involving multiple undesired sources in the antenna field-of-view bear some similarities. We discuss this approach and our results with the expectation that this principle may find greater application in related fields.

Using the new wideband capabilities of the Australia Telescope Compact Array (ATCA), we obtain spectra for PKS 1718-649, a well-known gigahertz-peaked spectrum radio source. The observations, between approximately 1 and 10 GHz over three epochs spanning approximately 21 months, reveal variability both above the spectral peak at ~3 GHz and below the peak. The combination of the low and high frequency variability cannot be easily explained using a single absorption mechanism, such as free-free absorption or synchrotron self-absorption. We find that the PKS 1718-649 spectrum and its variability are best explained by variations in the free-free optical depth on our line-of-sight to the radio source at low frequencies (below the spectral peak) and the adiabatic expansion of the radio source itself at high frequencies (above the spectral peak). The optical depth variations are found to be plausible when X-ray continuum absorption variability seen in samples of Active Galactic Nuclei is considered. We find that the cause of the peaked spectrum in PKS 1718-649 is most likely due to free-free absorption. In agreement with previous studies, we find that the spectrum at each epoch of observation is best fit by a free-free absorption model characterised by a power-law distribution of free-free absorbing clouds. This agreement is extended to frequencies below the 1 GHz lower limit of the ATCA by considering new observations with Parkes at 725 MHz and 199 MHz observations with the newly operational Murchison Widefield Array. These lower frequency observations argue against families of absorption models (both free-free and synchrotron self-absorption) that are based on simple homogenous structures.

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: $\alpha=-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.

We present new low-frequency observations of the nearby radio galaxy Fornax A at 154 MHz with the Murchison Widefield Array, microwave flux-density measurements obtained from WMAP and Planck data, and gamma-ray flux densities obtained from Fermi data. We also compile a comprehensive list of previously published images and flux-density measurements at radio, microwave and X-ray energies. A detailed analysis of the spectrum of Fornax A between 154 MHz and 1510 MHz reveals that both radio lobes have a similar spatially-averaged spectral index, and that there exists a steep-spectrum bridge of diffuse emission between the lobes. Taking the spectral index of both lobes to be the same, we model the spectral energy distribution of Fornax A across an energy range spanning eighteen orders of magnitude, to investigate the origin of the X-ray and gamma-ray emission. A standard leptonic model for the production of both the X-rays and gamma-rays by inverse-Compton scattering does not fit the multi-wavelength observations. Our results best support a scenario where the X-rays are produced by inverse-Compton scattering and the gamma-rays are produced primarily by hadronic processes confined to the filamentary structures of the Fornax A lobes.

We present the results of a survey for low frequency radio emission from 17 known exoplanetary systems with the Murchison Widefield Array. This sample includes 13 systems that have not previously been targeted with radio observations. We detected no radio emission at 154 MHz, and put 3 sigma upper limits in the range 15.2-112.5 mJy on this emission. We also searched for circularly polarised emission and made no detections, obtaining 3 sigma upper limits in the range 3.4-49.9 mJy. These are comparable with the best low frequency radio limits in the existing literature and translate to luminosity limits of between 1.2 x 10^14 W and 1.4 x 10^17 W if the emission is assumed to be 100% circularly polarised. These are the first results from a larger program to systematically search for exoplanetary emission with the MWA.

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.

The Murchison Widefield Array (MWA) is a new low frequency interferometric radio telescope, operating in the remote Murchison Radio Observatory in Western Australia. In this paper we present the first MWA observations of the well known radio relics in Abell 3667 (A3667) between 120 and 226 MHz. We clearly detect the radio relics in A3667 and present flux estimates and spectral indices for these features. The average spectral index of the north-west (NW) and south-east (SE) relics is -0.9 +/- 0.1 between 120 and 1400 MHz. We are able to resolve spatial variation in the spectral index of the NW relic from -1.7 to -0.4, which is consistent with results found at higher frequencies.

We present the first Murchison Widefield Array observations of the well-known cluster of galaxies Abell 3667 (A3667) between 105 and 241 MHz. A3667 is one of the best known examples of a galaxy cluster hosting a double radio relic and has been reported to contain a faint radio halo and bridge. The origins of radio halos, relics and bridges is still unclear, however galaxy cluster mergers seems to be an important factor. We clearly detect the North-West (NW) and South-East (SE) radio relics in A3667 and find an integrated flux density at 149 MHz of 28.1 +/- 1.7 and 2.4 +/- 0.1 Jy, respectively, with an average spectral index, between 120 and 1400 MHz, of -0.9 +/- 0.1 for both relics. We find evidence of a spatial variation in the spectral index across the NW relic steepening towards the centre of the cluster, which indicates an ageing electron population. These properties are consistent with higher frequency observations. We detect emission that could be associated with a radio halo and bridge. How- ever, due to the presence of poorly sampled large-scale Galactic emission and blended point sources we are unable to verify the exact nature of these features.

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.

Astronomical widefield imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new widefield interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependencies of CASA's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarisation correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2-12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square-Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.

The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged is not effective for solar imaging, because many of the features of interest are short lived and change dramatically over small fractional bandwidths. Building on the advances in radio frequency technology, digital signal processing and computing, the kind of instruments needed to simultaneously capture the evolution of solar emission in time, frequency, morphology and polarization over a large spectral span with the requisite imaging fidelity, and time and frequency resolution have only recently begun to appear. Of this class of instruments, the Murchison Widefield Array (MWA) is best suited for solar observations. The MWA has now entered a routine observing phase and here we present some early examples from MWA observations.

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.

We present a search for transient and variable radio sources at 154 MHz with the Murchison Widefield Array 32-tile prototype. Fifty-one images were obtained that cover a field of view of 1430 deg^2 centred on Hydra A. The observations were obtained over three days in 2010 March and three days in 2011 April and May. The mean cadence of the observations was 26 minutes and there was additional temporal information on day and year timescales. We explore the variability of a sample of 105 low frequency radio sources within the field. Four bright (S > 6 Jy) candidate variable radio sources were identified that displayed low levels of short timescale variability (26 minutes). We conclude that this variability is likely caused by simplifications in the calibration strategy or ionospheric effects. On the timescale of one year we find two sources that show significant variability. We attribute this variability to either refractive scintillation or intrinsic variability. No radio transients were identified and we place an upper limit on the surface density of sources rho < 7.5 x 10^-5 deg^-2 with flux densities > 5.5 Jy, and characteristic timescales of both 26 minutes and one year.

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.

We present new wide-field observations of Centaurus A (Cen A) and the surrounding region at 118 MHz with the Murchison Widefield Array (MWA) 32-tile prototype, with which we investigate the spectral-index distribution of Cen A's giant radio lobes. We compare our images to 1.4 GHz maps of Cen A and compute spectral indices using temperature-temperature plots and spectral tomography. We find that the morphologies at 118 MHz and 1.4 GHz match very closely apart from an extra peak in the southern lobe at 118 MHz, which provides tentative evidence for the existence of a southern counterpart to the northern middle lobe of Cen A. Our spatially-averaged spectral indices for both the northern and southern lobes are consistent with previous analyses, however we find significant spatial variation of the spectra across the extent of each lobe. Both the spectral-index distribution and the morphology at low radio frequencies support a scenario of multiple outbursts of activity from the central engine. Our results are consistent with inverse-Compton modelling of radio and gamma-ray data that supports a value for the lobe age of between 10 and 80 Myr.

The Murchison Widefield Array (MWA) is a new low frequency interferomeric radio telescope. The MWA is the low frequency precursor to the Square Kilometre Array (SKA) and is the first of three SKA precursors to be operational, supporting a varied science mission ranging from the attempted detection of the Epoch of Reionisation to the monitoring of solar flares and space weather. We explore the possibility that the MWA can be used for the purposes of Space Situational Awareness (SSA). In particular we propose that the MWA can be used as an element of a passive radar facility operating in the frequency range 87.5 - 108 MHz (the commercial FM broadcast band). In this scenario the MWA can be considered the receiving element in a bi-static radar configuration, with FM broadcast stations serving as non-cooperative transmitters. The FM broadcasts propagate into space, are reflected off debris in Earth orbit, and are received at the MWA. The imaging capabilities of the MWA can be used to simultaneously detect multiple pieces of space debris, image their positions on the sky as a function of time, and provide tracking data that can be used to determine orbital parameters. Such a capability would be a valuable addition to Australian and global SSA assets, in terms of southern and eastern hemispheric coverage. We provide a feasibility assessment of this proposal, based on simple calculations and electromagnetic simulations that shows the detection of sub-metre size debris should be possible (debris radius of >0.5 m to ~1000 km altitude). We also present a proof-of-concept set of observations that demonstrate the feasibility of the proposal, based on the detection and tracking of the International Space Station via reflected FM broadcast signals originating in south-west Western Australia. These observations broadly validate our calculations and simulations.

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 \pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(\alpha,\delta) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

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.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements (known as tiles) distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper concentrates on the capabilities of the MWA for solar science and summarises some of the solar science results to date, in advance of the initial operation of the final instrument in 2013.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper will focus on a brief discussion of the as-built MWA system, highlighting several novel characteristics of the instrument, and a brief progress report (as of June 2012) on the final construction phase. Practical completion of the MWA is expected in November 2012, with commissioning commencing from approximately August 2012 and operations commencing near mid 2013. A brief description of recent science results from the MWA prototype instrument is given.

A new generation of low frequency radio telescopes is seeking to observe the redshifted 21 cm signal from the Epoch of Reionization (EoR), requiring innovative methods of calibration and imaging to overcome the difficulties of widefield low frequency radio interferometry. Precise calibration will be required to separate the small expected EoR signal from the strong foreground emission at the frequencies of interest between 80 and 300 MHz. The Moon may be useful as a calibration source for detection of the EoR signature, as it should have a smooth and predictable thermal spectrum across the frequency band of interest. Initial observations of the Moon with the Murchison Widefield Array 32 tile prototype show that the Moon does exhibit a similar trend to that expected for a cool thermally emitting body in the observed frequency range, but that the spectrum is corrupted by reflected radio emission from Earth. In particular, there is an abrupt increase in the observed flux density of the Moon within the internationally recognised Frequency Modulated (FM) radio band. The observations have implications for future low frequency surveys and EoR detection experiments that will need to take this reflected emission from the Moon into account. The results also allow us to estimate the equivalent isotropic power emitted by the Earth in the FM band and to determine how bright the Earth might appear at metre wavelengths to an observer beyond our own solar system.

Using the final 128 antenna locations of the Murchison Widefield Array (MWA), we calculate its sensitivity to the Epoch of Reionization (EoR) power spectrum of red- shifted 21 cm emission for a fiducial model and provide the tools to calculate the sensitivity for any model. Our calculation takes into account synthesis rotation, chro- matic and asymmetrical baseline effects, and excludes modes that will be contaminated by foreground subtraction. For the fiducial model, the MWA will be capable of a 14{\sigma} detection of the EoR signal with one full season of observation on two fields (900 and 700 hours).

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.

The Murchison Widefield Array (MWA) is one of three Square Kilometre Array Precursor telescopes and is located at the Murchison Radio-astronomy Observatory in the Murchison Shire of the mid-west of Western Australia, a location chosen for its extremely low levels of radio frequency interference. The MWA operates at low radio frequencies, 80-300 MHz, with a processed bandwidth of 30.72 MHz for both linear polarisations, and consists of 128 aperture arrays (known as tiles) distributed over a ~3 km diameter area. Novel hybrid hardware/software correlation and a real-time imaging and calibration systems comprise the MWA signal processing backend. In this paper the as-built MWA is described both at a system and sub-system level, the expected performance of the array is presented, and the science goals of the instrument are summarised.

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.

Antenna layout is an important design consideration for radio interferometers because it determines the quality of the snapshot point spread function (PSF, or array beam). This is particularly true for experiments targeting the 21 cm Epoch of Reionization signal as the quality of the foreground subtraction depends directly on the spatial dynamic range and thus the smoothness of the baseline distribution. Nearly all sites have constraints on where antennas can be placed---even at the remote Australian location of the MWA (Murchison Widefield Array) there are rock outcrops, flood zones, heritages areas, emergency runways and trees. These exclusion areas can introduce spatial structure into the baseline distribution that enhance the PSF sidelobes and reduce the angular dynamic range. In this paper we present a new method of constrained antenna placement that reduces the spatial structure in the baseline distribution. This method not only outperforms random placement algorithms that avoid exclusion zones, but surprisingly outperforms random placement algorithms without constraints to provide what we believe are the smoothest constrained baseline distributions developed to date. We use our new algorithm to determine antenna placements for the originally planned MWA, and present the antenna locations, baseline distribution, and snapshot PSF for this array choice.

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.

The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 MHz to 300 MHz with an instantaneous field of view that is tens of degrees wide and a resolution of a few arcminutes. A 32-station prototype of the MWA has been recently commissioned and a set of observations taken that exercise the whole acquisition and processing pipeline. We present Stokes I, Q, and U images from two ~4 hour integrations of a field 20 degrees wide centered on Pictoris A. These images demonstrate the capacity and stability of a real-time calibration and imaging technique employing the weighted addition of warped snapshots to counter extreme wide field imaging distortions.

It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.