Source author record

M. Johnston-Hollitt

M. Johnston-Hollitt 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

84works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

84 published item(s)

preprint2022arXiv

Radio fossils, relics, and haloes in Abell 3266: cluster archaeology with ASKAP-EMU and the ATCA

Abell 3266 is a massive and complex merging galaxy cluster that exhibits significant substructure. We present new, highly sensitive radio continuum observations of Abell 3266 performed with the Australian Square Kilometre Array Pathfinder (0.8$-$1.1 GHz) and the Australia Telescope Compact Array (1.1$-$3.1 GHz). These deep observations provide new insights into recently-reported diffuse non-thermal phenomena associated with the intracluster medium, including a 'wrong-way' relic, a fossil plasma source, and an as-yet unclassified central diffuse ridge, which we reveal comprises the brightest part of a large-scale radio halo detected here for the first time. The 'wrong-way' relic is highly atypical of its kind: it exhibits many classical signatures of a shock-related radio relic, while at the same time exhibiting strong spectral steepening. While radio relics are generally consistent with a quasi-stationary shock scenario, the 'wrong-way' relic is not. We study the spectral properties of the fossil plasma source; it exhibits an ultra-steep and highly curved radio spectrum, indicating an extremely aged electron population. The larger-scale radio halo fills much of the cluster centre, and presents a strong connection between the thermal and non-thermal components of the intracluster medium, along with evidence of substructure. Whether the central diffuse ridge is simply a brighter component of the halo, or a mini-halo, remains an open question. Finally, we study the morphological and spectral properties of the multiple complex radio galaxies in this cluster in unprecedented detail, tracing their evolutionary history.

preprint2022arXiv

The merging galaxy cluster Abell 3266 at low radio frequencies

We present new low-frequency ($ν= 88$-216 MHz) observations of the complex merging galaxy cluster Abell 3266. These new observations are taken with the Murchison Widefield Array (MWA) in its Phase II 'extended', long-baseline configuration, offering the highest-resolution low-frequency view of the cluster to date. We report on the detection of four steep spectrum ($α\lesssim -1$ for $S_ν\propto ν^α$) extended radio sources within the cluster. We confirm the detection of a $\sim 570$ kpc radio relic to the south of the cluster, and a possible bridge of emission connecting the relic to the cluster core. We also detect two new ultra-steep-spectrum ($α\lesssim -1.7$) fossil plasma sources to the north and west of the cluster centre without associated compact radio emission. A previously detected radio galaxy in the cluster is also found to have a spectrally steepening tail with steep-spectrum components highlighted by the MWA. We do not detect a giant radio halo in the cluster. After simulating a range of radio haloes at 216 MHz we place upper limits on the radio luminosity corresponding to $\sim 7.2 \times 10^{24}$ W Hz$^{-1}$ at 150 MHz assuming the expected 500 kpc radius, up to a factor 5 lower than expected from scaling relations. Why a giant radio halo is undetected in Abell 3266 is unclear - the timeline of the merger and overall mass of the system as determined by optical and X-ray studies suggest the foundation for a hosting system.

preprint2022arXiv

Wide-Band Spectral Variability of Peaked Spectrum Sources

Characterising spectral variability of radio sources is a technique that offers the ability to determine the astrophysics of the intervening media, source structure, emission and absorption processes. We present broadband (0.072--10 GHz) spectral variability of 15 peaked-spectrum (PS) sources with the Australia Telescope Compact Array (ATCA) and the Murchison Widefield Array (MWA). These 15 PS sources were observed quasi-contemporaneously with ATCA and the MWA four to six times during 2020 with approximately a monthly cadence. Variability was not detected at 1--10GHz frequencies but 13 of the 15 targets show significant variability with the MWA at megahertz frequencies. We conclude the majority of variability seen at megahertz frequencies is due to refractive interstellar scintillation of a compact component ~25 mas across. We also identify four PS sources that show a change in their spectral shape at megahertz frequencies. Three of these sources are consistent with a variable optical depth from an inhomogeneous free-free absorbing cloud around the source. One PS source with a variable spectral shape at megahertz frequencies is consistent with an ejection travelling along the jet. We present spectral variability as a method for determining the physical origins of observed variability and for providing further evidence to support absorption models for PS sources where spectral modelling alone is insufficient.

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.

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

Discovery of a giant radio fossil in the Ophiuchus galaxy cluster

The Ophiuchus galaxy cluster exhibits a curious concave gas density discontinuity at the edge of its cool core. It was discovered in the Chandra X-ray image by Werner and collaborators, who considered a possibility of it being a boundary of an AGN-inflated bubble located outside the core, but discounted this possibility because it required much too powerful an AGN outburst. Using low-frequency (72-240 MHz) radio data from MWA GLEAM and GMRT, we found that the X-ray structure is, in fact, a giant cavity in the X-ray gas filled with diffuse radio emission with an extraordinarily steep radio spectrum. It thus appears to be a very aged fossil of the most powerful AGN outburst seen in any galaxy cluster ($pV\sim 5\times 10^{61}$ erg for this cavity). There is no apparent diametrically opposite counterpart either in X-ray or in the radio. It may have aged out of the observable radio band because of the cluster asymmetry. At present, the central AGN exhibits only a weak radio source, so it should have been much more powerful in the past to have produced such a bubble. The AGN is currently starved of accreting cool gas because the gas density peak is displaced by core sloshing. The sloshing itself could have been set off by this extraordinary explosion if it had occurred in an asymmetric gas core. This dinosaur may be an early example of a new class of sources to be uncovered by low-frequency surveys of galaxy clusters.

preprint2020arXiv

Radio observations of the merging galaxy cluster system Abell 3391-Abell 3395

The pre-merging system of galaxy clusters Abell 3391-Abell 3395 located at a mean redshift of 0.053 has been observed at 1 GHz in an ASKAP/EMU Early Science observation as well as in X-rays with eROSITA. The projected separation of the X-ray peaks of the two clusters is $\sim$50$'$ or $\sim$ 3.1 Mpc. Here we present an inventory of interesting radio sources in this field around this cluster merger. While the eROSITA observations provide clear indications of a bridge of thermal gas between the clusters, neither ASKAP nor MWA observations show any diffuse radio emission coinciding with the X-ray bridge. We derive an upper limit on the radio emissivity in the bridge region of $\langle J \rangle_{1\,{\rm GHz}}< 1.2 \times 10^{-44} {\rm W}\, {\rm Hz}^{-1} {\rm m}^{-3}$. A non-detection of diffuse radio emission in the X-ray bridge between these two clusters has implications for particle-acceleration mechanisms in cosmological large-scale structure. We also report extended or otherwise noteworthy radio sources in the 30 deg$^2$ field around Abell 3391-Abell 3395. We identified 20 Giant Radio Galaxies, plus 7 candidates, with linear projected sizes greater than 1 Mpc. The sky density of field radio galaxies with largest linear sizes of $>0.7$ Mpc is $\approx 1.7$ deg$^{-2}$, three times higher than previously reported. We find no evidence for a cosmological evolution of the population of Giant Radio Galaxies. Moreover, we find seven candidates for cluster radio relics and radio halos.

preprint2020arXiv

SPT-CL J2032-5627: a new Southern double relic cluster observed with ASKAP

We present a radio and X-ray analysis of the galaxy cluster SPT-CL J2032-5627. Investigation of public data from the Australian Square Kilometre Array Pathfinder (ASKAP) at 943 MHz shows two previously undetected radio relics at either side of the cluster. For both relic sources we utilise archival Australia Telescope Compact Array (ATCA) data at 5.5 GHz in conjunction with the new ASKAP data to determine that both have steep integrated radio spectra ($α_\mathrm{SE} = -1.52 \pm 0.10$ and $α_\mathrm{NW,full} = -1.18 \pm 0.10$ for the southeast and northwest relic sources, respectively). No shock is seen in XMM-Newton observations, however, the southeast relic is preceded by a cold front in the X-ray emitting intra-cluster medium. We suggest the lack of a detectable shock may be due to instrumental limitations, comparing the situation to the southeast relic in Abell 3667. We compare the relics to the population of double relic sources and find they are located below the current power-mass ($P$-$M$) scaling relation. We present an analysis of the low-surface brightness sensitivity of ASKAP and the ATCA, the excellent sensitivity of both allow the ability to find heretofore undetected diffuse sources, suggesting these low-power radio relics will become more prevalent in upcoming large-area radio surveys such as the Evolutionary Map of the Universe (EMU).

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 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 KAT-7 view of a low-mass sample of galaxy clusters

Radio observations over the last two decades have provided evidence that diffuse synchrotron emission in the form of megaparsec-scale radio halos in galaxy clusters is likely tracing regions of the intracluster medium where relativistic particles are accelerated during cluster mergers. In this paper we present results of a survey of 14 galaxy clusters carried out with the 7-element Karoo Array Telescope at 1.86 GHz, aimed to extend the current studies of radio halo occurrence to systems with lower masses (M$_{\rm 500} > 4\times10^{14}$ M${_\odot}$). We found upper limits at the $0.6 - 1.9 \times 10^{24}$ Watt Hz$^{-1}$ level for $\sim 50\%$ of the sample, confirming that bright radio halos in less massive galaxy clusters are statistically rare.

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

An Optical Analysis of the Merging Cluster Abell 3888

In this paper we present new AAOmega spectroscopy of 254 galaxies within a 30' radius around Abell 3888. We combine these data with the existing redshifts measured in a one degree radius around the cluster and performed a substructure analysis. We confirm 71 member galaxies within the core of A3888 and determine a new average redshift and velocity dispersion for the cluster of 0.1535 +\- 0.0009 and 1181 +\- 197 km/s, respectively. The cluster is elongated along an East-West axis and we find the core is bimodal along this axis with two sub-groups of 26 and 41 members detected. Our results suggest that A3888 is a merging system putting to rest the previous conjecture about the morphological status of the cluster derived from X-ray observations. In addition to the results on A3888 we also present six newly detected galaxy over-densities in the field, three of which we classify as new galaxy clusters.

preprint2016arXiv

An Overview of the SKA Science Analysis Pipeline

When completed the Square Kilometre Array (SKA) will feature an unprecedented rate of image generation. While previous generations of telescopes have relied on human expertise to extract scientifically interesting information from the images, the sheer data volume of the data will now make this impractical. Additionally, the rate at which data are accrued will not allow traditional imaging products to be stored indefinitely for later inspection meaning there is a strong imperative to discard uninteresting data in pseudo-real time. Here we outline components of the SKA science analysis pipeline being developed to produce a series of data products including continuum images, spectral cubes and Faraday depth spectra. We discuss a scheme to automatically extract value from these products and discard scientifically uninteresting data. This pipeline is thus expected to give both an increase in scientific productivity, and offers the possibility of reduced data archive size producing a considerable saving.

preprint2016arXiv

Another shock for the Bullet cluster, and the source of seed electrons for radio relics

With Australia Telescope Compact Array observations, we detect a highly elongated Mpc-scale diffuse radio source on the eastern periphery of the Bullet cluster 1E0657-55.8, which we argue has the positional, spectral and polarimetric characteristics of a radio relic. This powerful relic (2.3+/-0.1 x 10^25 W Hz^-1) consists of a bright northern bulb and a faint linear tail. The bulb emits 94% of the observed radio flux and has the highest surface brightness of any known relic. Exactly coincident with the linear tail we find a sharp X-ray surface brightness edge in the deep Chandra image of the cluster -- a signature of a shock front in the hot intracluster medium (ICM), located on the opposite side of the cluster to the famous bow shock. This new example of an X-ray shock coincident with a relic further supports the hypothesis that shocks in the outer regions of clusters can form relics via diffusive shock (re-)acceleration. Intriguingly, our new relic suggests that seed electrons for reacceleration are coming from a local remnant of a radio galaxy, which we are lucky to catch before its complete disruption. If this scenario, in which a relic forms when a shock crosses a well-defined region of the ICM polluted with aged relativistic plasma -- as opposed to the usual assumption that seeds are uniformly mixed in the ICM -- is also the case for other relics, this may explain a number of peculiar properties of peripheral relics.

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

CHIPS: The Cosmological HI Power Spectrum Estimator

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

preprint2016arXiv

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

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

preprint2016arXiv

Detecting Diffuse Sources in Astronomical Images

We present an algorithm capable of detecting diffuse, dim sources of any size in an astronomical image. These sources often defeat traditional methods for source finding, which expand regions around points of high intensity. Extended sources often have no bright points and are only detectable when viewed as a whole, so a more sophisticated approach is required. Our algorithm operates at all scales simultaneously by considering a tree of nested candidate bounding boxes, and inverts a hierarchical Bayesian generative model to obtain the probability of sources existing at given locations and sizes. This model naturally accommodates the detection of nested sources, and no prior knowledge of the distribution of a source, or even the background, is required. The algorithm scales nearly linear with the number of pixels making it feasible to run on large images, and requires minimal parameter tweaking to be effective. We demonstrate the algorithm on several types of astronomical and artificial images.

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 of noise and confusion in the MWA GLEAM year 1 survey

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a new relatively low resolution, contiguous 72-231 MHz survey of the entire sky south of declination +25 deg. In this paper, we outline one approach to determine the relative contribution of system noise, classical confusion and sidelobe confusion in GLEAM images. An understanding of the noise and confusion properties of GLEAM is essential if we are to fully exploit GLEAM data and improve the design of future low-frequency surveys. Our early results indicate that sidelobe confusion dominates over the entire frequency range, implying that enhancements in data processing have the potential to further reduce the noise.

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

Low Radio Frequency Observations and Spectral Modelling of the Remnant of Supernova 1987A

We present Murchison Widefield Array observations of the supernova remnant (SNR) 1987A between 72 and 230 MHz, representing the lowest frequency observations of the source to date. This large lever arm in frequency space constrains the properties of the circumstellar medium created by the progenitor of SNR 1987A when it was in its red supergiant phase. As of late-2013, the radio spectrum of SNR 1987A between 72 MHz and 8.64 GHz does not show any deviation from a non-thermal power-law with a spectral index of $-0.74 \pm 0.02$. This spectral index is consistent with that derived at higher frequencies, beneath 100 GHz, and with a shock in its adiabatic phase. A spectral turnover due to free-free absorption by the circumstellar medium has to occur below 72 MHz, which places upper limits on the optical depth of $\leq$ 0.1 at a reference frequency of 72 MHz, emission measure of $\lesssim$ 13,000 cm$^{-6}$ pc, and an electron density of $\lesssim$ 110 cm$^{-3}$. This upper limit on the electron density is consistent with the detection of prompt radio emission and models of the X-ray emission from the supernova. The electron density upper limit implies that some hydrodynamic simulations derived a red supergiant mass loss rate that is too high, or a wind velocity that is too low. The mass loss rate of $\sim 5 \times 10^{-6}$ $M_{\odot}$ yr$^{-1}$ and wind velocity of 10 km s$^{-1}$ obtained from optical observations are consistent with our upper limits, predicting a current turnover frequency due to free-free absorption between 5 and 60 MHz.

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

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

Three Tools to Aid Visualisation of FITS Files for Astronomy

Increasingly there is a need to develop astronomical visualisation and manipulations tools which allow viewers to interact with displayed data directly, in real time and across a range of platforms. In addition, increases in dynamic range available for astronomical images with next generation telescopes have led to a desire to develop enhanced visualisations capable of presenting information across a wide range of intensities. This paper describes three new tools for astronomical visualisation and image manipulation that are the result of a collaboration between software engineers and radio astronomers. The first tool, FITS3D, is a fast, interactive 3D data cube viewer designed to allow real-time interactive comparisons of multiple spectral line data cubes simultaneously. It features region specific selection manipulation including smoothing. The second tool, FITS2D, aids the visualisation and manipulation of 2D fits images. The tool supports the interactive creation of free-form masks which allow the user to extract any (potentially non-contiguous) subset of a fits image. It also supports annotations which can be placed without affecting the underlying data. The final tool is an R package for applying high dynamic range compression techniques to 2D fits images. This allows the full range of pixel brightness to be imaged in a single image, simultaneously showing the detail in bright sources while preserving the distinction of faint sources. Here we will present these three tools and demonstrate their capability using images from a range of astronomical images.

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

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

Cluster magnetic fields through the study of polarized radio halos in the SKA era

Galaxy clusters are unique laboratories to investigate turbulent fluid motions and large scale magnetic fields. Synchrotron radio halos at the center of merging galaxy clusters provide the most spectacular and direct evidence of the presence of relativistic particles and magnetic fields associated with the intracluster medium. The study of polarized emission from radio halos is extremely important to constrain the properties of intracluster magnetic fields and the physics of the acceleration and transport of the relativistic particles. However, detecting this polarized signal is a very hard task with the current radio facilities.We use cosmological magneto-hydrodynamical simulations to predict the expected polarized surface brightness of radio halos at 1.4 GHz. We compare these expectations with the sensitivity and the resolution reachable with the SKA1. This allows us to evaluate the potential for studying intracluster magnetic fields in the surveys planned for SKA1.

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

Frontier Fields Clusters: Chandra and JVLA View of the Pre-Merging Cluster MACS J0416.1-2403

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.

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

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

Tailed Radio Galaxies as Probes of Cluster Physics in the Square Kilometre Array Era

In recent years, the use of tailed radio galaxies as environmental probes has gained momentum as a method for galaxy cluster detection, examining the dynamics of individual clusters, measuring the density and velocity flows in the intra-cluster medium, and for probing cluster magnetic fields. To date instrumental limitations in terms of resolution and sensitivity have confined this research to the local (z < 0.7) Universe. The advent of SKA-1 surveys however will allow detection of well over 1 million tailed radio galaxies and their associated galaxy clusters out to redshifts of 2 or more. This is in fact ten times more than the current number of known clusters in the Universe. Such a substantial sample of tailed galaxies will provide an invaluable tool not only for detecting clusters, but also for characterizing their intra-cluster medium, magnetic fields and dynamical state as a function of cosmic time. In this paper we present an analysis of the usability of tailed radio galaxies as tracers of dense environments extrapolated from existing deep radio surveys such the Extended Chandra Deep Field-South.

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 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.

preprint2015arXiv

Unravelling the origin of large-scale magnetic fields in galaxy clusters and beyond through Faraday Rotation Measures with the SKA

We investigate the possibility for the SKA to detect and study the magnetic fields in galaxy clusters and in the less dense environments surrounding them using Faraday Rotation Measures. To this end, we produce 3-dimensional magnetic field models for galaxy clusters of different masses and in different stages of their evolution, and derive mock rotation measure observations of background radiogalaxies. According to our results, already in phase I, we will be able to infer the magnetic field properties in galaxy clusters as a function of the cluster mass, down to $10^{13}$ solar-masses. Moreover, using cosmological simulations to model the gas density, we have computed the expected rotation measure through shock-fronts that occur in the intra-cluster medium during cluster mergers. The enhancement in the rotation measure due to the density jump will permit to constraint the magnetic field strength and structure after the shock passage. SKA observations of polarised sources located behind galaxy clusters will answer several questions about the magnetic field strength and structure in galaxy clusters, and its evolution with cosmic time.

preprint2015arXiv

Using rotation measure grids to detect cosmological magnetic fields -- a Bayesian approach

Determining magnetic field properties in different environments of the cosmic large-scale structure as well as their evolution over redshift is a fundamental step toward uncovering the origin of cosmic magnetic fields. Radio observations permit the study of extragalactic magnetic fields via measurements of the Faraday depth of extragalactic radio sources. Our aim is to investigate how much different extragalactic environments contribute to the Faraday depth variance of these sources. We develop a Bayesian algorithm to distinguish statistically Faraday depth variance contributions intrinsic to the source from those due to the medium between the source and the observer. In our algorithm the Galactic foreground and the measurement noise are taken into account as the uncertainty correlations of the galactic model. Additionally, our algorithm allows for the investigation of possible redshift evolution of the extragalactic contribution. This work presents the derivation of the algorithm and tests performed on mock observations. With cosmic magnetism being one of the key science projects of the new generation of radio interferometers we have made predictions for the algorithm's performance on data from the next generation of radio interferometers. Applications to real data are left for future work.

preprint2015arXiv

Using Tailed Radio Galaxies to Probe the Environment and Magnetic Field of Galaxy Clusters in the SKA Era

The morphology of tailed radio galaxies is an invaluable source of environmental information, in which a history of the past interactions in the intra-cluster medium, such as complex galaxy motions and cluster merger shocks, are preserved. In recent years, the use of tailed radio galaxies as environmental probes has gained momentum as a method for galaxy cluster detection, examining the dynamics of individual clusters, measuring the density and velocity flows in the intra-cluster medium, and for probing cluster magnetic fields. To date instrumental limitations in terms of resolution and sensitivity have confined this research to the local (z < 0.7) Universe. The advent of SKA1 surveys however will allow detection of roughly 1,000,000 tailed radio galaxies and their associated galaxy clusters out to redshifts of 2 or more. This is in fact ten times more than the current number of known clusters in the Universe. Additionally between 50,000 and 100,000 tailed radio galaxies will be sufficiently polarized to allow characterization of the magnetic field of their parent cluster. Such a substantial sample of tailed galaxies will provide an invaluable tool not only for detecting clusters, but also for characterizing their intra-cluster medium, magnetic fields and dynamical state as a function of cosmic time. In this chapter we present an analysis of the usability of tailed radio galaxies as tracers of dense environments extrapolated from existing deep radio surveys.

preprint2014arXiv

Cluster Radio Halos at the crossroads between astrophysics and cosmology in the SKA era

Giant Radio Halos (RH) are diffuse, Mpc-sized, synchrotron radio sources observed in a fraction of merging galaxy clusters. The current scenario for the origin of RHs assumes that turbulence generated during cluster mergers re-accelerates pre-existing fossil and/or secondary electrons in the intra-cluster-medium (ICM) to the energies necessary to produce the observed radio emission. Moreover, more relaxed clusters could host diffuse "off state" halos produced by secondary electrons. In this Chapter we use Monte Carlo simulations, that combine turbulent-acceleration physics and the generation of secondaries in the ICM, to calculate the occurrence of RHs in the Universe, their spectral properties and connection with properties of the hosting clusters. Predictions for SKA1 surveys are presented at low (100-300 MHz) and mid (1-2 GHz) frequencies assuming the expected sensitivities and spatial resolutions of SKA1. SKA1 will step into an unexplored territory allowing us to study the formation and evolution of RHs in a totally new range of cluster masses and redshift, allowing firm tests of the current theoretical hypothesis. In particular, the combination of SKA1-LOW and SUR will allow the discovery of ~1000 ultrasteep- spectrum halos and to detect for the very first time "off state" RHs. We expect that at least ~2500 giant RHs will be discovered by SKA1-LOW surveys up to z~0.6. Remarkably these surveys will be sensitive to RHs in a cluster mass range (down to ~10^14 solar masses) and redshifts (up to ~1) that are unexplored by current observations. SKA1 surveys will be highly competitive with present and future SZ-surveys in the detection of high-redshift massive objects.

preprint2014arXiv

Discovery of large-scale diffuse radio emission and of a new galaxy cluster in the surroundings of MACSJ0520.7-1328

We report the discovery of large-scale diffuse radio emission South-East of the galaxy cluster MACS J0520.7-1328, detected through high sensitivity Giant Metrewave Radio Telescope 323 MHz observations. This emission is dominated by an elongated diffuse radio source and surrounded by other features of lower surface brightness. Patches of these faint sources are marginally detected in a 1.4 GHz image obtained through a re-analysis of archival NVSS data. Interestingly, the elongated radio source coincides with a previously unclassified extended X-ray source. We perform a multi-wavelength analysis based on archival infrared, optical and X-ray Chandra data. We find that this source is a low-temperature (~3.6 keV) cluster of galaxies, with indications of a disturbed dynamical state, located at a redshift that is consistent with the one of the main galaxy cluster MACS J0520.7-132 (z=0.336). We suggest that the diffuse radio emission is associated with the non-thermal components in the intracluster and intergalactic medium in and around the newly detected cluster. We are planning deeper multi-wavelength and multi-frequency radio observations to accurately investigate the dynamical scenario of the two clusters and to address more precisely the nature of the complex radio emission.

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

Multiwavelength Observations of AB Doradus

We have observed the bright, magnetically active multiple star AB Doradus in a multiwavelength campaign centring around two large facility allocations in November 2006 and January, 2007. Our observations have covered at least three large flares. These flares were observed to produce significant hardening of the X-ray spectra during their very initial stages. We monitored flare-related effects using the Suzaku X-ray satellite and the Australia Telescope Compact Array at 3.6 and 6 cm. Observations at 11 and 21 cm were also included, but they were compromised by interference. From our multiwavelength coverage we find that the observed effects can be mainly associated with a large active region near longitude zero. The second major X-ray and microwave flare of Jan 8, 2007 was observed with a favourable geometry that allowed its initial high-energy impulsive phase to be observed in the higher frequency range of Suzaku's XIS detectors. The fractional circular polarisation was measured for the complete runs, for 25 min integrations and, at 4.80 GHz, for 5 min integrations. Most of the full data sets showed circular polarisation fractions from AB Dor B that were significant at greater than the 3 sigma level. In several of the 5 min integrations at 4.80 and 8.64 GHz this fraction reached a significance level between 3 and 9 sigma. Lack of angular resolution prevented identification of these high V/I values with one or other of the two low-mass red-dwarf components of AB Dor B.

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 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.

preprint2014arXiv

Using radio jets of PKS J0334-3900 to probe the intra-cluster medium in A3135

We present a multi-wavelength study of the radio galaxy PKS J0334-3900, which resides at the centre of Abell 3135. Using Australia Telescope Compact Array (ATCA) observations at 1.4, 2.5, 4.6 & 8.6 GHz, we performed a detailed analysis of PKS J0334-3900. The morphology and spectral indices give physical parameters that constrain the dynamical history of the galaxy, which we use to produce a simulation of PKS J0334-3900. This simulation shows that the morphology can be generated by a wind in the intracluster medium (ICM), orbital motion caused by a companion galaxy, and precession of the black hole (BH).

preprint2014arXiv

Using the morphology and magnetic fields of tailed radio galaxies as environmental probes

Bent-tailed (BT) radio sources have long been known to trace over densities in the Universe up to z ~ 1 and there is increasing evidence this association persists out to redshifts of 2. The morphology of the jets in BT galaxies is primarily a function of the environment that they have resided in and so BTs provide invaluable clues as to their local conditions. Thus, not only can samples of BT galaxies be used as signposts of large-scale structure, but are also valuable for obtaining a statistical measurement of properties of the intra-cluster medium including the presence of cluster accretion shocks & winds, and as historical anemometers, preserving the dynamical history of their surroundings in their jets. We discuss the use of BTs to unveil large-scale structure and provide an example in which a BT was used to unlock the dynamical history of its host cluster. In addition to their use as density and dynamical indicators, BTs are useful probes of the magnetic field on their environment on scales which are inaccessible to other methods. Here we discuss a novel way in which a particular sub-class of BTs, the so-called `corkscrew' galaxies might further elucidate the coherence lengths of the magnetic fields in their vicinity. Given that BTs are estimated to make up a large population in next generation surveys we posit that the use of jets in this way could provide a unique source of environmental information for clusters and groups up to z = 2.

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

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.

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

On the reliability of polarization estimation using Rotation Measure Synthesis

We benchmark the reliability of the Rotation Measure (RM) synthesis algorithm using the 1005 Centaurus A field sources of Feain et al. (2009). The RM synthesis solutions are compared with estimates of the polarization parameters using traditional methods. This analysis provides verification of the reliability of RM synthesis estimates. We show that estimates of the polarization parameters can be made at lower S/N if the range of RMs is bounded, but reliable estimates of individual sources with unusual RMs require unconstrainted solutions and higher S/N. We derive from first principles the statistical properties of the polarization amplitude associated with RM synthesis in the presence of noise. The amplitude distribution depends explicitly on the amplitude of the underlying (intrinsic) polarization signal. Hence it is necessary to model the underlying polarization signal distribution in order to estimate the reliability and errors in polarization parameter estimates. We introduce a Bayesian method to derive the distribution of intrinsic amplitudes based on the distribution of measured amplitudes. The theoretically-derived distribution is compared with the empirical data to provide quantitative estimates of the probability that an RM synthesis solution is correct as a function of S/N. We provide quantitative estimates of the probability that any given RM synthesis solution is correct as a function of measured polarized amplitude and the intrinsic polarization amplitude compared to the noise.

preprint2012arXiv

Radio Halos in future surveys in the radio continuum

Giant radio halos (RH) are Mpc-scale synchrotron sources detected in a significant fraction of massive and merging galaxy clusters.Their statistical properties can be used to discriminate among various models for their origin. Theoretical predictions are important as new radio telescopes are about to begin to survey the sky at low and high frequencies with unprecedented sensitivity. We carry out Monte Carlo simulations to model the formation and evolution of RH in a cosmological framework by assuming that RH are either generated in turbulent merging clusters, or are purely hadronic sources generated in more relaxed clusters, "off-state" halos. The models predict that the luminosity function of RH at high radio luminosities is dominated by the contribution of RH generated in turbulent clusters. The generation of these RH becomes less efficient in less massive systems causing a flattening of the luminosity function at lower luminosities. This flattening is compensated by the contribution of "off-state" RH that dominate at lower luminosities. By restricting to clusters at z<0.6, we show that the planned EMU+WODAN surveys at 1.4 GHz have the potential to detect up to ~200 RH, increasing their number by one order of magnitude. A fraction of these sources will be "off-state" RH that should be found at flux level < 10 mJy, presently accessible only to deep pointed observations. We also explore the synergy between the Tier 1 LOFAR survey at 150 MHz and the EMU+WODAN surveys at 1.4 GHz. We predict a larger number of RH in the LOFAR survey due to the high LOFAR sensitivity, but also due to the existence of RH with very steep spectrum that glow up preferentially at lower frequencies. These RH are only predicted in the framework of turbulent re-acceleration models and should not have counterparts in the EMU+WODAN surveys, thus the combination of the two surveys will test theoretical models.

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

EMU: Evolutionary Map of the Universe

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

preprint2011arXiv

Very High Angular Resolution Science with the Square Kilometre Array

Preliminary specifications for the Square Kilometre Array (SKA) call for 25% of the total collecting area of the dish array to be located at distances greater than 180 km from the core, with a maximum baseline of at least 3000 km. The array will provide angular resolution ~ 40 - 2 mas at 0.5 - 10 GHz with image sensitivity reaching < 50 nJy/beam in an 8 hour integration with 500 MHz bandwidth. Given these specifications, the high angular resolution component of the SKA will be capable of detecting brightness temperatures < 200 K with milliarcsecond-scale angular resolution. The aim of this article is to bring together in one place a discussion of the broad range of new and important high angular resolution science that will be enabled by the SKA, and in doing so, address the merits of long baselines as part of the SKA. We highlight the fact that high angular resolution requiring baselines greater than 1000 km provides a rich science case with projects from many areas of astrophysics, including important contributions to key SKA science.

preprint2010arXiv

Current New Zealand Activities in Radio Astronomy: Building Capacity in Engineering & Science for the Square Kilometre Array

We present an update on the NZ-wide advances in the field of Radio Astronomy and Radio Engineering with a particular focus on contributions, not thus reported elsewhere, which hope to either directly or indirectly contribute to New Zealand's engagement with the international Square Kilometre Array (SKA) project. We discuss the status of the SKA project in New Zealand with particular reference to activities of the New Zealand Square Kilometre Array Research and Development Consortium.