Search

We present the first of two papers dedicated to verifying the Australian Epoch of Reionisation pipeline (AusEoRPipe) through simulation. The AusEoRPipe aims to disentangle 21-cm radiation emitted by gas surrounding the very first stars from contaminating foreground astrophysical sources, and has been in development for close to a decade. In this paper, we build an accurate 21-cm sky model that can be used by the WODEN simulation software to create visibilities containing a predictable 21-cm signal. We verify that the power spectrum estimator CHIPS can recover this signal in the absence of foregrounds. We also investigate how measurements in Fourier-space are correlated, and how their gridded density affects the power spectrum. We measure and fit for this effect using Gaussian-noise simulations of the MWA phase I layout. We find a gridding density correction factor of 2.651 appropriate for integrations equal to or greater than 30 minutes of data, which contain observations with multiple primary beam pointings and LSTs. Paper II of this series will use the results of this paper to test the AusEoRPipe in the presence of foregrounds and instrumental effects., Comment: 13 pages, 11 Figures, 1 table, accepted in the Publications of the Astronomical Society of Australia (PASA)

The cosmological 21cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially-dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal however is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heating and reionization, by calculating the skew spectrum (SS) of the 21cm signal using MERAXES, which couples a semi-analytic galaxy population with semi-numerical reionization simulations. The SS is the cross-spectrum of the quadratic temperature brightness field with itself. We generate a set of seven simulations from $z = 30$ to $z = 5$, varying the halo mass threshold for hosting star-formation, the X-ray luminosity per star-formation rate, and the minimum X-ray energy escaping host galaxies. We find the SS is predominantly negative as a function of redshift, transitioning to positive towards the start of reionization, and peaking during the midpoint of reionization. We do not see a negative dip in the SS during reionization, likely due to the specifics of modelling ionization sources. We normalise the SS by the PS during reionization isolating the non-Gaussianities. We find a trough ($k\sim\,0.1\,\textrm{Mpc}^{-1}$) and peak ($k\sim\,0.4-1\,\textrm{Mpc}^{-1}$) in the normalised SS during the mid to late periods of reionization. These correlate to the ionization topology, and neutral islands in the IGM. We calculate the cosmic variance of the normalised SS, and find these features are detectable in the absence of foregrounds with the SKA_LOW., Comment: 19 pages, 14 figures

We present deep upper limits from the 2014 Murchison Widefield Array (MWA) Phase I observing season, with a particular emphasis on identifying the spectral fingerprints of extremely faint radio frequency interference (RFI) contamination in the 21~cm power spectra (PS). After meticulous RFI excision involving a combination of the \textsc{SSINS} RFI flagger and a series of PS-based jackknife tests, our lowest upper limit on the Epoch of Reionization (EoR) 21~cm PS signal is $\Delta^2 \leq 1.61\cdot10^4 \text{ mK}^2$ at $k=0.258\text{ h Mpc}^{-1}$ at a redshift of 7.1 using 14.7 hours of data. By leveraging our understanding of how even fainter RFI is likely to contaminate the EoR PS, we are able to identify ultra-faint RFI signals in the cylindrical PS. Surprisingly this signature is most obvious in PS formed with less than an hour of data, but is potentially subdominant to other systematics in multiple-hour integrations. Since the total RFI budget in a PS detection is quite strict, this nontrivial integration behavior suggests a need to more realistically model coherently integrated ultra-faint RFI in PS measurements so that its potential contribution to a future detection can be diagnosed., Comment: Update acknowledgements, author metadata, and attach journal doi

Temporal broadening is a commonly observed property of fast radio bursts (FRBs), associated with turbulent media which cause radiowave scattering. Similarly to dispersion, scattering is an important probe of the media along the line of sight to an FRB source, such as the circum-burst or circum-galactic mediums (CGM). Measurements of characteristic scattering times alone are insufficient to constrain the position of the dominant scattering media along the line of sight. However, where more than one scattering screen exists, Galactic scintillation can be leveraged to form strong constraints. We quantify the scattering and scintillation in 10 FRBs with 1) known host galaxies and redshifts and 2) captured voltage data enabling high-time resolution analysis. We find strong evidence for two screens in three cases. For FRBs 20190608B and 20210320C, we find evidence for scattering screens less than approximately 16.7 and 3000 kpc respectively, from their sources, consistent with the scattering occurring in the circum-burst environment, the host ISM (inter-stellar medium) or the CGM. For FRB 20201124A we find a low modulation index that evolves over the burst's scattering tail, indicating the presence of a scattering screen $\approx9$ kpc from the host, and excluding the circum-burst environment from potential scattering sites. By assuming that pulse broadening is contributed by the host galaxy ISM or circum-burst environment, the lack of observed scintillation in four FRBs in our sample suggests that existing models may be poor estimators of scattering times associated with the Milky Way's ISM, similar to the anomalously low scattering observed for FRB 20201124A., Comment: 13 pages, 9 figures, 1 table. Accepted in MNRAS

The redshifted cosmological 21 cm signal emitted by neutral hydrogen during the first billion years of the universe is much fainter relative to other galactic and extragalactic radio emissions, posing a great challenge towards detection of the signal. Therefore, precise instrumental calibration is a vital prerequisite for the success of radio interferometers such as the Murchison Widefield Array (MWA), which aim for a 21 cm detection. Over the previous years, novel calibration techniques targeting the power spectrum paradigm of EoR science have been actively researched and where possible implemented. Using recently acquired computation resources for the MWA, we test the full capabilities of the state-of-the-art calibration techniques available for the MWA EoR project, with a focus on both direction dependent and direction independent calibration. Specifically, we investigate improvements that can be made in the vital calibration stages of sky modelling, ionospheric correction, and compact source foreground subtraction as applied in the hybrid foreground mitigation approach (one that combines both foreground subtraction and avoidance). Additionally, we investigate a method of ionospheric correction using interpolated ionospheric phase screens and assess its performance in the power spectrum space. Overall, we identify a refined RTS calibration configuration that leads to an at least 2 factor reduction of the EoR window power contamination at the $0.1 \; \text{hMpc}^{-1}$ scale. The improvement marks a step further towards detecting the 21 cm signal using the MWA and the forthcoming SKA low telescope., Comment: Accepted for publication by Publications of the Astronomical Society of Australia (PASA)

We investigate the contribution of extended radio sources such as Centaurus A, and Galactic supernova remnants (SNRs) to our ability to detect the statistical $21\,\rm{cm}$ signal from the Epoch of Reionisation (EoR) with the Murchison Widefield Array (MWA). These sources are typically ignored because they are in highly attenuated parts of the MWA primary beam, however in aggregate these sources have apparent flux densities of $10\,\rm{Jy}$ on angular scales we expect to detect the $21\,\rm{cm}$ signal. We create bespoke multi-component 2D Gaussian models for Galactic SNRs and for Centaurus A, and simulate the visibilities for two MWA snapshot observations. We grid those visibilities and then Fourier transform them with respect to frequency, averaging them both spherically and cylindrically to produce the 1D and 2D power spectra. We compare the simulated 1D power spectra to the expected $21\,\rm{cm}$ power spectrum. We find that although these extended sources are in highly attenuated parts of the MWA primary beam pattern, collectively they have enough power ($\sim10^4-10^5\,\rm{mK^2}\,\it{h^{-3}} \,\rm{Mpc^{3}}$) on EoR significant modes ($|\mathbf{k}| \leq 0.1 h \rm{Mpc}^{-1}$) to prohibit detection of the $21\,\rm{cm}$ signal ($10^4\,\rm{mK^2}\,\it{h^{-3}} \,\rm{Mpc^{3}}$). We find that $50-90\%$ of sources must be removed in order to reduce leakage to a level of $10-20\%$ of the $21\,\rm{cm}$ power spectrum on EoR significant modes. The effects of widefield extended sources will have implications on the detectability of the $21\,\rm{cm}$ signal for the MWA and with the future Square Kilometre Array (SKA)., Comment: 17 pages, 8 png figures

One of the principal systematic constraints on the Epoch of Reionisation (EoR) experiment is the accuracy of the foreground calibration model. Recent results have shown that highly accurate models of extended foreground sources, and including models for sources in both the primary beam and its sidelobes, are necessary for reducing foreground power. To improve the accuracy of the source models for the EoR fields observed by the Murchison Widefield Array (MWA), we conducted the MWA Long Baseline Epoch of Reionisation Survey (LoBES). This survey consists of multi-frequency observations of the main MWA EoR fields and their eight neighbouring fields using the MWA Phase II extended array. We present the results of the first half of this survey centred on the MWA EoR0 observing field (centred at RA(J2000) 0 h, Dec(J2000) -27 deg). This half of the survey covers an area of 3069 degrees$^2$, with an average rms of 2.1 mJy beam$^{-1}$. The resulting catalogue contains a total of 80824 sources, with 16 separate spectral measurements between 100 and 230 MHz, and spectral modelling for 78$\%$ of these sources. Over this region we estimate that the catalogue is 90$\%$ complete at 32 mJy, and 70$\%$ complete at 10.5~mJy. The overall normalised source counts are found to be in good agreement with previous low-frequency surveys at similar sensitivities. Testing the performance of the new source models we measure lower residual rms values for peeled sources, particularly for extended sources, in a set of MWA Phase I data. The 2-dimensional power spectrum of these data residuals also show improvement on small angular scales -- consistent with the better angular resolution of the LoBES catalogue. It is clear that the LoBES sky models improve upon the current sky model used by the Australian MWA EoR group for the EoR0 field., Comment: 25 pages, 14 figures, Accepted for publication in PASA

Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is instructive to explore different aspects of existing datasets. One such aspect is the choice of observing field. To date, MWA EoR observations have largely focused on the EoR0 field. In this work, we present a new detailed analysis of the EoR1 field. The EoR1 field is one of the coldest regions of the Southern radio sky, but contains the very bright radio galaxy Fornax-A. The presence of this bright extended source in the primary beam of the interferometer makes the calibration and analysis of EoR1 particularly challenging. We demonstrate the effectiveness of a recently developed shapelet model of Fornax-A in improving the results from this field. We also describe and apply a series of data quality metrics which identify and remove systematically contaminated data. With substantially improved source models, upgraded analysis algorithms and enhanced data quality metrics, we determine EoR power spectrum upper limits based on analysis of the best $\sim$14-hours data observed during 2015 and 2014 at redshifts 6.5, 6.8 and 7.1, with the lowest $2\sigma$ upper limit at z=6.5 of $\Delta^2 \leq (73.78 ~\mathrm{mK)^2}$ at $k=0.13~\mathrm{h~ Mpc^{-1}}$, improving on previous EoR1 measurement results., Comment: 18 pages, 17 figures, Accepted for publication in MNRAS

We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535-571, which underwent a prolonged outburst beginning on 2 September 2017. We monitored MAXI J1535-571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 20 September to 14 October 2017. The source was quasi-simultaneously observed over the frequency range 0.84-19 GHz by UTMOST (the upgraded Molonglo Observatory Synthesis Telescope), the Australian Square Kilometre Array Pathfinder, the Australia Telescope Compact Array (ATCA), and the Australian Long Baseline Array (LBA). Using the LBA observations from 23 September 2017, we measured the source size to be $34\pm1$ mas. During the brightest radio flare on 21 September 2017, the source was detected down to 119 MHz by the MWA, and the radio spectrum indicates a turnover between 250 and 500 MHz, which is most likely due to synchrotron self-absorption (SSA). By fitting the radio spectrum with a SSA model and using the LBA size measurement, we determined various physical parameters of the jet knot (identified in ATCA data), including the jet opening angle (= $4.5\pm1.2^{\circ}$) and the magnetic field strength (= $104^{+80}_{-78}$ mG). Our fitted magnetic field strength agrees reasonably well with that inferred from the standard equipartition approach, suggesting the jet knot to be close to equipartition. Our study highlights the capabilities of the Australian suite of radio telescopes to jointly probe radio jets in black hole X-ray binaries (BH-XRBs) via simultaneous observations over a broad frequency range, and with differing angular resolutions. This suite allows us to determine the physical properties of XRB jets. Finally, our study emphasizes the potential contributions that can be made by the low-frequency part of the Square Kilometre Array (SKA-Low) in the study of BH-XRBs., Comment: This article has been accepted for publication in PASA. This article has 16 pages, 7 figures and 1 table

Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration., Comment: 17 pages, 18 figures, accepted to MNRAS

The Epoch of Reionisation (EoR) is the period within which the neutral universe transitioned to an ionised one. This period remains unobserved using low-frequency radio interferometers which target the 21 cm signal of neutral hydrogen emitted in this era. The Murchison Widefield Array (MWA) radio telescope was built with the detection of this signal as one of its major science goals. One of the most significant challenges towards a successful detection is that of calibration, especially in the presence of the Earth's ionosphere. By introducing refractive source shifts, distorting source shapes and scintillating flux densities, the ionosphere is a major nuisance in low-frequency radio astronomy. We introduce SIVIO, a software tool developed for simulating observations of the MWA through different ionospheric conditions estimated using thin screen approximation models and propagated into the visibilities. This enables us to directly assess the impact of the ionosphere on observed EoR data and the resulting power spectra. We show that the simulated data captures the dispersive behaviour of ionospheric effects. We show that the spatial structure of the simulated ionospheric media is accurately reconstructed either from the resultant source positional offsets or from parameters evaluated during the data calibration procedure. In turn, this will inform on the best strategies of identifying and efficiently eliminating ionospheric contamination in EoR data moving into the Square Kilometre Array era., Comment: Accepted for publication in Publications of the Astronomical Society of Australia (PASA). 10 pages, 9 figures, 1 table

Aiming to fill a crucial gap in our observational knowledge of the early Universe, experiments around the world continue to attempt to verify the claimed detection of the redshifted 21-cm signal from Cosmic Dawn by the EDGES experiment. This sky-averaged or 'global' signal from neutral hydrogen should be detectable at low radio frequencies (50-200 MHz), but is difficult to measure due to bright foreground emission and difficulties in reaching the required levels of instrumental-calibration precision. In this paper we outline our progress toward using a novel new method to measure the global redshifted 21-cm signal. Motivated by the need to use alternative methods with very different systematic errors to EDGES for an independent result, we employ an array of closely-spaced antennas to measure the global sky signal interferometrically, rather than using the conventional approach with a single antenna. We use simulations to demonstrate our newly-developed methods and show that, for an idealised instrument, a 21-cm signal could theoretically be extracted from the visibilities of an array of closely-spaced dipoles. We verify that our signal-extraction methods work on real data using observations made with a Square Kilometre-Array-like prototype; the Engineering Development Array-2. Finally, we use the lessons learned in both our simulations and observations to lay out a clear plan for future work, which will ultimately lead to a new global redshifted 21-cm instrument: the All-Sky SignAl Short-Spacing INterferometer (ASSASSIN)., Comment: 18 pages, 13 figures

To make a power spectrum (PS) detection of the 21 cm signal from the Epoch of Reionisation (EoR), one must avoid/subtract bright foreground sources. Sources such as Fornax A present a modelling challenge due to spatial structures spanning from arc seconds up to a degree. We compare modelling with multi-scale (MS) CLEAN components to 'shapelets', an alternative set of basis functions. We introduce a new image-based shapelet modelling package, SHAMFI. We also introduce a new CUDA simulation code (WODEN) to generate point source, Gaussian, and shapelet components into visibilities. We test performance by modelling a simulation of Fornax A, peeling the model from simulated visibilities, and producing a residual PS. We find the shapelet method consistently subtracts large-angular-scale emission well, even when the angular-resolution of the data is changed. We find that when increasing the angular-resolution of the data, the MS CLEAN model worsens at large angular-scales. When testing on real MWA data, the expected improvement is not seen in real data because of the other dominating systematics still present. Through further simulation we find the expected differences to be lower than obtainable through current processing pipelines. We conclude shapelets are worthwhile for subtracting extended galaxies, and may prove essential for an EoR detection in the future, once other systematics have been addressed., Comment: 17 pages, 11 Figures, accepted for publication in Publications of the Astronomical Society of Australia (18/05/2020). "For the SHAMFI code, see: https://github.com/JLBLine/SHAMFI" . "For the SHAMFI documentation, see: https://shamfi.readthedocs.io/" . "For the WODEN code and documentation see: https://github.com/JLBLine/WODEN"

Detecting a signal from the Epoch of Reionisation (EoR) requires an exquisite understanding of galactic and extra-galactic foregrounds, low frequency radio instruments, instrumental calibration, and data analysis pipelines. In this work we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibration errors have the potential to inhibit EoR detections by introducing additional spectral features that mimic the structure of EoR signals. We present a straightforward way to estimate the impact of a wide variety of modelling residuals in EoR PS estimation. We apply this framework to the specific case of broken dipoles in Murchison Widefield Array (MWA) to understand its effect and estimate its impact on PS estimation. Combining an estimate of the percentage of MWA tiles that have at least one broken dipole (15%-40%) with an analytic description of beam errors induced by such dipoles, we compute the residuals of the foregrounds after calibration and source subtraction. We find that that incorrect beam modelling introduces bias in the 2D-PS on the order of $\sim 10^3\, \mathrm{mK}^2 \,h^{-3}\, \mathrm{Mpc}^{3}$. Although this is three orders of magnitude lower than current lowest limits, it is two orders of magnitude higher than the expected signal. Determining the accuracy of both current beam models and direction dependent calibration pipelines is therefore crucial in our search for an EoR signal., Comment: 14 pages, 8 figures. Accepted for publication in MNRAS December 2019. Erratum submitted on 20 February 2020

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}}$ $

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., Comment: 38 pages, 12 figures, accepted for publication in PASA

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., Comment: 14 pages, 5 figures, accepted for publication in PASA

Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $\Delta^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude., Comment: 21 pages, 10 figures, published in ApJ, updated title and text

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., Comment: 09 pages, 4 figures, 2 tables

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., Comment: 19 pages, 7 figures, 4 tables, accepted for publication in ApJ; minor updates to keep in-line with published version

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., Comment: 9 pages, 3 figures , 2 tables (accepted for publication in ApJL)

We present early results from a project to measure the sky-averaged (global), redshifted $21\,$cm signal from the Epoch of Reionisation (EoR), using the Murchison Widefield Array (MWA) telescope. Because interferometers are not sensitive to a spatially-invariant global average, they cannot be used to detect this signal using standard techniques. However, lunar occultation of the radio sky imprints a spatial structure on the global signal, allowing us to measure the average brightness temperature of the patch of sky immediately surrounding the Moon. In this paper we present one night of Moon observations with the MWA between 72 - 230 MHz and verify our techniques to extract the background sky temperature from measurements of the Moon's flux density. We improve upon previous work using the lunar occultation technique by using a more sophisticated model for reflected `earthshine' and by employing image differencing to remove imaging artefacts. We leave the Moon's (constant) radio brightness temperature as a free parameter in our fit to the data and as a result, measure $T_{\rm{moon}} = 180 \pm 12 $ K and a Galactic synchrotron spectral index of $-2.64\pm0.14$, at the position of the Moon. Finally, we evaluate the prospects of the lunar occultation technique for a global EoR detection and map out a way forward for future work with the MWA., Comment: 13 pages, 7 figures, accepted in MNRAS

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., Comment: 13 pages, 7 figures Accepted to PASA

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., Comment: 23 pages, 10 Figures, 4 Tables, Accepted to ApJS

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., Comment: 20 pages, 11 figures. Accepted to ApJ

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., Comment: Accepted for Publication in MNRAS on 11 April 2018, 11 pages, 11 figures, 5 tables

Observation of the 21cm line signal from neutral hydrogen during the Epoch of Reionization is challenging due to extremely bright Galactic and extragalactic foregrounds and complicated instrumental calibration. A reasonable approach for mitigating these problems is the cross correlation with other observables. In this work, we present the first results of the cross power spectrum (CPS) between radio images observed by the Murchison Widefield Array and the cosmic microwave background (CMB), measured by the Planck experiment. We study the systematics due to the ionospheric activity, the dependence of CPS on group of pointings, and frequency. The resulting CPS is consistent with zero because the error is dominated by the foregrounds in the 21cm observation. Additionally, the variance of the signal indicates the presence of unexpected systematics error at small scales. Furthermore, we reduce the error by one order of magnitude with application of a foreground removal using a polynomial fitting method. Based on the results, we find that the detection of the 21cm-CMB CPS with the MWA Phase I requires more than 99.95% of the foreground signal removed, 2000 hours of deep observation and 50% of the sky fraction coverage., Comment: 15 pages, 16 figures, accepted to MNRAS

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., Comment: Nine pages, four figures. Accepted for publication in ApJ

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., Comment: 26 pages, 9 figures, 15 tables

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment., Comment: 18 pages, 11 figures, accepted for publication in PASA

We detail new techniques for analysing ionospheric activity, using Epoch of Reionisation (EoR) datasets obtained with the Murchison Widefield Array (MWA), calibrated by the `Real-Time System' (RTS). Using the high spatial- and temporal-resolution information of the ionosphere provided by the RTS calibration solutions over 19 nights of observing, we find four distinct types of ionospheric activity, and have developed a metric to provide an `at a glance' value for data quality under differing ionospheric conditions. For each ionospheric type, we analyse variations of this metric as we reduce the number of pierce points, revealing that a modest number of pierce points is required to identify the intensity of ionospheric activity; it is possible to calibrate in real-time, providing continuous information of the phase screen. We also analyse temporal correlations, determine diffractive scales, examine the relative fractions of time occupied by various types of ionospheric activity, and detail a method to reconstruct the total electron content responsible for the ionospheric data we observe. These techniques have been developed to be instrument agnostic, useful for application on LOFAR and SKA-Low., Comment: Accepted by MNRAS on 14 July 2017

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., Comment: 9 figs, 2 tables. Accepted for publication in PASA

The current generation of experiments aiming to detect the neutral hydrogen signal from the Epoch of Reionisation (EoR) is likely to be limited by systematic effects associated with removing foreground sources from target fields. In this paper we develop a model for the compact foreground sources in one of the target fields of the MWA's EoR key science experiment: the `EoR1' field. The model is based on both the MWA's GLEAM survey and GMRT 150 MHz data from the TGSS survey, the latter providing higher angular resolution and better astrometric accuracy for compact sources than is available from the MWA alone. The model contains 5049 sources, some of which have complicated morphology in MWA data, Fornax A being the most complex. The higher resolution data show that 13% of sources that appear point-like to the MWA have complicated morphology such as double and quad structure, with a typical separation of 33~arcsec. We derive an analytic expression for the error introduced into the EoR two-dimensional power spectrum due to peeling close double sources as single point sources and show that for the measured source properties, the error in the power spectrum is confined to high $k_\bot$ modes that do not affect the overall result for the large-scale cosmological signal of interest. The brightest ten mis-modelled sources in the field contribute 90% of the power bias in the data, suggesting that it is most critical to improve the models of the brightest sources. With this hybrid model we reprocess data from the EoR1 field and show a maximum of 8% improved calibration accuracy and a factor of two reduction in residual power in $k$-space from peeling these sources. Implications for future EoR experiments including the SKA are discussed in relation to the improvements obtained., Comment: 16 pages, 8 figures, 3 tables. Accepted for publication in PASA

We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis)., Comment: Accepted for publication in Astrophysical Journal on 06 February 2017

We present new software to cross-match low-frequency radio catalogues: the Positional Update and Matching Algorithm (PUMA). PUMA combines a positional Bayesian probabilistic approach with spectral matching criteria, allowing for confusing sources in the matching process. We go on to create a radio sky model using PUMA based on the Murchison Widefield Array Commissioning Survey, and are able to automatically cross-match 98.5% of sources. Using the characteristics of this sky model, we create simple simulated mock catalogues on which to test PUMA, and find that PUMA can reliably find the correct spectral indices of sources, along with being able to recover ionospheric offsets. Finally, we use this sky model to calibrate and remove foreground sources from simulated interferometric data, generated using OSKAR (the Oxford University visibility generator). We demonstrate that there is a substantial improvement in foreground source removal when using higher frequency and higher resolution source positions, even when correcting positions by an average of 0.3 given a synthesized beam-width of 2.3., Comment: Accepted for publication in Publications of the Astronomical Society of Australia (PASA) - 34 pages - 20 figures - PUMA github repository: https://github.com/JLBLine/PUMA

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., Comment: 17 pages, 10 figures, accepted for publication in The Astrophysical Journal

We present and evaluate several strategies to search for prompt, low-frequency radio emission associated with gravitational wave transients using the Murchison Widefield Array (MWA). As we are able to repoint the MWA on timescales of tens of seconds, we can search for the dispersed radio signal that has been predicted to originate along with or shortly after a neutron star-neutron star merger. We find that given the large, 600 deg^2 instantaneous field-of-view of the MWA we can cover a significant fraction of the predicted gravitational wave error region, although due to the complicated geometry of the latter we only cover >50% of the error region for approximately 5% of events, and roughly 15% of events will be located <10 deg from the MWA pointing center such that they will be covered in the radio images. For optimal conditions our limiting flux density for a 10-s long transient would be 0.1 Jy, increasing to about 1 Jy for a wider range of events. This corresponds to luminosity limits of 1e38-1e39 erg/s based on expectations for the distances of the gravitational wave transients, which should be sufficient to detect or significantly constrain a range of models for prompt emission., Comment: 9 pages, 5 figures. Accepted for publication in PASA

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., Comment: 21 pages, 16 figures, accepted to ApJ

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., Comment: 32 pages, 20 figures, accepted to ApJ

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., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal on June 30, 2016

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., Comment: accepted to ApJ

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., Comment: 31 pages, 23 figures, corrected spelling errors in metadata

We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the centre of the image is 4-5 mJy/beam. The MWA counts are in excellent agreement with counts from other instruments and are the most precise ever derived in the flux density range 30-200 mJy due to the sky area covered. Using the deepest available source count data, we find that the MWA image is affected by sidelobe confusion noise at the ~3.5 mJy/beam level, due to incompletely-peeled and out-of-image sources, and classical confusion becomes apparent at ~1.7 mJy/beam. This work highlights that (i) further improvements in ionospheric calibration and deconvolution imaging techniques would be required to probe to the classical confusion limit and (ii) the shape of low-frequency source counts, including any flattening towards lower flux densities, must be determined from deeper ~150 MHz surveys as it cannot be directly inferred from higher frequency data., Comment: 13 pages, 13 figures, accepted for publication in MNRAS

We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field., Comment: MNRAS Accepted, 17 pages, 10 figures

Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the MWA EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by PSF sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source., Comment: Accepted for publication in MNRAS. 15 pages, 16 figures. Catalogue of sources externally available as CSV file

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., Comment: 13 pages, 9 figures, matches version accepted to ApJ

We present the results of a pilot study search for Fast Radio Bursts (FRBs) using the Murchison Widefield Array (MWA) at low frequencies (139 - 170 MHz). We utilised MWA data obtained in a routine imaging mode from observations where the primary target was a field being studied for Epoch of Reionisation detection. We formed images with 2 second time resolution and 1.28~MHz frequency resolution for 10.5 hours of observations, over 400 square degrees of the sky. We de-dispersed the dynamic spectrum in each of 372,100 resolution elements of 2$\times$2 arcmin$^{2}$, between dispersion measures of 170 and 675~pc~cm$^{-3}$. Based on the event rate calculations in Trott, Tingay & Wayth (2013), which assumes a standard candle luminosity of $8\times10^{37}$ Js$^{-1}$, we predict that with this choice of observational parameters, the MWA should detect ($\sim10$,$\sim2$,$\sim0$) FRBs with spectral indices corresponding to ($-$2, $-$1, 0), based on a 7$\sigma$ detection threshold. We find no FRB candidates above this threshold from our search, placing an event rate limit of $<700$ above 700 Jy.ms per day per sky and providing evidence against spectral indices $\alpha<-1.2$ ($S\propto\nu^{\alpha}$). We compare our event rate and spectral index limits with others from the literature. We briefly discuss these limits in light of recent suggestions that supergiant pulses from young neutron stars could explain FRBs. We find that such supergiant pulses would have to have much flatter spectra between 150 and 1400 MHz than have been observed from Crab giant pulses to be consistent with the FRB spectral index limit we derive., Comment: 14 pages, 5 figures

We present observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to >1e6 km. The serendipitous night-time nature of these detections illustrates the new regime that the MWA has opened for IPS studies with sensitive night-time, wide-field, low-frequency observations. This regime complements traditional dedicated strategies for observing IPS and can be utilized in real-time to facilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact (arcsec) structures in low-frequency radio sources despite the 2 arcmin resolution of the array., Comment: 7 pages, 4 figures. ApJ Letters, in press

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., Comment: Accepted for publication in ApJ Letters, 6 pages, 2 figures