Search

The locations of Ly-$\alpha$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for the spatial distribution of the gas, thereby providing prior information unavailable to power spectrum measurements. We use a test set of 12 hours of observations from the Murchison Widefield Array (MWA) in extended array configuration, to constrain the neutral hydrogen signature of 58 LAEs, detected with the Subaru Hypersuprime Cam in the \textit{Silverrush} survey, centred on $z$=6.58. We assume that detectable emitters reside in the centre of ionised HII bubbles during the end of reionization, and predict the redshifted neutral hydrogen signal corresponding to the remaining neutral regions using a set of different ionised bubble radii. A prewhitening matched filter detector is introduced to assess detectability. We demonstrate the ability to detect, or place limits upon, the amplitude of brightness temperature fluctuations, and the characteristic HII bubble size. With our limited data, we constrain the brightness temperature of neutral hydrogen to $\Delta{\rm T}_B, Comment: 9 pages, 6 figures, accepted for publication in MNRAS

The far-field pattern of a geometrically large and complex antenna used in low-frequency radio astronomy is computationally expensive to simulate on electromagnetic simulators, such as FEKO. For example, one station of the Square Kilometer Array, which consists of 256 log-periodic antenna elements, will take years to simulate using the full CAD model for the full operational frequency band. This paper focuses on reducing the simulation time for a single antenna element by simplifying the simulation model, thus decreasing the number of unknowns that have to be solved in a simulation. An iterative process for optimizing the simplification of such an element is described, while keeping the reflection coefficient within 1 dB absolute mean deviation of the measured data. After four iterations, the amount of unknowns to be solved, which includes the number of triangles and segments, was reduced from 29,307 to 11,991. This decreased the computation time by 86.5%, making array simulations feasible. Using the techniques described in the paper, other antenna constructions can benefit from it and be simulated more efficiently.

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)., 18 pages, 13 figures

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., This article has been accepted for publication in PASA. This article has 16 pages, 7 figures and 1 table

We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year., 19 pages, 11 figures

In radio astronomy, holography is a commonly used technique to create an image of the electric field distribution in the aperture of a dish antenna. The image is used to detect imperfections in the reflector surface. Similarly, holography can be applied to phased array telescopes, in order to measure the complex gains of the receive paths of individual antennas. In this paper, a holographic technique is suggested to calibrate the digital beamformer of a phased array telescope. The effectiveness of the technique was demonstrated by applying it on data from the Engineering Development Array 2, one of the prototype stations of the low frequency component of the Square Kilometre Array. The calibration method is very quick and requires few resources. In contrast to holography for dish antennas, it works without a reference antenna. We demonstrate the utility of this technique for initial station commissioning and verification as well as for routine station calibration., 14 pages, 10 figures. Accepted for publication in Radio Science

The sensitivity of a single station is a crucial performance characteristic of the low-frequency component of the Square Kilometre Array (SKA-Low). In this contribution, the sensitivity of a prototype station of the SKA-Low, the Aperture Array Verification System Version 2 (AAVS2), measured at frequencies 70.3, 110.2, 159.4, 229.7 and 320.3MHz is presented and compared with electromagnetic simulations. The measured sensitivity data are in good agreement with the predictions of electromagnetic simulations, confirming good understanding of the instrument and its accurate representation in the simulations.

We present the South Galactic Pole (SGP) data release from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. These data combine both years of GLEAM observations at 72-231 MHz conducted with the Murchison Widefield Array (MWA) and cover an area of 5,113 $\mathrm{deg}^2$ centred on the SGP at 20$^\mathrm{h}$40$^\mathrm{m}$ < RA < 05$^\mathrm{h}$04$^\mathrm{m}$ and $-48��$ < Dec < $-2��$. At 216 MHz, the typical rms noise is $\approx 5$ mJy/beam and the angular resolution $\approx 2$ arcmin. The source catalogue contains a total of 108,851 components above $5��$, of which 77 per cent have measured spectral indices between 72 and 231 MHz. Improvements to the data reduction in this release include the use of the GLEAM Extragalactic catalogue as a sky model to calibrate the data, a more efficient and automated algorithm to deconvolve the snapshot images, and a more accurate primary beam model to correct the flux scale. This data release enables more sensitive large-scale studies of extragalactic source populations as well as spectral variability studies on a one-year timescale., 13 pages, 7 figures, accepted for publication in PASA. Added link to images and source catalogue from GLEAM SGP data release in section 5. Corrected some mistakes in column numbers in section 4.5 describing catalogue. The catalogue is now also available on VizieR

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

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

The GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) is a radio continuum survey at 76-227 MHz of the entire southern sky (Declination $ 55$ mJy across an area of $\approx 16,700~\mathrm{deg}^{2}$. From the optical spectra, we identify the dominant physical process responsible for the radio emission from each galaxy: 73 per cent are fuelled by an active galactic nucleus (AGN) and 27 per cent by star formation. We present the local radio luminosity function for AGN and star-forming galaxies at 200 MHz and characterise the typical radio spectra of these two populations between 76 MHz and $\sim 1$ GHz. For the AGN, the median spectral index between 200 MHz and $\sim 1$ GHz, $\alpha_{\mathrm{high}}$, is $-0.600 \pm 0.010$ (where $S \propto \nu^{\alpha}$) and the median spectral index within the GLEAM band, $\alpha_{\mathrm{low}}$, is $-0.704 \pm 0.011$. For the star-forming galaxies, the median value of $\alpha_{\mathrm{high}}$ is $-0.650 \pm 0.010$ and the median value of $\alpha_{\mathrm{low}}$ is $-0.596 \pm 0.015$. Among the AGN population, flat-spectrum sources are more common at lower radio luminosity, suggesting the existence of a significant population of weak radio AGN that remain core-dominated even at low frequencies. However, around 4 per cent of local radio AGN have ultra-steep radio spectra at low frequencies ($\alpha_{\mathrm{low}} < -1.2$). These ultra-steep-spectrum sources span a wide range in radio luminosity, and further work is needed to clarify their nature., Comment: 28 pages, 18 figures, accepted for publication in PASA

In the context of the Radio over Fiber downlinks belonging to the Aperture Array Verification System (AAVS) of the low frequency Square Kilometre Array (SKA-LOW) radio telescope, the combination of modulation of the optical carrier by Radio Frequency Interfering Tones (RFIs) and fiber-induced Rayleigh Backscattering (RBS) determines the presence of undesired spurious frequency terms. A successful countermeasure to the phenomenon has been proposed, based on the additional modulation with a low frequency tone (called dithering tone) which mitigates such distortions exploiting the laser chirp effect. In the present work, through a rigorous theoretical and experimental study, the optimal design parameters of such solution are put into evidence, showing that an appropriate choice of the dithering tone frequency allows to reduce to acceptable levels the RBS induced nonlinearities, while minimizing the power consumption of the corresponding introduced devices.

Here we present stringent low-frequency 185MHz limits on coherent radio emission associated with a short gamma-ray burst (SGRB). Our observations of the short GRB 180805A were taken with the upgraded Murchison Widefield Array (MWA) rapid-response system, which triggered within 20s of receiving the transient alert from Swift, corresponding to 83.7s post-burst. The SGRB was observed for 30m, resulting in a 3sigma persistent flux density upper-limit of 40.2mJy/beam. Transient searches were conducted at the Swift position of this GRB on 0.5s, 5s, 30s, and 2m timescales, resulting in 3sigma limits of 570-1830, 270-630, 200-420, and 100-200mJy/beam, respectively. We also performed a dedispersion search for prompt signals at the position of the SGRB with a temporal and spectral resolution of 0.5s and 1.28MHz, resulting in a 6sigma fluence upper-limit range from 570Jyms at DM=3000pc/cm^3 (z~2.5) to 1750Jyms at DM=200pc/cm^3 (z~0.1). We compare the fluence prompt emission limit and the persistent upper-limit to SGRB coherent emission models assuming the merger resulted in a stable magnetar. Our observations were not sensitive enough to detect prompt emission associated with the alignment of magnetic fields of a binary neutron star just prior to the merger, from the interaction between the relativistic jet and the interstellar medium or persistent pulsar-like emission from the spin-down of the magnetar. However, in the case of a more powerful SGRB (a gamma-ray fluence an order of magnitude higher than GRB 180805A and/or a brighter X-ray counterpart), our MWA observations may be sensitive enough to detect coherent radio emission from the jet-ISM interaction and/or the magnetar remnant. Finally, we demonstrate that of all current low-frequency radio telescopes, only the MWA has the sensitivity and response times capable of probing prompt emission models associated with the initial SGRB merger event., Comment: 20 pages, 11 figures, accepted for publication in PASA

We report the discovery of the first new pulsar with the Murchison Widefield Array (MWA), PSR J0036$-$1033, a long-period (0.9 s) nonrecycled pulsar with a dispersion measure (DM) of 23.1 ${\rm pc\,cm^{-3}}$. It was found after processing only a small fraction ($\sim$1%) of data from an ongoing all-sky pulsar survey. Follow-up observations have been made with the MWA, the upgraded Giant Metrewave Radio Telescope (uGMRT), and the Parkes 64 m telescopes, spanning a frequency range from $\sim$150 MHz to 4 GHz. The pulsar is faint, with an estimated flux density ($S$) of $\sim$1 mJy at 400 MHz and a spectrum $S(\nu)\,\propto\,\nu^{-2.0 \pm 0.2}$, where $\nu$ is frequency. The DM-derived distance implies that it is also a low-luminosity source ($\sim$ 0.1 ${\rm mJy\,kpc^2}$ at 1400 MHz). The analysis of archival MWA observations reveals that the pulsar's mean flux density varies by up to a factor of $\sim$5-6 on timescales of several weeks to months. By combining MWA and uGMRT data, the pulsar position was determined to arcsecond precision. We also report on polarization properties detected in the MWA and Parkes bands. The pulsar's nondetection in previous pulsar and continuum imaging surveys, the observed high variability, and its detection in a small fraction of the survey data searched to date, all hint at a larger population of pulsars that await discovery in the southern hemisphere, with the MWA and the future low-frequency Square Kilometre Array., Comment: 12 pages, 6 figures, 2 tables, Accepted for publication in The ApJ Letters

We present the first southern-hemisphere all-sky imager and radio-transient monitoring system implemented on two prototype stations of the low-frequency component of the Square Kilometre Array. Since its deployment the system has been used for real-time monitoring of the recorded commissioning data. Additionally, a transient searching algorithm has been executed on the resulting all-sky images. It uses a difference imaging technique, and has enabled identification of a wide variety of transient classes, ranging from human-made radio-frequency interference to genuine astrophysical events. Observations at the frequency 159.4 MHz and higher in a single coarse channel (0.926 MHz) were made with 2s time resolution, and multiple nights were analysed. Despite having modest sensitivity (~few Jy/beam), using a single coarse channel and 2-s imaging, the system detected bright transients from PSR B0950+08, proving that it can be used to detect bright transients of an astrophysical origin. The unusual, extreme activity of the pulsar PSR B0950+08 (up to ~155 Jy/beam) was initially detected in a "blind" search in the 2020-04-10/11 data and later assigned to this specific pulsar. The limitations of our data, however, prevent use from making firm conclusions of the effect being due to a combination of refractive and diffractive scintillation or intrinsic emission mechanisms. The system can routinely collect data over many days without interruptions; the large amount of recorded data at 159.4 and 229.7 MHz allowed us to determine a preliminary transient surface density upper limit of $1.32 \times 10^{-9} \text{deg}^{-2}$ for a timescale and limiting flux density of 2s and 42 Jy, respectively. We plan to extend the observing bandwidth to tens of MHz and improve time resolution to tens of milliseconds in order to increase the sensitivity and enable detections of Fast Radio Bursts below 300 MHz., Comment: 22 pages, 11 figures, 3 tables, Accepted for publication in PASA

Radio interferometer arrays with non-homogeneous element patterns are more difficult to calibrate compared to the more common homogeneous array. In particular, the non-homogeneity of the patterns has significant implications on the computational tractability of evaluating the calibration solutions. We apply the A-stacking technique to this problem and explore the trade-off to be made between the calibration accuracy and computational complexity. Through simulations, we show that this technique can be favourably applied in the context of an SKA-Low station. We show that the minimum accuracy requirements can be met at a significantly reduced computational cost, and this cost can be reduced even further if the station calibration time-scale is relaxed from 10 min to several hours. We demonstrate the impact antenna designs with differing levels of non-homogeneity have on the overall computational complexity in addition to some cases where calibration performs poorly.

The design and development process for the Square Kilometre Array (SKA) radio telescope’s Low Frequency Aperture Array component was progressed during the SKA pre-construction phase by an international consortium, with the goal of meeting requirements for a critical design review. As part of the development process a full-sized prototype SKA Low ‘station’ was deployed – the Aperture Array Verification System 1 (AAVS1). We provide a system overview and describe the commissioning results of AAVS1, which is a low frequency radio telescope with 256 dual-polarisation log-periodic dipole antennas working as a phased array. A detailed system description is provided, including an in-depth overview of relevant sub-systems, ranging from hardware, firmware, software, calibration, and control sub-systems. Early commissioning results cover initial bootstrapping, array calibration, stability testing, beam-forming, and on-sky sensitivity validation. Lessons learned are presented, along with future developments.

The Square Kilometre Array telescope at low-frequency (SKA-Low) will be a phased array telescope supporting a wide range of science cases covering the frequency band 50 - 350 MHz, while at the same time asking for high sensitivity and excellent characteristics. These extremely challenging requirements resulted in a design using 512 groups of 256 log periodic dual polarized antennas each (where each group is called "station"), for a total of 131072 antennas. The 512 stations are randomly distributed mostly within a dense area around the centre of the SKA-Low, and then in 3 arms having 16 station clusters each. In preparation for the SKA Phase 1 (SKA1) System Critical Design Review (CDR), prototype stations were deployed at the Murchison Radio-astronomy Observatory (MRO) site (Western Australia) near the Murchison Widefield Array (MWA) radio telescope. The project involved multiple parties in an International collaboration building and testing different prototypes of the SKA1-Low station near the actual site. This resulted in both organisational and logistic challenges typical of a deployment of the actual telescope. The test set-up involved a phased build-up of the complex station of log-periodic antennas, by starting from the deployment of 48 antennas and related station signal processing (called AAVS1.5, where AAVS stands for Aperture Array Verification System), followed by expansion to a full station (AAVS2.0). As reference a station with dipole antennas EDA2 (EDA: Engineering Development Array) was deployed. This test set-up was used for an extensive test and evaluation programme. All test antenna configurations were simulated in detail by electromagnetic (EM) models, and the prediction of the models was further verified by appropriate tests with a drone-based radio frequency source. Astronomical observations on Sun and galaxy transit were performed with calibrated stations of both EDA2, AAVS1.5 and AAVS2.0. All 3 configurations were calibrated. EM modelling and calibration results for the full station AAVS2.0 and EM verification for the AAVS1.5 station are presented. The comparisons between the behaviour of the log-periodic antennas and the dipoles have advanced our understanding the calibration quality and the technological maturity of the future SKA1-Low array.

We present the first survey of radio frequency interference (RFI) at the future site of the low frequency Square Kilometre Array (SKA), the Murchison Radio-astronomy Observatory (MRO), that both temporally and spatially resolves the RFI. The survey is conducted in a 1 MHz frequency range within the FM band, designed to encompass the closest and strongest FM transmitters to the MRO (located in Geraldton, approximately 300 km distant). Conducted over approximately three days using the second iteration of the Engineering Development Array in an all-sky imaging mode, we find a range of RFI signals. We are able to categorise the signals into: those received directly from the transmitters, from their horizon locations; reflections from aircraft (occupying approximately 13% of the observation duration); reflections from objects in Earth orbit; and reflections from meteor ionisation trails. In total we analyse 33,994 images at 7.92 s time resolution in both polarisations with angular resolution of approximately 3.5 deg., detecting approximately forty thousand RFI events. This detailed breakdown of RFI in the MRO environment will enable future detailed analyses of the likely impacts of RFI on key science at low radio frequencies with the SKA., Accepted for publication in PASA. 15 pages, 16 figures

The computational electromagnetic modelling of a large radio telescope prototype array for the Square Kilometer Array is described. The numerical models, using the Method of Moments, are characterised by a very large number of unknowns, requiring the use of fast solution methods and high performance computing platforms. Good agreement has been obtained between results obtained on two different commercial codes. Results for both embedded element patterns and the station beam are shown. The use of the computed embedded element patterns for array calibration is briefly addressed.

The entire southern sky (Declination, $\delta $ 4 Jy) in the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg). We refer to these 1,863 sources as the GLEAM 4-Jy (G4Jy) Sample, and use radio images (of $\leq$ 45-arcsec resolution), and multi-wavelength information, to assess their morphology and identify the galaxy that is hosting the radio emission (where appropriate). Details of how to access all of the overlays used for this work are available at https://github.com/svw26/G4Jy. Alongside this we conduct further checks against the literature, which we document in this paper for individual sources. Whilst the vast majority of the G4Jy Sample are active galactic nuclei with powerful radio-jets, we highlight that it also contains a nebula, two nearby, star-forming galaxies, a cluster relic, and a cluster halo. There are also three extended sources for which we are unable to infer the mechanism that gives rise to the low-frequency emission. In the G4Jy catalogue we provide mid-infrared identifications for 86% of the sources, and flag the remainder as: having an uncertain identification (129 sources), having a faint/uncharacterised mid-infrared host (126 sources), or it being inappropriate to specify a host (2 sources). For the subset of 129 sources, there is ambiguity concerning candidate host-galaxies, and this includes four sources (B0424$-$728, B0703$-$451, 3C 198, and 3C 403.1) where we question the existing identification., Comment: 37 pages (24 MB in size), 23 figures, 3 tables, accepted for publication in PASA, and now updated to match the final proof. Full-resolution images will be used for the published version, available through the journal

A number of natural phenomena occurring at the Western Australian site of the Square Kilometre Array (SKA), such as lightning or whirlwinds, can cause damage to electronic parts of antenna array elements. We consider an SKA precursor, the Murchison Widefield Array (MWA), in order to evaluate the consequences of faulty low-noise amplifiers in an array. Using network analysis methods, we predict exactly the field patterns of each array element under changing loading conditions. The values of load impedances used for MWA field calculations were based on measured data. It was observed that the tile pattern of MWA is robust to occasional low-noise amplifier damage. The ability to predict the embedded element patterns without the need to repeat numerical simulations is especially useful for such large-scale projects as the SKA.

In this paper we propose the use of a parallel-plate waveguide simulator as a useful design verification step of very large phased-array systems. We base the derivation of the theoretical concept upon the wideband capacitively connected-dipole array. It is shown to be key to correctly terminate the cavity to the free-space boundary with the use of electromagnetic absorbers to minimize reflections.

Optical links based on the Radio over Fiber (RoF) Technology are, and are going to be, utilized within the antenna downlinks of important Radioastronomic facilities. In this context, the presence of spurious frequency terms has been observed despite the relatively low power levels of the transmitted of RF signals. The present work gives the theoretical explanation of the onset of these nonlinearities, and rigorously demonstrates an effective countermeasure to the undesired phenomenon.

Diffuse, non-thermal emission in galaxy clusters is increasingly being detected in low-frequency radio surveys and images. We present a new diffuse, steep-spectrum, non-thermal radio source within the cluster Abell 1127 found in survey data from the Murchison Widefield Array (MWA). We perform follow-up observations with the 'extended' configuration MWA Phase II with improved resolution to better resolve the source and measure its low-frequency spectral properties. We use archival Very Large Array S-band data to remove the discrete source contribution from the MWA data, and from a power law model fit we find a spectral index of $-1.83\pm0.29$ broadly consistent with relic-type sources. The source is revealed by the Giant Metrewave Radio Telescope (GMRT) at 150 MHz to have an elongated morphology, with a projected linear size of 850 kpc as measured in the MWA data. Using Chandra observations we derive morphological estimators and confirm quantitatively that the cluster is in a disturbed dynamical state, consistent with the majority of phoenices and relics being hosted by merging clusters. We discuss the implications of relying on morphology and low-resolution imaging alone for the classification of such sources and highlight the usefulness of the MHz to GHz radio spectrum in classifying these types of emission. Finally, we discuss the benefits and limitations of using the MWA Phase II in conjunction with other instruments for detailed studies of diffuse, steep-spectrum, non-thermal radio emission within galaxy clusters., Comment: 17 pages, 16 figures, accepted for publication in Publications of the Astronomical Society of Australia (PASA)

The Murchison Widefield Array (MWA) has observed the entire southern sky (Declination, $\delta $ 10 deg) to define the GLEAM 4-Jy (G4Jy) Sample. This is a complete sample of the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy), the majority of which are active galactic nuclei with powerful radio-jets. Crucially, low-frequency observations allow the selection of such sources in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). We then use higher-resolution radio images, and information at other wavelengths, to morphologically classify the brightest components in GLEAM. We also conduct cross-checks against the literature, and perform internal matching, in order to improve sample completeness (which is estimated to be $>$ 95.5%). This results in a catalogue of 1,863 sources, making the G4Jy Sample over 10 times larger than that of the revised Third Cambridge Catalogue of Radio Sources (3CRR; $S_{\mathrm{178MHz}} >$ 10.9 Jy). Of these G4Jy sources, 78 are resolved by the MWA (Phase-I) synthesised beam ($\sim$2 arcmin at 200 MHz), and we label 67% of the sample as 'single', 26% as 'double', 4% as 'triple', and 3% as having 'complex' morphology at $\sim$1 GHz (45-arcsec resolution). Alongside this, our value-added catalogue provides mid-infrared source associations (subject to 6-arcsec resolution at 3.4 micron) for the radio emission, as identified through visual inspection and thorough checks against the literature. As such, the G4Jy Sample can be used as a reliable training set for cross-identification via machine-learning algorithms. [Abstract abridged for arXiv submission.], Comment: 57 pages (30 MB in size), 23 figures, 16 tables, accepted for publication in PASA, and now updated to match the final proof. Full-resolution images will be used for the published version, available through the journal. To keep up-to-date regarding the G4Jy Sample, bookmark https://github.com/svw26/G4Jy

Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in \cite{Li_2018} and \cite{Wenyang_2019} studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the power spectrum from tandem calibration are significant. To understand this result, we analyze both the calibration solutions themselves and the effects on the power spectrum over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model-incompleteness error., Comment: 20 pages, 11 figures, accepted for publication in PASA

The search for emission from weakly interacting massive particle (WIMP) dark matter annihilation and decay has become a multi-pronged area of research not only targeting a diverse selection of astrophysical objects, but also taking advantage of the entire electromagnetic spectrum. The decay of WIMP particles into standard model particles has been suggested as a possible channel for synchrotron emission to be detected at low radio frequencies. Here, we present the stacking analysis of a sample of 33 dwarf spheroidal (dSph) galaxies with low-frequency (72 - 231 MHz) radio images from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. We produce radial surface brightness profiles of images centred upon each dSph galaxy with background radio sources masked. We remove ten fields from the stacking due to contamination from either poorly subtracted, bright radio sources or strong background gradients across the field. The remaining 23 dSph galaxies are stacked in an attempt to obtain a statistical detection of any WIMP-induced synchrotron emission in these systems. We find that the stacked radial brightness profile does not exhibit a statistically significant detection above the 95% confidence level of $\sim$1.5 mJy beam$^{-1}$. This novel technique shows the potential of using low-frequency radio images to constrain fundamental properties of particle dark matter., Comment: 11 pages, 7 Figures, Accepted for publication in MNRAS

We present a calibration component for the Murchison Widefield Array All-Sky Virtual Observatory (MWA ASVO) utilising a newly developed PostgreSQL database of calibration solutions. Since its inauguration in 2013, the MWA has recorded over thirty-four petabytes of data archived at the Pawsey Supercomputing Centre. According to the MWA Data Access policy, data become publicly available eighteen months after collection. Therefore, most of the archival data are now available to the public. Access to public data was provided in 2017 via the MWA ASVO interface, which allowed researchers worldwide to download MWA uncalibrated data in standard radio astronomy data formats (CASA measurement sets or UV FITS files). The addition of the MWA ASVO calibration feature opens a new, powerful avenue for researchers without a detailed knowledge of the MWA telescope and data processing to download calibrated visibility data and create images using standard radio-astronomy software packages. In order to populate the database with calibration solutions from the last six years we developed fully automated pipelines. A near-real-time pipeline has been used to process new calibration observations as soon as they are collected and upload calibration solutions to the database, which enables monitoring of the interferometric performance of the telescope. Based on this database we present an analysis of the stability of the MWA calibration solutions over long time intervals., 12 pages, 9 figures, Accepted for publication in PASA

We derive the synchrotron distribution in the Milky Way disk from HII region absorption observations over -40{\deg} < l < 40{\deg} at six frequencies of 76.2, 83.8, 91.5, 99.2, 106.9, and 114.6 MHz with the GaLactic and Extragalactic All-sky Murchison widefield array survey (GLEAM). We develop a new method of emissivity calculation by taking advantage of the Haslam et al., (1981) map and known spectral indices, which enable us to simultaneously derive the emissivity and the optical depth of HII regions at each frequency. We show our derived synchrotron emissivities based on 152 absorption features of HII regions using both the method previously adopted in the literature and our improved method. We derive the synchrotron emissivity from HII regions to the Galactic edge along the line of sight and, for the first time, derive the emissivity from HII regions to the Sun. These results provide direct information on the distribution of the Galactic magnetic field and cosmic-ray electrons for future modelling., Comment: Accepted for publication in MNRAS, 20 pages, 9 figures, 2 tables

We present new radio and optical images of the nearest radio galaxy Centaurus A and its host galaxy NGC 5128. We focus our investigation on the northern transition region, where energy is transported from the ~5 kpc (~5 arcmin) scales of the Northern Inner Lobe (NIL) to the ~30 kpc (~30 arcmin) scales of the Northern Middle Lobe (NML). Our Murchison Widefield Array observations at 154 MHz and our Parkes radio telescope observations at 2.3 GHz show diffuse radio emission connecting the NIL to the NML, in agreement with previous Australia Telescope Compact Array observations at 1.4 GHz. Comparison of these radio data with our widefield optical emission line images show the relationship between the NML radio emission and the ionised filaments that extend north from the NIL, and reveal a new ionised filament to the east, possibly associated with a galactic wind. Our deep optical images show clear evidence for a bipolar outflow from the central galaxy extending to intermediate scales, despite the non-detection of a southern radio counterpart to the NML. Thus, our observational overview of Centaurus A reveals a number of features proposed to be associated with AGN feedback mechanisms, often cited as likely to have significant effects in galaxy evolution models. As one of the closest galaxies to us, Centaurus A therefore provides a unique laboratory to examine feedback mechanisms in detail., 19 pages, 9 figures, accepted by MNRAS

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., Accepted by MNRAS on 14 July 2017

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., 14 pages, 8 figures. Accepted for publication in MNRAS December 2019. Erratum submitted on 20 February 2020

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 Engineering Development Array (EDA) is a single test station for Square Kilometre Array (SKA) precursor technology. We have used the EDA to detect low-frequency radio recombination lines (RRLs) from the Galactic Centre region. Low-frequency RRLs are an area of interest for future low-frequency SKA work as these lines provide important information on the physical properties of the cold neutral medium. In this project we investigate the EDA, its bandpass and the radio frequency interference environment for low-frequency spectroscopy. We present line spectra from 30 to 325 MHz for the Galactic Centre region. The decrease in sensitivity for the EDA at the low end of the receiver prevents carbon and hydrogen RRLs to be detected below 40 and 60 MHz respectively. RFI strongly affects frequencies in the range 276-292, 234-270, 131-138, 95-102 and below 33 MHz. Cn$\alpha$ RRLs were detected in absorption for quantum levels n = 378 to 550 (39-121 MHz) and in emission for n = 272 to 306 (228-325 MHz). Cn$\beta$ lines were detected in absorption for n = 387 to 696 (39-225 MHz). Hn$\alpha$ RRLs were detected in emission for n = 272 to 480 (59-325 MHz). Hn$\beta$ lines were detected for n = 387 to 453 (141-225 MHz). The stacked Hn$\alpha$ detection at 63 MHz is the lowest frequency detection made for hydrogen RRLs and shows that a cold (partially) ionized medium exists along the line of sight to the Galactic Centre region. The size and velocity of this cold H$^{+}$ gas indicates that it is likely associated with the nearby Riegel-Crutcher cloud., Comment: 18 pages, 6 figures and 5 tables

We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights., 19 pages, 18 figures, accepted for publication in MNRAS

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$., Accepted for publication in A&A

We have detected 27 new supernova remnants (SNRs) using a new data release of the GLEAM survey from the Murchison Widefield Array telescope, including the lowest surface brightness SNR ever detected, G 0.1 – 9.7. Our method uses spectral fitting to the radio continuum to derive spectral indices for 26/27 candidates, and our low-frequency observations probe a steeper spectrum population than previously discovered. None of the candidates have coincidentWISEmid-IR emission, further showing that the emission is non-thermal. Using pulsar associations we derive physical properties for six candidate SNRs, finding G 0.1 – 9.7 may be younger than 10 kyr. Sixty per cent of the candidates subtend areas larger than 0.2 deg2on the sky, compared to < 25% of previously detected SNRs. We also make the first detection of two SNRs in the Galactic longitude range 220°–240°.

The numerical modelling of the Square Kilometre Array (SKA) Aperture Array Verification System Version 1.5 is discussed. The role of this SKA prototype is placed within the perspective of the overall SKA project. The dual-polarized log-periodic SKALA4.1 antenna elements comprising a station are briefly described, along with some considerations regarding the station. Computational simulation aspects are considered, and preliminary results shown. An appraisal of the implications for station-level calibration concludes the paper.

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