Source author record

S. J. Tingay

S. J. Tingay appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

108works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

108 published item(s)

preprint2022arXiv

Search for non-standard neutrino interactions with 10 years of ANTARES data

Non-standard interactions of neutrinos arising in many theories beyond the Standard Model can significantly alter matter effects in atmospheric neutrino propagation through the Earth. In this paper, a search for deviations from the prediction of the standard 3-flavour atmospheric neutrino oscillations using the data taken by the ANTARES neutrino telescope is presented. Ten years of atmospheric neutrino data collected from 2007 to 2016, with reconstructed energies in the range from $\sim$16 GeV to $100$ GeV, have been analysed. A log-likelihood ratio test of the dimensionless coefficients $\varepsilon_{μτ}$ and $\varepsilon_{ττ} - \varepsilon_{μμ}$ does not provide clear evidence of deviations from standard interactions. For normal neutrino mass ordering, the combined fit of both coefficients yields a value 1.7$σ$ away from the null result. However, the 68% and 95% confidence level intervals for $\varepsilon_{μτ}$ and $\varepsilon_{ττ} - \varepsilon_{μμ}$, respectively, contain the null value. Best fit values, one standard deviation errors and bounds at the 90% confidence level for these coefficients are given for both normal and inverted mass orderings. The constraint on $\varepsilon_{μτ}$ is among the most stringent to date and it further restrains the strength of possible non-standard interactions in the $μ- τ$ sector.

preprint2022arXiv

Search for secluded dark matter towards the Galactic Centre with the ANTARES neutrino telescope

Searches for dark matter (DM) have not provided any solid evidence for the existence of weakly interacting massive particles in the GeV-TeV mass range. Coincidentally, the scale of new physics is being pushed by collider searches well beyond the TeV domain. This situation strongly motivates the exploration of DM masses much larger than a TeV. Secluded scenarios contain a natural way around the unitarity bound on the DM mass, via the early matter domination induced by the mediator of its interactions with the Standard Model. High-energy neutrinos constitute one of the very few direct accesses to energy scales above a few TeV. An indirect search for secluded DM signals has been performed with the ANTARES neutrino telescope using data from 2007 to 2015. Upper limits on the DM annihilation cross section for DM masses up to 6 PeV are presented and discussed.

preprint2022arXiv

Search for solar atmospheric neutrinos with the ANTARES neutrino telescope

Solar Atmospheric Neutrinos (SA$ν$s) are produced by the interaction of cosmic rays with the solar medium. The detection of SA$ν$s would provide useful information on the composition of primary cosmic rays as well as the solar density. These neutrinos represent an irreducible source of background for indirect searches for dark matter towards the Sun and the measurement of their flux would allow for a better assessment of the uncertainties related to these searches. In this paper we report on the analysis performed, based on an unbinned likelihood maximisation, to search for SA$ν$s with the ANTARES neutrino telescope. After analysing the data collected over 11 years, no evidence for a solar atmospheric neutrino signal has been found. An upper limit at 90\% confidence level on the flux of solar atmospheric neutrinos has been obtained, equal to 7$\times$$10^{-11}$ [TeV$^{-1}$cm$^{-2}$s$^{-1}$] at E$_ν=$ 1 TeV for the reference cosmic ray model assumed.

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2022arXiv

What is the SKA-Low Sensitivity for Your Favourite Radio Source ?

The SKA will be the largest radio astronomy observatory ever built, providing unprecedented sensitivity over a very broad frequency (50 MHz to 15.3 GHz). The SKA-Low (50 - 350 MHz), will be built at the MRO in Western Australia. It will consist of 512 stations each composed of 256 dual-polarised antennas, and the sensitivity of an individual station is pivotal to the performance of the entire SKA-Low telescope. The answer to the question in the title is, it depends. The sensitivity of a low frequency array, such as an SKA-Low station, depends strongly on the pointing direction of the digitally formed station beam and the local sidereal time (LST), and is different for the two orthogonal polarisations of the antennas. The accurate prediction of the SKA-Low sensitivity in an arbitrary direction in the sky is crucial for future observation planning. We present here a sensitivity calculator for the SKA-Low radio telescope, using a database of pre-computed sensitivity values for two realisations of an SKA-Low station architecture. One realisation uses the log-periodic antennas selected for SKA-Low. The second uses a known benchmark, in the form of the bowtie dipoles of the MWA. Data collected by both stations (deployed at the MRO in 2019) were used to measure their sensitivity at selected frequencies and over at least 24 h intervals, and were compared to the predictions described in this paper. The sensitivity values stored in the SQLite database were pre-computed for the X, Y and Stokes I polarisations in 10 MHz frequency steps, 0.5 hour LST intervals, and 5 degree resolution in pointing directions. The database allows users to estimate the sensitivity of SKA-Low for their favourite object using interactive web-based or command line interface, which can also calculate the sensitivity for arbitrary pointing directions, frequencies, and times without interpolations.

preprint2021arXiv

Early-time Searches for Coherent Radio Emission from Short GRBs with the Murchison Widefield Array

Here we present a low frequency (170-200MHz) search for coherent radio emission associated with nine short GRBs detected by the Swift and/or Fermi satellites using the Murchison Widefield Array (MWA) rapid-response observing mode. The MWA began observing these events within 30 to 60s of their high-energy detection, enabling us to capture any dispersion delayed signals emitted by short GRBs for a typical range of redshifts. We conducted transient searches at the GRB positions on timescales of 5s, 30s and 2min, resulting in the most constraining flux density limits on any associated transient of 0.42, 0.29, and 0.084Jy, respectively. We also searched for dispersed signals at a temporal and spectral resolution of 0.5s and 1.28MHz but none were detected. However, the fluence limit of 80-100Jy ms derived for GRB 190627A is the most stringent to date for a short GRB. We compared the fluence and persistent emission limits to short GRB coherent emission models, placing constraints on key parameters including the radio emission efficiency of the nearly merged neutron stars ($\lesssim10^{-4}$), the fraction of magnetic energy in the GRB jet ($\lesssim2\times10^{-4}$), and the radio emission efficiency of the magnetar remnant ($\lesssim10^{-3}$). Comparing the limits derived for our full GRB sample to the same emission models, we demonstrate that our 30-min flux density limits were sensitive enough to theoretically detect the persistent radio emission from magnetar remnants up to a redshift of $z\sim0.6$. Our non-detection of this emission could imply that some GRBs in the sample were not genuinely short or did not result from a binary neutron star merger, the GRBs were at high redshifts, these mergers formed atypical magnetars, the radiation beams of the magnetar remnants were pointing away from Earth, or the majority did not form magnetars but rather collapse directly into black holes.

preprint2021arXiv

Remnant Radio Galaxies Discovered in a Multi-frequency Survey

The remnant phase of a radio galaxy begins when the jets launched from an active galactic nucleus are switched off. To study the fraction of radio galaxies in a remnant phase, we take advantage of a $8.31$\,deg$^2$ sub-region of the GAMA~23~field which comprises of surveys covering the frequency range 0.1--9\,GHz. We present a sample of 104 radio galaxies compiled from observations conducted by the Murchison Wide-field Array (216\,MHz), the Australia Square Kilometer Array Pathfinder (887\,MHz), and the Australia Telescope Compact Array (5.5\,GHz). We adopt an `absent radio core' criterion to identify 10 radio galaxies showing no evidence for an active nucleus. We classify these as new candidate remnant radio galaxies. Seven of these objects still display compact emitting regions within the lobes at 5.5\,GHz; at this frequency the emission is short-lived, implying a recent jet switch-off. On the other hand, only three show evidence of aged lobe plasma by the presence of an ultra-steep spectrum ($α<-1.2$) and a diffuse, low surface-brightness radio morphology. The predominant fraction of young remnants is consistent with a rapid fading during the remnant phase. Within our sample of radio galaxies, our observations constrain the remnant fraction to $4\%\lesssim f_{\mathrm{rem}} \lesssim 10\%$; the lower limit comes from the limiting case in which all remnant candidates with hotspots are simply active radio galaxies with faint, undetected radio cores. Finally, we model the synchrotron spectrum arising from a hotspot to show they can persist for 5--10\,Myr at 5.5\,GHz after the jets switch off -- radio emission arising from such hotspots can therefore be expected in an appreciable fraction of genuine remnants.

preprint2021arXiv

Search for neutrino counterparts of gravitational-wave events detected by LIGO and Virgo during run O2 with the ANTARES telescope

An offline search for a neutrino counterpart to gravitational-wave (GW) events detected during the second observation run (O2) of Advanced-LIGO and Advanced-Virgo performed with ANTARES data is presented. In addition to the search for long tracks induced by $ν_μ$ ($\barν_μ$) charged current interactions, a search for showering events induced by interactions of neutrinos of any flavour is conducted. The severe spatial and time coincidence provided by the gravitational-wave alert allows regions above the detector horizon to be probed, extending the ANTARES sensitivity over the entire sky. The results of this all-neutrino-flavour and all-sky time dependent analysis are presented. The search for prompt neutrino emission within $\pm$500~s around the time of six GW events yields no neutrino counterparts. Upper limits on the neutrino spectral fluence and constraints on the isotropic radiated energy are set for each GW event analysed.

preprint2021arXiv

Search for neutrinos from the tidal disruption events AT2019dsg and AT2019fdr with the ANTARES telescope

On October 1, 2019, the IceCube Collaboration detected a muon track neutrino with high probability of being of astrophysical origin, IC191001A. After a few hours, the tidal disruption event (TDE) AT2019dsg, observed by the Zwicky Transient Facility (ZTF), was indicated as the most likely counterpart of the IceCube track. More recently, the follow-up campaign of the IceCube alerts by ZTF suggested a second TDE, AT2019fdr, as a promising counterpart of another IceCube muon track candidate, IC200530A, detected on May 30, 2020. These are the second and third associations between astrophysical sources and high-energy neutrinos after the compelling identification of the blazar TXS 0506+056. Here, the search for ANTARES neutrinos from the directions of AT2019dsg and AT2019fdr using a time-integrated approach is presented. As no significant evidence for space clustering is found in the ANTARES data, upper limits on the one-flavour neutrino flux and fluence are set.

preprint2020arXiv

ANTARES and IceCube Combined Search for Neutrino Point-like and Extended Sources in the Southern Sky

A search for point-like and extended sources of cosmic neutrinos using data collected by the ANTARES and IceCube neutrino telescopes is presented. The data set consists of all the track-like and shower-like events pointing in the direction of the Southern Sky included in the nine-year ANTARES point-source analysis, combined with the through-going track-like events used in the seven-year IceCube point-source search. The advantageous field of view of ANTARES and the large size of IceCube are exploited to improve the sensitivity in the Southern Sky by a factor $\sim$2 compared to both individual analyses. In this work, the Southern Sky is scanned for possible excesses of spatial clustering, and the positions of preselected candidate sources are investigated. In addition, special focus is given to the region around the Galactic Centre, whereby a dedicated search at the location of SgrA* is performed, and to the location of the supernova remnant RXJ 1713.7-3946. No significant evidence for cosmic neutrino sources is found and upper limits on the flux from the various searches are presented.

preprint2020arXiv

Calibration database for the Murchison Widefield Array All-Sky Virtual Observatory

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.

preprint2020arXiv

Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations

We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$σ$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

preprint2020arXiv

High Cadence Optical Transient Searches using Drift Scan Imaging II: Event Rate Upper Limits on Optical Transients of Duration <21 ms and Magnitude <6.6

We have realised a simple prototype system to perform searches for short timescale optical transients, utilising the novel drift scan imaging technique described by Tingay (2020). We used two coordinated and aligned cameras, with an overlap field-of-view of approximately 3.7 sq. deg., to capture over 34000 X 5 second images during approximately 24 hours of observing. The system is sensitive to optical transients, due to an effective exposure time per pixel of 21 ms, brighter than a V magnitude of 6.6. In our 89.7 sq. deg. hr of observations we find no candidate astronomical transients, giving an upper limit to the rate of these transients of 0.8 per square degree per day, competitive with other experiments of this type. The system is triggered by reflections from satellites and various instrumental effects, which are easily identifiable due to the two camera system. The next step in the development of this promising technique is to move to a system with larger apertures and wider fields of view.

preprint2020arXiv

Monte Carlo simulations for the ANTARES underwater neutrino telescope

Monte Carlo simulations are a unique tool to check the response of a detector and to monitor its performance. For a deep-sea neutrino telescope, the variability of the environmental conditions that can affect the behaviour of the data acquisition system must be considered, in addition to a reliable description of the active parts of the detector and of the features of physics events, in order to produce a realistic set of simulated events. In this paper, the software tools used to produce neutrino and cosmic ray signatures in the telescope and the strategy developed to represent the time evolution of the natural environment and of the detector efficiency are described.

preprint2020arXiv

Observation of the cosmic ray shadow of the Sun with the ANTARES neutrino telescope

The ANTARES detector is an undersea neutrino telescope in the Mediterranean Sea. The search for point-like neutrino sources is one of the main goals of the ANTARES telescope, requiring a reliable method to evaluate the detector angular resolution and pointing accuracy. This work describes the study of the Sun "shadow" effect with the ANTARES detector. The shadow is the deficit in the atmospheric muon flux in the direction of the Sun caused by the absorption of the primary cosmic rays. This analysis is based on the data collected between 2008 and 2017 by the ANTARES telescope. The observed statistical significance of the Sun shadow detection is $3.7σ$, with an estimated angular resolution of $0.59^\circ\pm0.10^\circ$ for downward-going muons. The pointing accuracy is found to be consistent with the expectations and no evidence of systematic pointing shifts is observed.

preprint2020arXiv

Sensitivity of a Low-Frequency Polarimetric Radio Interferometer

Aims: This paper aims to derive an expression for the sensitivity of a polarimetric radio interferometer that is valid for all-sky observations of arbitrarily polarized sources, with neither a restriction on FoV nor with any a priori assumption regarding the polarization state of the source. We verify the resulting formula with an all-sky observation using the Murchison Widefield Array (MWA) telescope. Methods: The sensitivity expression is developed from first principles by applying the concept of System Equivalent Flux Density (SEFD) to a polarimetric radio interferometer not by computing $A_e/T_{sys}$. The SEFD is calculated from the standard deviation of the noisy flux density estimate for a target source due to system noise. Results: The SEFD for a polarimetric radio interferometer is generally not $1/\sqrt{2}$ of a single-polarized interferometer as is often assumed for narrow FoV. This assumption can lead to significant errors for a dual-polarized dipole based system, which is common in low-frequency radio astronomy: up to $\sim 15\%$ for a zenith angle (ZA) coverage of $45^\circ$, and up to $\sim45\%$ for $60^\circ$ coverage. The worst case errors occur in the diagonal planes of the dipole for very wide FoV. This is demonstrated through theory, simulation and observations. Furthermore, using the resulting formulation, calculation of the off-zenith sensitivity is straightforward and unambiguous. Conclusions: For wide FoV observations pertinent to low-frequency radio interferometer such as the SKA-Low, the narrow FoV and the single-polarized sensitivity expressions are not correct and should be replaced by the formula derived in this paper.

preprint2020arXiv

The SKA Particle Array Prototype: The First Particle Detector at the Murchison Radio-astronomy Observatory

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

preprint2019arXiv

A VOEvent based automatic trigger system for the Murchison Widefield Array

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.

preprint2019arXiv

Science with the Murchison Widefield Array: Phase I Results and Phase II Opportunities

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.

preprint2016arXiv

A Fast Radio Burst Host Galaxy

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $Ω_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

preprint2016arXiv

A High Reliability Survey of Discrete Epoch of Reionization Foreground Sources in the MWA EoR0 Field

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.

preprint2016arXiv

A Large Scale, Low Frequency Murchison Widefield Array Survey of Galactic HII regions between 260< l <\340

We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spectral indices of these HII regions. By identifying the point at which HII regions transition from the optically thin to thick regime we derive the physical properties including the electron density, ionised gas mass and ionising photon flux, towards 61 HII regions. This catalogue of HII regions represents the most extensive and uniform low frequency survey of HII regions in the Galaxy to date.

preprint2016arXiv

A multi-resolution, multi-epoch low radio frequency survey of the Kepler K2 mission Campaign 1 field

We present the first dedicated radio continuum survey of a Kepler K2 mission field, Field 1 covering the North Galactic Cap. The survey is wide field, contemporaneous, multi-epoch, and multi-resolution in nature and was conducted at low radio frequencies between 140 and 200 MHz. The multi-epoch and ultra wide field (but relatively low resolution) part of the survey was provided by 15 nights of observation with the Murchison Widefield Array (MWA) over a period of approximately a month, contemporaneous with K2 observations of the field. The multi-resolution aspect of the survey was provided by the low resolution (4') MWA imaging, complemented by non-contemporaneous but much higher resolution (20") observations using the Giant Metrewave Radio Telescope (GMRT). The survey is therefore sensitive to the details of radio structures across a wide range of angular scales. Consistent with other recent low radio frequency surveys, no significant radio transients or variables were detected in the survey. The resulting source catalogs consist of 1,085 and 1,468 detections in the two MWA observation bands (centered at 154 and 185 MHz, respectively) and 7,445 detections in the GMRT observation band (centered at 148 MHz), over 314 square degrees. The survey is presented as a significant resource for multi-wavelength investigations of the more than 21,000 target objects in the K2 field. We briefly examine our survey data against K2 target lists for dwarf star types (stellar types M and L) that have been known to produce radio flares.

preprint2016arXiv

An Opportunistic Search for ExtraTerrestrial Intelligence (SETI) with the Murchison Widefield Array

A spectral line image cube generated from 115 minutes of MWA data that covers a field of view of 400 sq. deg. around the Galactic Centre is used to perform the first Search for ExtraTerrestrial Intelligence (SETI) with the Murchison Widefield Array. Our work constitutes the first modern SETI experiment at low radio frequencies, here between 103 and 133 MHz, paving the way for large-scale searches with the MWA and, in the future, the low frequency Square Kilometre Array. Limits of a few hundred mJy/beam for narrow band emission (10 kHz) are derived from our data, across our 400 sq. deg. field of view. Within this field, 45 exoplanets in 38 planetary systems are known. We extract spectra at the locations of these systems from our image cube, to place limits on the presence of narrow line emission from these systems. We then derive minimum isotropic transmitter powers for these exoplanets; a small handful of the closest objects (10s of pc) yield our best limits of order $10^{14}$ W (Equivalent Isotropic Radiated Power: EIRP). These limits lie above the highest power directional transmitters near these frequencies currently operational on Earth. A SETI experiment with the MWA covering the full accessible sky and its full frequency range would require approximately one month of observing time. The MWA frequency range, its Southern Hemisphere location on an extraordinarily radio quiet site, its very large field of view, and its high sensitivity make it a unique facility for SETI.

preprint2016arXiv

Beamforming Errors in Murchison Widefield Array Phased Array Antennas and their effects on Epoch of Reionization Science

Accurate antenna beam models are critical for radio observations aiming to isolate the redshifted 21cm spectral line emission from the Dark Ages and the Epoch of Reionization and unlock the scientific potential of 21cm cosmology. Past work has focused on characterizing mean antenna beam models using either satellite signals or astronomical sources as calibrators, but antenna-to-antenna variation due to imperfect instrumentation has remained unexplored. We characterize this variation for the Murchison Widefield Array (MWA) through laboratory measurements and simulations, finding typical deviations of order +/- 10-20% near the edges of the main lobe and in the sidelobes. We consider the ramifications of these results for image- and power spectrum-based science. In particular, we simulate visibilities measured by a 100m baseline and find that using an otherwise perfect foreground model, unmodeled beamforming errors severely limit foreground subtraction accuracy within the region of Fourier space contaminated by foreground emission (the "wedge"). This region likely contains much of the cosmological signal, and accessing it will require measurement of per-antenna beam patterns. However, unmodeled beamforming errors do not contaminate the Fourier space region expected to be free of foreground contamination (the "EOR window"), showing that foreground avoidance remains a viable strategy.

preprint2016arXiv

Capturing the electromagnetic counterparts of binary neutron star mergers through low latency gravitational wave triggers

We investigate the prospects for joint low-latency gravitational wave (GW) detection and prompt electromagnetic (EM) follow-up observations of coalescing binary neutron stars (BNSs). For BNS mergers associated with short duration gamma-ray bursts (SGRBs), we for the first time evaluate the feasibility of rapid EM follow-ups to capture the prompt emission, early engine activity or reveal any potential by-products such as magnetars or fast radio bursts. To achieve our goal, we first simulate a population of coalescing BNSs using realistic distributions of source parameters and estimate the detectability and localisation efficiency at different times before merger. We then use a selection of facilities with GW follow-up agreements in place, from low-frequency radio to high energy $γ$-ray to assess the prospects of prompt follow-up. We quantify our assessment using observational SGRB flux data extrapolated to be within the horizon distances of the advanced GW interferometric detectors LIGO and Virgo and to the prompt phase immediately following the binary merger. Our results illustrate that while challenging, breakthrough multi-messenger science is possible with EM follow-up facilities with fast responses and wide fields-of-view. We demonstrate that the opportunity to catch the prompt stage ($<$ 5s) of SGRBs, can be enhanced by speeding up the detection pipelines of both GW observatories and EM follow-up facilities. We further show that the addition of an Australian instrument to the optimal detector network could possibly improve the angular resolution by a factor of two and thereby contribute significantly to GW-EM multi-messenger astronomy.

preprint2016arXiv

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

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

preprint2016arXiv

Fast response electromagnetic follow-ups from low latency GW triggers

We investigate joint low-latency gravitational wave (GW) detection and prompt electromagnetic (EM) follow-up observations of coalescing binary neutron stars (BNSs). Assuming that BNS mergers are associated with short duration gamma ray bursts (SGRBs), we evaluate if rapid EM follow-ups can capture the prompt emission, early engine activity or reveal any potential by-products such as magnetars or fast radio bursts. To examine the expected performance of extreme low-latency search pipelines, we simulate a population of coalescing BNSs and use these to estimate the detectability and localisation efficiency at different times before merger. Using observational SGRB flux data corrected to the range of the advanced GW interferometric detectors, we determine what EM observations could be achieved from low-frequency radio up to high energy $γ$-ray. We show that while challenging, breakthrough multi-messenger science is possible through low latency pipelines.

preprint2016arXiv

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

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

preprint2016arXiv

First Season MWA EoR Power Spectrum Results at Redshift 7

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

preprint2016arXiv

High-energy sources at low radio frequency: the Murchison Widefield Array view of Fermi blazars

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

preprint2016arXiv

Ionospheric Modelling using GPS to Calibrate the MWA. II: Regional ionospheric modelling using GPS and GLONASS to estimate ionospheric gradients

We estimate spatial gradients in the ionosphere using the Global Positioning System (GPS) and GLONASS (Russian global navigation system) observations, utilising data from multiple GPS stations in the vicinity of Murchison Radio-astronomy Observatory (MRO). In previous work the ionosphere was characterised using a single-station to model the ionosphere as a single layer of fixed height and this was compared with ionospheric data derived from radio astronomy observations obtained from the Murchison Widefield Array (MWA). Having made improvements to our data quality (via cycle slip detection and repair) and incorporating data from the GLONASS system, we now present a multi-station approach. These two developments significantly improve our modelling of the ionosphere. We also explore the effects of a variable-height model. We conclude that modelling the small-scale features in the ionosphere that have been observed with the MWA will require a much denser network of Global Navigation Satellite System (GNSS) stations than is currently available at the MRO.

preprint2016arXiv

Limits on Fast Radio Bursts and other transient sources at 182 MHz using the Murchison Widefield Array

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.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Low frequency observations of linearly polarized structures in the interstellar medium near the south Galactic pole

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.

preprint2016arXiv

Murchison Widefield Array Limits on Radio Emission from ANTARES Neutrino Events

We present a search, using the Murchison Widefield Array (MWA), for electromagnetic counterparts to two candidate high energy neutrino events detected by the ANTARES neutrino telescope in 2013 November and 2014 March. These events were selected by ANTARES because they are consistent, within 0.4 degrees, with the locations of galaxies within 20 Mpc of Earth. Using MWA archival data at frequencies between 118 and 182 MHz, taken ~20 days prior to, at the same time as, and up to a year after the neutrino triggers, we look for transient or strongly variable radio sources consistent with the neutrino positions. No such counterparts are detected, and we set a 5 sigma upper limit for low-frequency radio emission of ~1E37 erg/s for progenitors at 20 Mpc. If the neutrino sources are instead not in nearby galaxies, but originate in binary neutron star coalescences, our limits place the progenitors at z > 0.2. While it is possible, due to the high background from atmospheric neutrinos, that neither event is astrophysical, the MWA observations are nevertheless among the first to follow up neutrino candidates in the radio, and illustrate the promise of wide-field instruments like MWA to detect electromagnetic counterparts to such events.

preprint2016arXiv

Optical and radio astrometry of the galaxy associated with FRB150418

A fading radio source, coincident in time and position with the fast radio burst FRB150418, has been associated with the galaxy WISE J071634.59-190039.2. Subsequent observations of this galaxy have revealed that it contains a persistent, but variable, radio source. We present e-MERLIN, VLBA, and ATCA radio observations and Subaru optical observations of WISE J071634.59-190039.2 and find that the persistent radio source is unresolved and must be compact (<0.01 kpc), and that its location is consistent with the optical centre of the galaxy. We conclude that it is likely that WISE J071634.59-190039.2 contains a weak radio AGN.

preprint2016arXiv

Parametrising Epoch of Reionization foregrounds: A deep survey of low-frequency point-source spectra with the MWA

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.

preprint2016arXiv

Radio light curve of the galaxy possibly associated with FRB150418

We present observations made with the Australia Telescope Compact Array (ATCA), the Jansky Very Large Array (JVLA) and the Giant Metre-Wave Telescope of the radio source within the galaxy WISE~J071634.59-190039.2, claimed to be host of FRB~150418 by Keane et al. (2016). We have established a common flux density scale between the ATCA and JVLA observations, the main result of which is to increase the flux densities obtained by Keane et al. At a frequency of 5.5 GHz, the source has a mean flux density of 140uJy and is variable on short timescales with a modulation index of 0.36. Statistical analysis of the flux densities shows that the variations seen are consistent with refractive interstellar scintillation of the weak active galactic nucleus at the centre of the galaxy. It may therefore be the case that the FRB and the galaxy are not associated. However, taking into account the rarity of highly variable sources in the radio sky, and our lack of knowledge of the progenitors of FRBs as a class, the association between WISE~J071634.59-190039.2 and FRB~150418 remains a possibility.

preprint2016arXiv

SN 1978K: An evolved supernova outside our Local Group detected at millimetre wavelengths

Supernova 1978K is one of the oldest-known examples of the class of Type IIn supernovae that show evidence for strong interaction between the blast wave and a dense, pre-existing circumstellar medium. Here we report detections of SN 1978K at both 34 GHz and 94 GHz, making it only the third extragalactic supernova (after SN 1987A and SN 1996cr) to be detected at late-times at these frequencies. We find SN 1978K to be >400 times more luminous than SN 1987A at millimetre wavelengths in spite of the roughly nine year difference in ages, highlighting the risk in adopting SN 1987A as a template for the evolution of core-collapse supernovae in general. Additionally, from new VLBI observations at 8.4 GHz, we measure a deconvolved diameter for SN 1978K of ~5 milli-arcsec, and a corresponding average expansion velocity of <1500 km/s. These observations provide independent evidence of an extremely dense circumstellar medium surrounding the progenitor star.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

The 154 MHz radio sky observed by the Murchison Widefield Array: noise, confusion and first source count analyses

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.

preprint2016arXiv

The Importance of Wide-field Foreground Removal for 21 cm Cosmology: A Demonstration With Early MWA Epoch of Reionization Observations

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.

preprint2016arXiv

The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology

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

preprint2016arXiv

Time-domain and spectral properties of pulsars at 154 MHz

We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < δ < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715.

preprint2015arXiv

A Deep Search for Prompt Radio Emission from the Short GRB 150424A With The Murchison Widefield Array

We present a search for prompt radio emission associated with the short-duration gamma-ray burst (GRB) 150424A using the Murchison Widefield Array (MWA) at frequencies from 80-133 MHz. Our observations span delays of 23 s-30 min after the GRB, corresponding to dispersion measures of 100-7700 pc/cm^3. We see no excess flux in images with timescales of 4 s, 2 min, or 30 min, and set a 3 sigma flux density limit of 3.0 Jy at 132 MHz on the shortest timescales: some of the most stringent limits to date on prompt radio emission from any type of GRB. We use these limits to constrain a number of proposed models for coherent emission from short-duration GRBs, although we show that our limits are not particularly constraining for fast radio bursts because of reduced sensitivity for this pointing. Finally, we discuss the prospects for using the MWA to search for prompt radio emission from gravitational wave transients and find that while the flux density and luminosity limits are likely to be very constraining, the latency of the gravitational wave alert may limit the robustness of any conclusions.

preprint2015arXiv

A search for Fast Radio Bursts at low frequencies with Murchison Widefield Array high time resolution imaging

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$σ$ 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 $α<-1.2$ ($S\proptoν^α$). 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.

preprint2015arXiv

An analysis of the halo and relic radio emission from Abell 3376 from Murchison Widefield Array observations

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

preprint2015arXiv

Antenna array characterization via radio interferometry observation of astronomical sources

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

preprint2015arXiv

BIGHORNS - Broadband Instrument for Global HydrOgen ReioNisation Signal

The redshifted 21cm line of neutral hydrogen (HI), potentially observable at low radio frequencies (~50-200 MHz), should be a powerful probe of the physical conditions of the inter-galactic medium during Cosmic Dawn and the Epoch of Reionisation (EoR). The sky-averaged HI signal is expected to be extremely weak (~100 mK) in comparison to the foreground of up to 10000 K at the lowest frequencies of interest. The detection of such a weak signal requires an extremely stable, well characterised system and a good understanding of the foregrounds. Development of a nearly perfectly (~mK accuracy) calibrated total power radiometer system is essential for this type of experiment. We present the BIGHORNS (Broadband Instrument for Global HydrOgen ReioNisation Signal) experiment which was designed and built to detect the sky-averaged HI signal from the EoR at low radio frequencies. The BIGHORNS system is a mobile total power radiometer, which can be deployed in any remote location in order to collect radio-interference (RFI) free data. The system was deployed in remote, radio quiet locations in Western Australia and low RFI sky data have been collected. We present a description of the system, its characteristics, details of data analysis and calibration. We have identified multiple challenges to achieving the required measurement precision, which triggered two major improvements for the future system.

preprint2015arXiv

Broadband Spectral Modeling of the Extreme Gigahertz-Peaked Spectrum Radio Source PKS B0008-421

We present broadband observations and spectral modeling of PKS B0008-421, and identify it as an extreme gigahertz-peaked spectrum (GPS) source. PKS B0008-421 is characterized by the steepest known spectral slope below the turnover, close to the theoretical limit of synchrotron self-absorption, and the smallest known spectral width of any GPS source. Spectral coverage of the source spans from 0.118 to 22 GHz, which includes data from the Murchison Widefield Array and the wide bandpass receivers on the Australia Telescope Compact Array. We have implemented a Bayesian inference model fitting routine to fit the data with various absorption models. We find that without the inclusion of a high-frequency exponential break the absorption models can not accurately fit the data, with significant deviations above and below the peak in the radio spectrum. The addition of a high-frequency break provides acceptable spectral fits for the inhomogeneous free-free absorption and double-component synchrotron self-absorption models, with the inhomogeneous free-free absorption model statistically favored. The requirement of a high-frequency spectral break implies that the source has ceased injecting fresh particles. Additional support for the inhomogeneous free-free absorption model as being responsible for the turnover in the spectrum is given by the consistency between the physical parameters derived from the model fit and the implications of the exponential spectral break, such as the necessity of the source being surrounded by a dense ambient medium to maintain the peak frequency near the gigahertz region. The discovery of PKS B0008-421 suggests that the next generation of low radio frequency surveys could reveal a large population of GPS sources that have ceased activity, and that a portion of the ultra-steep spectrum source population could be composed of these GPS sources in a relic phase.

preprint2015arXiv

Characterization of a Low-Frequency Radio Astronomy Prototype Array in Western Australia

We report characterization results for an engineering prototype of a next-generation low-frequency radio astronomy array. This prototype, which we refer to as the Aperture Array Verification System 0.5 (AAVS0.5), is a sparse pseudo-random array of 16 log-periodic antennas designed for 70-450 MHz. It is co-located with the Murchison Widefield Array (MWA) at the Murchison Radioastronomy Observatory (MRO) near the Australian Square Kilometre Array (SKA) core site. We characterize the AAVS0.5 using two methods: in-situ radio interferometry with astronomical sources and an engineering approach based on detailed full-wave simulation. In-situ measurement of the small prototype array is challenging due to the dominance of the Galactic noise and the relatively weaker calibration sources easily accessible in the southern sky. The MWA, with its 128 "tiles" and up to 3 km baselines, enabled in-situ measurement via radio interferometry. We present array sensitivity and beam pattern characterization results and compare to detailed full-wave simulation. We discuss areas where differences between the two methods exist and offer possibilities for improvement. Our work demonstrates the value of the dual astronomy-simulation approach in upcoming SKA design work.

preprint2015arXiv

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

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

preprint2015arXiv

Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

preprint2015arXiv

Foregrounds in Wide-Field Redshifted 21 cm Power Spectra

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

preprint2015arXiv

Ionospheric modelling using GPS to calibrate the MWA. 1: Comparison of first order ionospheric effects between GPS models and MWA observations

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

preprint2015arXiv

Measuring Phased-Array Antenna Beampatterns with High Dynamic Range for the Murchison Widefield Array using 137 MHz ORBCOMM Satellites

Detection of the fluctuations in 21 cm line emission from neutral hydrogen during the Epoch of Reionization in thousand hour integrations poses stringent requirements on calibration and image quality, both of which necessitate accurate primary beam models. The Murchison Widefield Array (MWA) uses phased array antenna elements which maximize collecting area at the cost of complexity. To quantify their performance, we have developed a novel beam measurement system using the 137 MHz ORBCOMM satellite constellation and a reference dipole antenna. Using power ratio measurements, we measure the {\it in situ} beampattern of the MWA antenna tile relative to that of the reference antenna, canceling the variation of satellite flux or polarization with time. We employ angular averaging to mitigate multipath effects (ground scattering), and assess environmental systematics with a null experiment in which the MWA tile is replaced with a second reference dipole. We achieve beam measurements over 30 dB dynamic range in beam sensitivity over a large field of view (65\% of the visible sky), far wider and deeper than drift scans through astronomical sources allow. We verify an analytic model of the MWA tile at this frequency within a few percent statistical scatter within the full width at half maximum. Towards the edges of the main lobe and in the sidelobes, we measure tens of percent systematic deviations. We compare these errors with those expected from known beamforming errors.

preprint2015arXiv

Murchison Widefield Array Observations of Anomalous Variability: A Serendipitous Night-time Detection of Interplanetary Scintillation

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.

preprint2015arXiv

Power spectrum analysis of ionospheric fluctuations with the Murchison Widefield Array

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

preprint2015arXiv

Quantifying ionospheric effects on time-domain astrophysics with the Murchison Widefield Array

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

preprint2015arXiv

Real-time imaging of density ducts between the plasmasphere and ionosphere

Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures.

preprint2015arXiv

Scintillation arcs in low-frequency observations of the timing-array millisecond pulsar J0437-4715

Low-frequency observations of pulsars provide a powerful means for probing the microstructure in the turbulent interstellar medium (ISM). Here we report on high-resolution dynamic spectral analysis of our observations of the timing-array millisecond pulsar J0437-4715 with the Murchison Widefield Array (MWA), enabled by our recently commissioned tied-array beam processing pipeline for voltage data recorded from the high time resolution mode of the MWA. A secondary spectral analysis reveals faint parabolic arcs, akin to those seen in high-frequency observations of pulsars with the Green Bank and Arecibo telescopes. Data from Parkes observations at a higher frequency of 732 MHz reveal a similar parabolic feature, with a curvature that scales approximately as the square of the observing wavelength ($λ^2$) to the MWA's frequency of 192 MHz. Our analysis suggests that scattering toward PSR J0437-4715 predominantly arises from a compact region about 115 pc from the Earth, which matches well with the expected location of the edge of the Local Bubble that envelopes the local Solar neighborhood. As well as demonstrating new and improved pulsar science capabilities of the MWA, our analysis underscores the potential of low-frequency pulsar observations for gaining valuable insights into the local ISM and for characterizing the ISM toward timing-array pulsars.

preprint2015arXiv

Simultaneous Observations of Giant Pulses from the Crab Pulsar, with the Murchison Widefield Array and Parkes Radio Telescope: Implications for the Giant Pulse Emission Mechanism

We report on observations of giant pulses from the Crab pulsar performed simultaneously with the Parkes radio telescope and the incoherent combination of the Murchison Widefield Array (MWA) antenna tiles. The observations were performed over a duration of approximately one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center frequency of 193 MHz with 15 MHz bandwidth at the MWA. Our analysis has led to the detection of 55 giant pulses at the MWA and 2075 at Parkes above a threshold of 3.5$σ$ and 6.5$σ$ respectively. We detected 51$\%$ of the MWA giant pulses at the Parkes radio telescope, with spectral indices in the range of $-3.6>α> -4.9$ ($S_{\rm ν} \propto ν^α$). We present a Monte Carlo analysis supporting the conjecture that the giant pulse emission in the Crab is intrinsically broadband, the less than $100\%$ correlation being due to the relative sensitivities of the two instruments and the width of the spectral index distribution. Our observations are consistent with the hypothesis that the spectral index of giant pulses is drawn from normal distribution of standard deviation 0.6, but with a mean that displays an evolution with frequency from -3.00 at 1382 MHz, to -2.85 at 192 MHz.

preprint2015arXiv

The High Time and Frequency Resolution Capabilities of the Murchison Widefield Array

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

preprint2015arXiv

The low-frequency environment of the Murchison Widefield Array: radio-frequency interference analysis and mitigation

The Murchison Widefield Array (MWA) is a new low-frequency interferometric radio telescope built in Western Australia at one of the locations of the future Square Kilometre Array (SKA). We describe the automated radio-frequency interference (RFI) detection strategy implemented for the MWA, which is based on the AOFlagger platform, and present 72-231-MHz RFI statistics from 10 observing nights. RFI detection removes 1.1% of the data. RFI from digital TV (DTV) is observed 3% of the time due to occasional ionospheric or atmospheric propagation. After RFI detection and excision, almost all data can be calibrated and imaged without further RFI mitigation efforts, including observations within the FM and DTV bands. The results are compared to a previously published Low-Frequency Array (LOFAR) RFI survey. The remote location of the MWA results in a substantially cleaner RFI environment compared to LOFAR's radio environment, but adequate detection of RFI is still required before data can be analysed. We include specific recommendations designed to make the SKA more robust to RFI, including: the availability of sufficient computing power for RFI detection; accounting for RFI in the receiver design; a smooth band-pass response; and the capability of RFI detection at high time and frequency resolution (second and kHz-scale respectively).

preprint2015arXiv

The multi-frequency parsec-scale structure of PKS 2254-367 (IC 1459): a luminosity-dependent break in morphology for the precursors of radio galaxies?

We present the first multi-frequency VLBI images of PKS 2254-367, a Giga-hertz Peaked Spectrum (GPS) radio source hosted by the nearby galaxy IC 1459 (D=20.5 Mpc). PKS 2254-367 and the radio source in NGC 1052 (PKS 0238-084; D=17.2 Mpc) are the two closest GPS radio sources to us, far closer than the next closest example, PKS 1718-649 (D=59 Mpc). As such, IC 1459 and NGC 1052 offer opportunities to study the details of the pc-scale radio sources as well as the environments that the radio sources inhabit, across the electromagnetic spectrum. Given that some models for the origin and evolution of GPS radio sources require a strong connection between the radio source morphology and the gaseous nuclear environment, such opportunities for detailed study are important. Our VLBI images of PKS 2254-367 show that the previously identified similarities between IC 1459 and NGC 1052 continue onto the pc-scale. Both compact radio sources appear to have symmetric jets of approximately the same luminosity, much lower than typically noted in compact double GPS sources. Similarities between PKS 2254-367 and NGC 1052, and differences with respect to other GPS galaxies, lead us to speculate that a sub-class of GPS radio sources, with low luminosity and with jet-dominated morphologies, exists and would be largely absent from radio source surveys with ~1 Jy flux density cutoffs. We suggest that this possible low-luminosity, jet-dominated population of GPS sources could be an analog of the FR-I radio galaxies, with the higher luminosity lobe-dominated GPS sources being the analog of the FR-II radio galaxies.

preprint2015arXiv

The Murchison Widefield Array Correlator

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

preprint2014arXiv

First Look Murchison Widefield Array observations of Abell 3667

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

preprint2014arXiv

Limits on low frequency radio emission from southern exoplanets with the Murchison Widefield Array

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

preprint2014arXiv

Modelling of the Spectral Energy Distribution of Fornax A: Leptonic and Hadronic Production of High Energy Emission from the Radio Lobes

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

preprint2014arXiv

Observing the Sun with the Murchison Widefield Array

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

preprint2014arXiv

The first Murchison Widefield Array low frequency radio observations of cluster scale non-thermal emission: the case of Abell 3667

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

preprint2014arXiv

The importance of jet bending in gamma-ray AGN - revisited

We investigate the hypothesis that γ-ray-quiet AGN have a larger tendency for jet bending than γ-ray-loud AGN, revisiting the analysis of Tingay, Murphy & Edwards (1998). We perform a statistical analysis using a large sample of 351 radio-loud AGN along with γ-ray identifications from the Fermi Large Area Telescope (LAT). Our results show no statistically significant differences in jet-bending properties between γ-ray-loud and γ-ray-quiet populations, indicating that jet bending is not a significant factor for γ-ray detection in AGN.

preprint2014arXiv

The low-frequency characteristics of PSR J0437-4715 observed with the Murchison Widefield Array

We report on the detection of the millisecond pulsar PSR J0437-4715 with the Murchison Widefield Array (MWA) at a frequency of 192 MHz. Our observations show rapid modulations of pulse intensity in time and frequency that arise from diffractive scintillation effects in the interstellar medium (ISM), as well as prominent drifts of intensity maxima in the time-frequency plane that arise from refractive effects. Our analysis suggests that the scattering screen is located at a distance of $\sim$80-120 pc from the Sun, in disagreement with a recent claim that the screen is closer ($\sim$10 pc). Comparisons with higher frequency data from Parkes reveals a dramatic evolution of the pulse profile with frequency, with the outer conal emission becoming comparable in strength to that from the core and inner conal regions. As well as demonstrating high time resolution science capabilities currently possible with the MWA, our observations underscore the potential to conduct low-frequency investigations of timing-array millisecond pulsars, which may lead to increased sensitivity for the detection of nanoHertz gravitational waves via the accurate characterisation of ISM effects.

preprint2014arXiv

The Spectral Variability of the GHz-Peaked Spectrum Radio Source PKS 1718-649 and a Comparison of Absorption Models

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

preprint2013arXiv

A 189 MHz, 2400 square degree polarization survey with the Murchison Widefield Array 32-element prototype

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

preprint2013arXiv

A survey for transients and variables with the Murchison Widefield Array 32-tile prototype at 154 MHz

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

preprint2013arXiv

On the detection and tracking of space debris using the Murchison Widefield Array. I. Simulations and test observations demonstrate feasibility

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

preprint2013arXiv

Performing a stellar autopsy using the radio-bright remnant of SN 1996cr

We present archival radio observations of SN1996cr from ATCA and MOST, and model its radio lightcurves using X-ray constrained hydrodynamical simulations of the interaction between the SN ejecta and the CSM. The early radio data show signs of free-free absorption (FFA), which decreases gradually and is minimal above 1.4 GHz after day ~3000. FFA optical depth constraints provide estimates of the CSM free electron density, which allows insight into the ionisation of SN1996cr's CSM and offers a test on the density distribution adopted by the hydrodynamical simulation. The intrinsic spectral index of the radiation shows evidence for spectral flattening, which is characterised by alpha = 0.852 +/- 0.002 at day 3000 and a decay rate of d_alpha = -0.014 +/- 0.001 yr^-1. The similarity in the spectral flattening of SN1987A, SN1993J, and SN1996cr suggests this may be a relatively common feature of SNe/CSM shocks. We adopt this spectral index variation to model the synchrotron radio emission of the shock, and consider several scalings that relate the parameters of the hydrodynamical simulation to the magnetic field and electron distribution. The simulated light curves match the large-scale features of the observed light curves, but fail to match certain tightly constraining sections. This suggests that simple energy density scalings may not be able to account for the complexities of the true physical processes at work, or alternatively, that the parameters of the simulation require modification in order to accurately represent SN1996cr's CSM.

preprint2013arXiv

Radio observations of GRB 100418a: Test of an energy injection model explaining long-lasting GRB afterglows

We present the results of our radio observational campaign on GRB 100418a, for which we used the Australia Telescope Compact Array (ATCA), Very Large Array (VLA) and the Very Long Baseline Array (VLBA). GRB 100418a was a peculiar GRB with unusual X-ray and optical afterglow profiles featuring a plateau phase with a very shallow rise. This observed plateau phase was believed to be due to a continued energy injection mechanism which powered the forward shock, giving rise to an unusual and long-lasting afterglow. The radio afterglow of GRB 100418a was detectable several weeks after the prompt emission. We conducted long-term monitoring observations of the afterglow and attempted to test the energy injection model advocating that the continuous energy injection is due to shells of material moving at a wide range of Lorentz factors. We obtained an upper limit of $γ$ $<$ 7 for the expansion rate of the GRB 100418a radio afterglow, indicating that the range-of-Lorentz factor model could only be applicable for relatively slow moving ejecta. A preferred explanation could be that continued activity of the central engine may have powered the long-lasting afterglow.

preprint2013arXiv

The AuScope Geodetic VLBI Array

The AuScope geodetic Very Long Baseline Interferometry array consists of three new 12 m radio telescopes and a correlation facility in Australia. The telescopes at Hobart (Tasmania), Katherine (Northern Territory) and Yarragadee (Western Australia) are co-located with other space geodetic techniques including Global Navigation Satellite Systems (GNSS) and gravity infrastructure, and in the case of Yarragadee, Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) facilities. The correlation facility is based in Perth (Western Australia). This new facility will make significant contributions to improving the densification of the International Celestial Reference Frame in the Southern Hemisphere, and subsequently enhance the International Terrestrial Reference Frame through the ability to detect and mitigate systematic error. This, combined with the simultaneous densification of the GNSS network across Australia will enable the improved measurement of intraplate deformation across the Australian tectonic plate. In this paper we present a description of this new infrastructure and present some initial results, including telescope performance measurements and positions of the telescopes in the International Terrestrial Reference Frame. We show that this array is already capable of achieving centimetre precision over typical long-baselines and that network and reference source systematic effects must be further improved to reach the ambitious goals of VLBI2010.

preprint2013arXiv

The giant lobes of Centaurus A observed at 118 MHz with the Murchison Widefield Array

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

preprint2013arXiv

The Murchison Widefield Array: solar science with the low frequency SKA Precursor

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

preprint2013arXiv

The unusual afterglow of the Gamma-Ray Burst 100621A

In order to constrain the broad-band spectral energy distribution of the afterglow of GRB 100621A, dedicated observations were performed in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at the 2.2m MPG/ESO telescope, in the sub-millimeter band with the large bolometer array LABOCA at APEX, and at radio frequencies with ATCA. Utilizing also Swift X-ray observations, we attempt an interpretation of the observational data within the fireball scenario. The afterglow of GRB 100621A shows a very complex temporal as well as spectral evolution. We identify three different emission components, the most spectacular one causing a sudden intensity jump about one hour after the prompt emission. The spectrum of this component is much steeper than the canonical afterglow. We interpret this component using the prescription of Vlasis et al. (2011) for a two-shell collision after the first shell has been decelerated by the circumburst medium. We use the fireball scenario to derive constraints on the microphysical parameters of the first shell. Long-term energy injection into a narrow jet seems to provide an adequate description. Another noteworthy result is the large ($A_V$ = 3.6 mag) line-of-sight host extinction of the afterglow in an otherwise extremely blue host galaxy.

preprint2013arXiv

The unusual multiwavelength properties of the gamma-ray source PMNJ1603-4904

We investigate the nature and classification of PMNJ1603-4904, a bright radio source close to the Galactic plane, which is associated with one of the brightest hard-spectrum gamma-ray sources detected by Fermi/LAT. It has previously been classified as a low-peaked BL Lac object based on its broadband emission and the absence of optical emission lines. Optical measurements, however, suffer strongly from extinction and the absence of pronounced short-time gamma-ray variability over years of monitoring is unusual for a blazar. We are combining new and archival multiwavelength data in order to reconsider the classification and nature of this unusual gamma-ray source. For the first time, we study the radio morphology at 8.4GHz and 22.3GHz, and its spectral properties on milliarcsecond (mas) scales, based on VLBI observations from the TANAMI program. We combine the resulting images with multiwavelength data in the radio, IR, optical/UV, X-ray, and gamma-ray regimes. PMNJ1603-4904 shows a symmetric brightness distribution at 8.4GHz on mas-scales, with the brightest, and most compact component in the center of the emission region. The morphology is reminiscent of a Compact Symmetric Object (CSO). Such objects have been predicted to produce gamma-ray emission but have not been detected as a class by Fermi/LAT so far. Sparse (u, v)-coverage at 22.3GHz prevents an unambiguous modeling of the source morphology. IR measurements reveal an excess in the spectral energy distribution (SED), which can be modeled with a blackbody with a temperature of about 1600K, and which is usually not present in blazar SEDs. The VLBI data and the shape of the SED challenge the current blazar classification. PMNJ1603-4904 seems to be either a highly peculiar BL Lac object or a misaligned jet source. In the latter case, the intriguing VLBI structure opens room for a possible classification as a gamma-ray bright CSO.

preprint2013arXiv

VLBI imaging throughout the primary beam using accurate UV shifting

For Very Long Baseline Interferometry (VLBI), the fringe spacing is extremely narrow compared to the field of view imposed by the primary beam of each element. This means that an extremely large number of resolution units can potentially be imaged from a single observation. We implement and test a technique for efficiently and accurately imaging large VLBI datasets. The DiFX software correlator is used to generate a dataset with extremely high time and frequency resolution. This large dataset is then transformed and averaged multiple times to generate many smaller datasets, each with a phase centre located at a different area of interest. Results of an 8.4 GHz four-station VLBI observation of a field containing multiple sources are presented. Observations of the calibrator 3C345 were used for preliminary tests of accuracy of the shifting algorithm. A high level of accuracy was achieved, making the method suitable even for the most demanding astrometric VLBI observations. One target source (1320+299A) was detected and was used as a phase-reference calibrator in searching for further detections. An image containing 13 billion pixels was constructed by independently imaging 782 visibility datasets covering the entire primary beam of the array. Current implementations of this algorithm and possible future developments in VLBI data analysis are discussed.

preprint2012arXiv

A new layout optimization technique for interferometric arrays, applied to the MWA

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

preprint2012arXiv

Low Frequency Observations of the Moon with the Murchison Widefield Array

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

preprint2012arXiv

Realisation of a low frequency SKA Precursor: The Murchison Widefield Array

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

preprint2012arXiv

The EoR Sensitivity of the Murchison Widefield Array

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

preprint2012arXiv

The First Very Long Baseline Interferometric SETI Experiment

The first Search for Extra-Terrestrial Intelligence (SETI) conducted with Very Long Baseline Interferometry (VLBI) is presented. By consideration of the basic principles of interferometry, we show that VLBI is efficient at discriminating between SETI signals and human generated radio frequency interference (RFI). The target for this study was the star Gliese 581, thought to have two planets within its habitable zone. On 2007 June 19, Gliese 581 was observed for 8 hours at 1230-1544 with the Australian Long Baseline Array. The dataset was searched for signals appearing on all interferometer baselines above five times the noise limit. A total of 222 potential SETI signals were detected and by using automated data analysis techniques, were ruled out as originating from the Gliese 581 system. From our results we place an upper limit of 7 MW/Hz on the power output of any isotropic emitter located in the Gliese 581 system, within this frequency range. This study shows that VLBI is ideal for targeted SETI, including follow-up observations. The techniques presented are equally applicable to next-generation interferometers, such as the long baselines of the Square Kilometre Array (SKA).

preprint2012arXiv

The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies

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

preprint2011arXiv

A VLBA search for binary black holes in active galactic nuclei with double-peaked optical emission line spectra

We have examined a subset of 11 active galactic nuclei (AGN) drawn from a sample of 87 objects that possess double-peaked optical emission line spectra, as put forward by Wang et al. (2009a) and are detectable in the FIRST survey at radio wavelengths. The double-peaked nature of the optical emission line spectra has been suggested as evidence for the existence of binary black holes in these AGN, although this interpretation is controversial. We make a simple suggestion, that direct evidence of binary black holes in these objects could be searched for in the form of dual sources of compact radio emission associated with the AGN. To explore this idea, we have used the Very Long Baseline Array to observe these 11 objects from the Wang et al. (2009a) sample. Of the 11 objects, we detect compact radio emission from two, SDSS J151709+335324 and SDSS J160024+264035. Both objects show single components of compact radio emission. The morphology of SDSS J151709+335324 is consistent with a recent comprehensive multi-wavelength study of this object by Rosario et al. (2010). Assuming that the entire sample consists of binary black holes, we would expect of order one double radio core to be detected, based on radio wavelength detection rates from FIRST and VLBI surveys. We have not detected any double cores, thus this work does not substantially support the idea that AGN with double-peaked optical emission lines contain binary black holes. However, the study of larger samples should be undertaken to provide a more secure statistical result, given the estimated detection rates.

preprint2011arXiv

Dual-frequency VLBI study of Centaurus A on sub-parsec scales

Centaurus A is the closest active galactic nucleus. High resolution imaging using Very Long Baseline Interferometry (VLBI) enables us to study the spectral and kinematic behavior of the radio jet-counterjet system on sub-parsec scales, providing essential information for jet emission and formation models. Our aim is to study the structure and spectral shape of the emission from the central-parsec region of Cen A. As a target of the Southern Hemisphere VLBI monitoring program TANAMI (Tracking Active Galactic Nuclei with Milliarcsecond Interferometry), VLBI observations of Cen A are made regularly at 8.4 and 22.3 GHz with the Australian Long Baseline Array (LBA) and associated telescopes in Antarctica, Chile, and South Africa. The first dual-frequency images of this source are presented along with the resulting spectral index map. An angular resolution of 0.4 mas x 0.7 mas is achieved at 8.4 GHz, corresponding to a linear scale of less than 0.013 pc. Hence, we obtain the highest resolution VLBI image of Cen A, comparable to previous space-VLBI observations. By combining with the 22.3 GHz image, which has been taken without contributing transoceanic baselines at somewhat lower resolution, we present the corresponding dual-frequency spectral index distribution along the sub-parsec scale jet revealing the putative emission regions for recently detected gamma-rays from the core region by Fermi/LAT. We resolve the innermost structure of the milliarcsecond scale jet and counterjet system of Cen A into discrete components. The simultaneous observations at two frequencies provide the highest resolved spectral index map of an AGN jet allowing us to identify multiple possible sites as the origin of the high energy emission.

preprint2011arXiv

e-VLBI observations of Circinus X-1: monitoring of the quiescent and flaring radio emission on AU scales

A recent detection of the peculiar neutron star X-ray binary Circinus X-1 with electronic very long baseline interferometry (e-VLBI) prompted the suggestion that compact, non-variable radio emission persists through the entire 16.6-day orbit of the binary system. We present the results of a high angular resolution monitoring campaign conducted with the Australian Long Baseline Array in real-time e-VLBI mode. e-VLBI observations of Circinus X-1 were made on alternate days over a period of 20 days covering the full binary orbit. A compact radio source associated with Circinus X-1 was clearly detected at orbital phases following periastron passage but no compact radio emission was detected at any other orbital phase, ruling out the presence of a persistent, compact emitting region at our sensitivity levels. The jet was not resolved at any epoch of our 1.4-GHz monitoring campaign, suggesting that the ultrarelativistic flow previously inferred to exist in this source is likely to be dark. We discuss these findings within the context of previous radio monitoring of Circinus X-1.

preprint2011arXiv

Prospects for accurate distance measurements of pulsars with the SKA: enabling fundamental physics

Parallax measurements of pulsars allow for accurate measurements of the interstellar electron density and contribute to accurate tests of general relativity using binary systems. The Square Kilometre Array (SKA) will be an ideal instrument for measuring the parallax of pulsars, because it has a very high sensitivity, as well as baselines extending up to several thousands of kilometres. We performed simulations to estimate the number of pulsars for which the parallax can be measured with the SKA and the distance to which a parallax can be measured. We compare two different methods. The first method measures the parallax directly by utilising the long baselines of the SKA to form high angular resolution images. The second method uses the arrival times of the radio signals of pulsars to fit a transformation between time coordinates in the terrestrial frame and the comoving pulsar frame directly yielding the parallax. We find that with the first method a parallax with an accuracy of 20% or less can be measured up to a maximum distance of 13 kpc, which would include 9,000 pulsars. By timing pulsars with the most stable arrival times for the radio emission, parallaxes can be measured for about 3,600 millisecond pulsars up to a distance of 9 kpc with an accuracy of 20%.

preprint2011arXiv

The first resolved imaging of milliarcsecond-scale jets in Circinus X-1

We present the first resolved imaging of the milliarcsecond-scale jets in the neutron star X-ray binary Circinus X-1, made using the Australian Long Baseline Array. The angular extent of the resolved jets is ~20 milliarcseconds, corresponding to a physical scale of ~150 au at the assumed distance of 7.8 kpc. The jet position angle is relatively consistent with previous arcsecond-scale imaging with the Australia Telescope Compact Array. The radio emission is symmetric about the peak, and is unresolved along the minor axis, constraining the opening angle to be less than 20 degrees. We observe evidence for outward motion of the components between the two halves of the observation. Constraints on the proper motion of the radio-emitting components suggest that they are only mildly relativistic, although we cannot definitively rule out the presence of the unseen, ultra-relativistic (Lorentz factor >15) flow previously inferred to exist in this system.

preprint2010arXiv

Evolution of the pc-scale structure of PKS 1934-638 revisited: first science with the ASKAP and New Zealand telescopes

We have studied the archetypal Gigahertz Peaked Spectrum radio galaxy, PKS 1934-638, using the Australian Long Baseline Array, augmented with two new telescopes that greatly improve the angular resolution of the array. These VLBI observations represent the first scientific results from a new antenna in NZ and the first antenna of the Australian SKA Pathfinder (ASKAP). A compact double radio source, PKS 1934-638, has been monitored over a period of 40 years, and the observation described here provides the latest datum, eight years after the previous observation, to aid in the study of the long-term evolution of the source structure. We take advantage of these new long baselines to probe PKS 1934-638 at the relatively low frequency of 1.4 GHz, in order to examine the effects of optical depth on the structure of the radio source. Optical depth effects, resulting in the observation of frequency dependent structure, may have previously been interpreted in terms of an expansion of the source as a function of time. Expansion and frequency dependent effects are important to disentangle in order to estimate the age of PKS 1934-638. We show that frequency dependent structure effects are likely to be important in PKS 1934-638 and present a simple two-dimensional synchrotron source model in which opacity effects due to synchrotron self-absorption are taken into account. Evidence for expansion of the radio source over 40 years is therefore weak, with consequences for the estimated age of the radio source.

preprint2010arXiv

Interferometric imaging with the 32 element Murchison Wide-field Array

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

preprint2010arXiv

TANAMI - Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry

We present a summary of the observation strategy of TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry), a monitoring program to study the parsec-scale structure and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Australian Long Baseline Array (LBA) and the trans-oceanic antennas Hartebeesthoek, TIGO, and O'Higgins. TANAMI is focusing on extragalactic sources south of -30 degrees declination with observations at 8.4 GHz and 22 GHz every ~2 months at milliarcsecond resolution. The initial TANAMI sample of 43 sources has been defined before the launch of the Fermi Gamma Ray Space Telescope to include the most promising candidates for bright gamma-ray emission to be detected with its Large Area Telescope (LAT). Since November 2008, we have been adding new sources to the sample, which now includes all known radio- and gamma-ray bright AGN of the Southern Hemisphere. The combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN and the nature of relativistic jets.

preprint2010arXiv

TANAMI: Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry I. First-Epoch 8.4 GHz Images

We introduce the TANAMI program (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) which is monitoring an initial sample of 43 extragalactic jets located south of -30 degrees declination at 8.4 GHz and 22 GHz since 2007. All aspects of the program are discussed. First epoch results at 8.4 GHz are presented along with physical parameters derived therefrom. We present first epoch images for 43 sources, some observed for the first time at milliarcsecond resolution. Parameters of these images as well as physical parameters derived from them are also presented and discussed. These and subsequent images from the TANAMI survey are available at http://pulsar.sternwarte.uni-erlangen.de/tanami/ We obtain reliable, high dynamic range images of the southern hemisphere AGN. All the quasars and BL Lac objects in the sample have a single-sided radio morphology. Galaxies are either double-sided, single-sided or irregular. About 28% of the TANAMI sample has been detected by LAT during its first three months of operations. Initial analysis suggests that when galaxies are excluded, sources detected by LAT have larger opening angles than those not detected by LAT. Brightness temperatures of LAT detections and non-detections seem to have similar distributions. The redshift distributions of the TANAMI sample and sub-samples are similar to those seen for the bright gamma-ray AGN seen by LAT and EGRET but none of the sources with a redshift above 1.8 have been detected by LAT.

preprint2009arXiv

Parsec-scale shocks in the kiloparsec-scale jet of Centaurus A

High angular resolution VLBI observations of Centaurus A have been undertaken that allow access to a wide field-of-view, encompassing both the well-studied pc-scale jet and the inner part of the kpc-scale jet. The VLBI observations have detected compact regions of synchrotron emission in the kpc-scale jet that coincide with three stationary features identified from previous VLA monitoring observations. Each of these stationary features is associated with strong localised X-ray emission. The VLBI results strengthen arguments made by previous authors suggesting that the stationary features may be the result of stellar objects or gas clouds traversing the jet flow, intercepting the jet and causing strong shocks. The VLBI data show that the most strongly shocked regions in these features are resolved but have extents no larger than a few pc, reducing the required mass of the typical intercepting object by a factor of ~10 relative to previous estimates, making explanations based on high mass loss stars or low density gas clouds more plausible.

preprint2009arXiv

Precision southern hemisphere VLBI pulsar astrometry II: Measurement of seven parallaxes

Accurate measurement of pulsar distances via astrometry using very long baseline interferometry enables the improvement of Galactic electron density distribution models, improving distance estimates for the vast majority of pulsars for which parallax measurements are unavailable. However, pulsars at southern declinations have been under-represented in previous interferometric astrometry campaigns. In order to redress this imbalance, we have conducted a two-year astrometric campaign targeting eight southern pulsars with the Australian Long Baseline Array. The program summarized in this paper has resulted in the measurement of seven new pulsar parallaxes, with success on objects down to a mean flux density of 0.8 mJy at 1600 MHz. Our results highlight the substantial uncertainties that remain when utilizing free electron density models for individual pulsar distances. Until this study, PSR J0630-2834 was believed to convert 16% of its spin-down energy into x-rays, but our measured parallax distance of 332 (+52 -40) pc has revised this value to <1%. In contrast, PSR J0108-1431 was found to be almost a factor of two more distant than previously thought, making its conversion of spin-down energy to x-rays the most efficient known (>1%). The 8.5 second radio pulsar J2144-3933 was found to be closer than previously predicted, making its apparent 1400 MHz radio luminosity the lowest of any known pulsar (20 microJy kpc^2). We have examined the growing population of neutron stars with accurate parallaxes to determine the effect of distance errors on the underlying neutron star velocity distribution, and find that typical distance errors may be biasing the estimated mean pulsar velocity upwards by 5%, and are likely to exaggerate the distribution's high-velocity tail.

preprint2009arXiv

TANAMI: Milliarcsecond Resolution Observations of Extragalactic Gamma-ray Sources

The TANAMI (Tracking AGN with Austral Milliarcsecond Interferometry) and associated programs provide comprehensive radio monitoring of extragalactic gamma-ray sources south of declination -30 degrees. Joint quasi-simultaneous observations between the Fermi Gamma-ray Space Telescope and ground based observatories allow us to discriminate between competing theoretical blazar emission models. High resolution VLBI observations are the only way to spatially resolve the sub-parsec level emission regions where the high-energy radiation originates. The gap from radio to gamma-ray energies is spanned with near simultaneous data from the Swift satellite and ground based optical observatories. We present early results from the TANAMI program in the context of this panchromatic suite of observations.

preprint2009arXiv

The TANAMI Program

TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) is a monitoring program to study the parsec-scale structures and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Long Baseline Array and associated telescopes. Extragalactic jets south of -30 degrees declination are observed at 8.4 GHz and 22 GHz every two months at milliarcsecond resolution. The initial TANAMI sample is a hybrid radio and gamma-ray selected sample since the combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN.

preprint2008arXiv

The Sub-parsec Scale Radio Properties of Southern Starburst Galaxies. II. Supernova Remnants, the Supernova Rate, and the Ionised Medium in the NGC 4945 Starburst

Wide-field VLBI observations of the nearby starburst galaxy NGC 4945, obtained with the Australian LBA, have produced 2.3 GHz images over two epochs with a maximum spatial resolution of 0.3 pc. 15 sources were detected, 13 of which correspond to sources identified in higher frequency (3 cm and 12 mm) ATCA images. Four of the sources are resolved into shell-like structures ranging between 1.1 to 2.1 pc in diameter. From these data the spectra of 13 compact radio sources in NGC 4945 were modelled; nine were found to be consistent with free-free absorbed power laws and four with a simple power law spectrum. The free-free opacity is highest toward the nucleus but varies significantly throughout the nuclear region, implying that the overall structure of the ionised medium is clumpy. Of the 13 sources, 10 have steep intrinsic spectra associated with synchrotron emission from supernova remnants, the remaining sources have flat intrinsic spectra which may be associated with thermal radio emission. A non-thermal source with a jet-like morphology is detected ~1" from the assumed location of the AGN. A type II supernova (SN) rate upper limit of 15.3/yr is determined for the inner 250 pc region of the galaxy at the 95% confidence level, based on the lack of detection of new sources in observations spanning 1.9 years and a simple model for the evolution of supernova remnants. A type II SN rate of >0.1(v/1e4)/yr is implied from estimates of supernova remnant source counts, sizes and expansion rates, where v is the radial expansion velocity of the supernova remnant in km/s. A SFR of 2.4(v/1e4)<SFR(M>=5Msun)<370 Msun/yr has been estimated directly from the SN rate limits and is of the same order of magnitude as rates determined from integrated FIR and radio luminosities.

preprint2001arXiv

The star-formation rate in the host of GRB 990712

We have observed the host galaxy of GRB 990712 at 1.4 GHz with the Australia Telescope Compact Array, to obtain an estimate of its total star-formation rate. We do not detect a source at the position of the host. The 2 sigma upper limit of 70 microJy implies that the total star-formation rate is lower than 100 Msun/yr, using conservative values for the spectral index and cosmological parameters. This upper limit is in stark contrast with recent reports of radio/submillimeter-determined star-formation rates of roughly 500 Msun/yr for two other GRB host galaxies. Our observations present the deepest radio-determined star-formation rate limit on a GRB host galaxy yet, and show that also from the unobscured radio point-of-view, not every GRB host galaxy is a vigorous starburst.