Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.

This paper is concerned with algorithms for calibration of direction dependent effects (DDE) in aperture synthesis radio telescopes (ASRT). After correction of Direction Independent Effects (DIE) using self-calibration, imaging performance can be limited by the imprecise knowledge of the forward gain of the elements in the array. In general, the forward gain pattern is directionally dependent and varies with time due to a number of reasons. Some factors, such as rotation of the primary beam with Parallactic Angle for Azimuth-Elevation mount antennas are known a priori. Some, such as antenna pointing errors and structural deformation/projection effects for aperture-array elements cannot be measured {\em a priori}. Thus, in addition to algorithms to correct for DD effects known a priori, algorithms to solve for DD gains are required for high dynamic range imaging. Here, we discuss a mathematical framework for antenna-based DDE calibration algorithms and show that this framework leads to computationally efficient optimal algorithms which scale well in a parallel computing environment. As an example of an antenna-based DD calibration algorithm, we demonstrate the Pointing SelfCal algorithm to solve for the antenna pointing errors. Our analysis show that the sensitivity of modern ASRT is sufficient to solve for antenna pointing errors and other DD effects. We also discuss the use of the Pointing SelfCal algorithm in real-time calibration systems and extensions for antenna Shape SelfCal algorithm for real-time tracking and corrections for pointing offsets and changes in antenna shape.

We describe the performance of the Boolardy Engineering Test Array (BETA), the prototype for the Australian Square Kilometre Array Pathfinder telescope ASKAP. BETA is the first aperture synthesis radio telescope to use phased array feed technology, giving it the ability to electronically form up to nine dual-polarization beams. We report the methods developed for forming and measuring the beams, and the adaptations that have been made to the traditional calibration and imaging procedures in order to allow BETA to function as a multi-beam aperture synthesis telescope. We describe the commissioning of the instrument and present details of BETA's performance: sensitivity, beam characteristics, polarimetric properties and image quality. We summarise the astronomical science that it has produced and draw lessons from operating BETA that will be relevant to the commissioning and operation of the final ASKAP telescope.

The Boolardy Engineering Test Array is a 6 x 12 m dish interferometer and the prototype of the Australian Square Kilometre Array Pathfinder (ASKAP), equipped with the first generation of ASKAP's phased array feed (PAF) receivers. These facilitate rapid wide-area imaging via the deployment of simultaneous multiple beams within a 30 square degree field of view. By cycling the array through 12 interleaved pointing positions and using 9 digitally formed beams we effectively mimic a traditional 1 hour x 108 pointing survey, covering 150 square degrees over 711 - 1015 MHz in 12 hours of observing time. Three such observations were executed over the course of a week. We verify the full bandwidth continuum imaging performance and stability of the system via self-consistency checks and comparisons to existing radio data. The combined three epoch image has arcminute resolution and a 1-sigma thermal noise level of 375 micro-Jy per beam, although the effective noise is a factor 3 higher due to residual sidelobe confusion. From this we derive a catalogue of 3,722 discrete radio components, using the 35 percent fractional bandwidth to measure in-band spectral indices for 1,037 of them. A search for transient events reveals one significantly variable source within the survey area. The survey covers approximately two-thirds of the Spitzer South Pole Telescope Deep Field. This pilot project demonstrates the viability and potential of using PAFs to rapidly and accurately survey the sky at radio wavelengths.

We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The absorption line is detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKS B1740$-$517 and demonstrates ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. Optical spectroscopy and imaging using the Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the radio source. The narrow OIII emission lines show clear double-peaked structure, indicating either large-scale outflow or rotation of the ionized gas. Archival data from the \emph{XMM-Newton} satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the active galactic nucleus. The HI absorption profile is complex, with four distinct components ranging in width from 5 to 300 km s$^{-1}$ and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circumnuclear disc or ring structure aligned with the radio jet, we find evidence for a possible broad outflow of neutral gas moving at a radial velocity of $v \sim 300$ km s$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$ yr) is cocooned within a dense medium and may be driving circumnuclear neutral gas in an outflow of $\sim$ 1 $\mathrm{M}_{\odot}$ yr$^{-1}$.

We present HI imaging of the galaxy group IC 1459 carried out with six antennas of the Australian SKA Pathfinder equipped with phased-array feeds. We detect and resolve HI in eleven galaxies down to a column density of $\sim10^{20}$ cm$^{-2}$ inside a ~6 deg$^2$ field and with a resolution of ~1 arcmin on the sky and ~8 km/s in velocity. We present HI images, velocity fields and integrated spectra of all detections, and highlight the discovery of three HI clouds -- two in the proximity of the galaxy IC 5270 and one close to NGC 7418. Each cloud has an HI mass of $10^9$ M$_\odot$ and accounts for ~15% of the HI associated with its host galaxy. Available images at ultraviolet, optical and infrared wavelengths do not reveal any clear stellar counterpart of any of the clouds, suggesting that they are not gas-rich dwarf neighbours of IC 5270 and NGC 7418. Using Parkes data we find evidence of additional extended, low-column-density HI emission around IC 5270, indicating that the clouds are the tip of the iceberg of a larger system of gas surrounding this galaxy. This result adds to the body of evidence on the presence of intra-group gas within the IC 1459 group. Altogether, the HI found outside galaxies in this group amounts to several times $10^9$ M$_\odot$, at least 10% of the HI contained inside galaxies. This suggests a substantial flow of gas in and out of galaxies during the several billion years of the group's evolution.

This paper describes the system architecture of a newly constructed radio telescope - the Boolardy Engineering Test Array, which is a prototype of the Australian Square Kilometre Array Pathfinder telescope. Phased array feed technology is used to form multiple simultaneous beams per antenna, providing astronomers with unprecedented survey speed. The test array described here is a 6-antenna interferometer, fitted with prototype signal processing hardware capable of forming at least 9 dual-polarisation beams simultaneously, allowing several square degrees to be imaged in a single pointed observation. The main purpose of the test array is to develop beamforming and wide-field calibration methods for use with the full telescope, but it will also be capable of limited early science demonstrations.

We describe a method performing w-projection using the fast Gauss transform of Strain (1991). We derive the theoretical performance, and simulate the actual performance for a range of w for a canonical array. While our implementation is dominated by overheads, we argue that this approach could for the basis of a higher-performing algorithms with particular application to the Square Kilometer Array.

Wide-field radio interferometric telescopes such as the Square Kilometre Array now being designed are subject to a number of aberrations. One particularly pernicious aberration is that due to non-coplanar baselines whereby long baselines incur a quadratic image-plane phase error. There are numerous algorithms for dealing with the non-coplanar baselines effect. As a result of our experience with developing processing software for the Australian Square Kilometre Array Pathfinder, we advocate the use of a hybrid algorithm, called w snapshots, based on a combination of w projection and snapshot imaging. This hybrid overcomes some of the deficiencies of each and has advantages from both. Compared to pure w projection, w snapshots uses less memory and execution time, and compared to pure snapshot imaging, w snapshots uses less memory and is more accurate. At the asymptotes, w snapshots devolves to w projection and to snapshots.

Searching for dispersed radio pulses in interferometric data is of great scientific interest, but poses a formidable computational burden. Here we present two efficient, new antenna-coherent solutions: The Chirpolator and The Chimageator. We describe the equations governing both techniques and propose a number of novel optimizations. We compare the implementation costs of our techniques with classical methods using three criteria: the operations rates (1) before and (2) after the integrate-and-dump stage, and (3) the data rate directly after the integrate-and-dump stage. When compared with classical methods, our techniques excel in the regime of sparse arrays, where they both require substantially lower data rates, and The Chirpolator requires a much lower post-integrator operations rate. In general, our techniques require more pre-integrator operations than the classical ones. We argue that the data and operations rates required by our techniques are better matched to future supercomputer architectures, where the arithmetic capability is outstripping the bandwidth capability. Our techniques are, therefore, viable candidates for deploying on future interferometers such as the Square Kilometer Array.

A 45 deg2 radio continuum imaging campaign of the nearest radio galaxy, Centaurus A, is reported. Using the Australia Telescope Compact Array and the Parkes 64m radio telescope at 1.4 GHz, the spatial resolution of the resultant image is ~600 pc, resolving the ~500 kpc giant radio lobes with approximately five times better physical resolution compared to any previous image, and making this the most detailed radio continuum image of any radio galaxy to date. In this paper, we present these new data and discuss briefly some of the most interesting morphological features that we have discovered in the images... abstract truncated.

Aims : We describe MS-MFS, a multi-scale multi-frequency deconvolution algorithm for wide-band synthesis-imaging, and present imaging results that illustrate the capabilities of the algorithm and the conditions under which it is feasible and gives accurate results. Methods : The MS-MFS algorithm models the wide-band sky-brightness distribution as a linear combination of spatial and spectral basis functions, and performs image-reconstruction by combining a linear-least-squares approach with iterative $\chi^2$ minimization. This method extends and combines the ideas used in the MS-CLEAN and MF-CLEAN algorithms for multi-scale and multi-frequency deconvolution respectively, and can be used in conjunction with existing wide-field imaging algorithms. We also discuss a simpler hybrid of spectral-line and continuum imaging methods and point out situations where it may suffice. Results : We show via simulations and application to multi-frequency VLA data and wideband EVLA data, that it is possible to reconstruct both spatial and spectral structure of compact and extended emission at the continuum sensitivity level and at the angular resolution allowed by the highest sampled frequency.

Faraday rotation measure (RM) synthesis is an important tool to study and analyze galactic and extra-galactic magnetic fields. Since there is a Fourier relation between the Faraday dispersion function and the polarized radio emission, full reconstruction of the dispersion function requires knowledge of the polarized radio emission at both positive and negative square wavelengths $\lambda^2$. However, one can only make observations for $\lambda^2 > 0$. Furthermore observations are possible only for a limited range of wavelengths. Thus reconstructing the Faraday dispersion function from these limited measurements is ill-conditioned. In this paper, we propose three new reconstruction algorithms for RM synthesis based upon compressive sensing/sampling (CS). These algorithms are designed to be appropriate for Faraday thin sources only, thick sources only, and mixed sources respectively. Both visual and numerical results show that the new RM synthesis methods provide superior reconstructions of both magnitude and phase information than RM-CLEAN

Compressive sampling is a new paradigm for sampling, based on sparseness of signals or signal representations. It is much less restrictive than Nyquist-Shannon sampling theory and thus explains and systematises the widespread experience that methods such as the H\"ogbom CLEAN can violate the Nyquist-Shannon sampling requirements. In this paper, a CS-based deconvolution method for extended sources is introduced. This method can reconstruct both point sources and extended sources (using the isotropic undecimated wavelet transform as a basis function for the reconstruction step). We compare this CS-based deconvolution method with two CLEAN-based deconvolution methods: the H\"ogbom CLEAN and the multiscale CLEAN. This new method shows the best performance in deconvolving extended sources for both uniform and natural weighting of the sampled visibilities. Both visual and numerical results of the comparison are provided.

The Parkes pulsar data archive currently provides access to 144044 data files obtained from observations carried out at the Parkes observatory since the year 1991. Around 10^5 files are from surveys of the sky, the remainder are observations of 775 individual pulsars and their corresponding calibration signals. Survey observations are included from the Parkes 70cm and the Swinburne Intermediate Latitude surveys. Individual pulsar observations are included from young pulsar timing projects, the Parkes Pulsar Timing Array and from the PULSE@Parkes outreach program. The data files and access methods are compatible with Virtual Observatory protocols. This paper describes the data currently stored in the archive and presents ways in which these data can be searched and downloaded.

This paper summarizes some of the major calibration and image reconstruction techniques used in radio interferometry and describes them in a common mathematical framework. The use of this framework has a number of benefits, ranging from clarification of the fundamentals, use of standard numerical optimization techniques, and generalization or specialization to new algorithms.

We present the results of an Australia Telescope Compact Array 1.4 GHz spectropolarimetric aperture synthesis survey of 34 square degrees centred on Centaurus A - NGC 5128. A catalogue of 1005 extragalactic compact radio sources in the field to a continuum flux density of 3 mJy/beam is provided along with a table of Faraday rotation measures (RMs) and linear polarised intensities for the 28 percent of sources with high signal-to-noise in linear polarisation. We use the ensemble of 281 background polarised sources as line-of-sight probes of the structure of the giant radio lobes of Centaurus A. This is the first time such a method has been applied to radio galaxy lobes and we explain how it differs from the conventional methods that are often complicated by depth and beam depolarisation effects. We use an RM structure function analysis and report the detection of a turbulent RM signal, with rms of 17 rad/m^2 and scale size 0.3 degrees, associated with the southern giant lobe. We cannot verify whether this signal arises from turbulent structure throughout the lobe or only in a thin skin (or sheath) around the edge, although we favour the latter. The RM signal is modelled as possibly arising from a thin skin with a thermal plasma density equivalent to the Centaurus intragroup medium density and a coherent magnetic field that reverses its sign on a spatial scale of 20 kpc.

A practical evaluation of the Multi-Scale CLEAN algorithm is presented. The data used in the comparisons are taken from The HI Nearby Galaxy Survey (THINGS). The implementation of Multi-Scale CLEAN in the CASA software package is used, although comparisons are made against the very similar Multi-Resolution CLEAN algorithm implemented in AIPS. Both are compared against the classical CLEAN algorithm (as implemented in AIPS). The results of this comparison show that several of the well-known characteristics and issues of using classical CLEAN are significantly lessened (or eliminated completely) when using the Multi-Scale CLEAN algorithm. Importantly, Multi-Scale CLEAN reduces significantly the effects of the clean `bowl' caused by missing short-spacings, and the `pedestal' of low-level un-cleaned flux (which affects flux scales and resolution). Multi-Scale CLEAN can clean down to the noise level without the divergence suffered by classical CLEAN. We discuss practical applications of the added contrast provided by Multi-Scale CLEAN using two selected astronomical examples: HI holes in the interstellar medium and anomalous gas structures outside the main galactic disk.

We consider a troublesome form of non-isoplanatism in synthesis radio telescopes: non-coplanar baselines. We present a novel interpretation of the non-coplanar baselines effect as being due to differential Fresnel diffraction in the neighborhood of the array antennas. We have developed a new algorithm to deal with this effect. Our new algorithm, which we call "W-projection", has markedly superior performance compared to existing algorithms. At roughly equivalent levels of accuracy, W-projection can be up to an order of magnitude faster than the corresponding facet-based algorithms. Furthermore, the precision of result is not tightly coupled to computing time. W-projection has important consequences for the design and operation of the new generation of radio telescopes operating at centimeter and longer wavelengths.

Radio synthesis imaging is dependent upon deconvolution algorithms to counteract the sparse sampling of the Fourier plane. These deconvolution algorithms find an estimate of the true sky brightness from the necessarily incomplete sampled visibility data. The most widely used radio synthesis deconvolution method is the CLEAN algorithm of Hogbom. This algorithm works extremely well for collections of point sources and surprisingly well for extended objects. However, the performance for extended objects can be improved by adopting a multi-scale approach. We describe and demonstrate a conceptually simple and algorithmically straightforward extension to CLEAN that models the sky brightness by the summation of components of emission having different size scales. While previous multiscale algorithms work sequentially on decreasing scale sizes, our algorithm works simultaneously on a range of specified scales. Applications to both real and simulated data sets are given.

Astronomical imaging using aperture synthesis telescopes requires deconvolution of the point spread function as well as calibration of instrumental and atmospheric effects. In general, such effects are time-variable and vary across the field of view as well, resulting in direction-dependent (DD), time-varying gains. Most existing imaging and calibration algorithms assume that the corruptions are direction independent, preventing even moderate dynamic range full-beam, full-Stokes imaging. We present a general framework for imaging algorithms which incorporate DD errors. We describe as well an iterative deconvolution algorithm that corrects known DD errors due to the antenna power patterns and pointing errors for high dynamic range full-beam polarimetric imaging. Using simulations we demonstrate that errors due to realistic primary beams as well as antenna pointing errors will limit the dynamic range of upcoming higher sensitivity instruments and that our new algorithm can be used to correct for such errors. We have applied this algorithm to VLA 1.4 GHz observations of a field that contains two ``4C'' sources and have obtained Stokes-I and -V images with systematic errors that are one order of magnitude lower than those obtained with conventional imaging tools. Our simulations show that on data with no other calibration errors, the algorithm corrects pointing errors as well as errors due to known asymmetries in the antenna pattern.

We have conducted sensitive (1 sigma<30 uJy) 1.4 GHz radio observations with the Australia Telescope Compact Array of a field largely coincident with infrared observations of the Spitzer Wide-Area Extragalactic Survey. The field is centred on the European Large Area ISO Survey S1 region and has a total area of 3.9 deg. We describe the observations and calibration, source extraction, and cross-matching to infrared sources. Two catalogues are presented; one of the radio components found in the image and one of radio sources with counterparts in the infrared and extracted from the literature. 1366 radio components were grouped into 1276 sources, 1183 of which were matched to infrared sources. We discover 31 radio sources with no infrared counterpart at all, adding to the class of Infrared-Faint Radio Sources.

Co-addition of deep (rms about 30 microJy) 20 cm data obtained with the Australia Telescope Compact Array at the location of Spitzer Wide field survey (SWIRE) sources has yielded statistics of radio source counterparts to faint 24 micron sources in stacked images with rms < 1 microJy. We confirm that the infrared-radio correlation extends to f(24micron) = 100 microJy but with a significantly lower coefficient, f (20cm) = 0.039 f(24micron) (or q24 = 1.39) than hitherto reported. We postulate that this may be due to a change in the mean q24 value ratio for objects with f(24micron) < 1mJy.

Deconvolution of the telescope Point Spread Function (PSF) is necessary for even moderate dynamic range imaging with interferometric telescopes. The process of deconvolution can be treated as a search for a model image such that the residual image is consistent with the noise model. For any search algorithm, a parameterized function representing the model such that it fundamentally separates signal from noise will give optimal results. In this paper, the first in a series of forthcoming papers, we argue that in general, spatial correlation length (a measure of the scale of emission) is a stronger separator of the signal from the noise, compared to the strength of the signal alone. Consequently scale sensitive deconvolution algorithms result into more noise-like residuals. We present a scale-sensitive deconvolution algorithm for radio interferometric images, which models the image as a collection of Adaptive Scale Pixels (Asp). Some attempts at optimizing the runtime performance are also presented.

We present deep images of a region around the Crab nebula made with the VLA, utilizing new imaging and deconvolution algorithms in a search for a faint radio shell. The existence of a high-velocity, hydrogen-rich envelope has been predicted to account for the low total mass and kinetic energy of the observed nebula. No radio emission was detected from an extended source outside the Crab nebula. Our limits on surface brightness are sufficiently low to rule out the existence of a shell around the Crab whose brightness is at least two orders of magnitude below SN\thinspace{1006}, the faintest historical shell-type supernova remnant. We consider models for the progenitor star and the pre-supernova environment and conclude that if a fast, outer shock exists then it has a sharply reduced efficiency at accelerating relativistic particles from the kinetic energy of the blast wave. We also looked for a steepening of the spectral index along the boundary of the Crab nebula itself, the signature of an outer shock. However, contrary to previous claims, no such steepening was found. The absence of any evidence at radio wavelengths that either the Crab nebula or a hypothetical shell is interacting with the ambient medium leads to an interpretation that the supernova of 1054 AD was a peculiar low energy event.