Breakthrough Listen is a ten-year initiative to search for signatures of technologies created by extraterrestrial civilizations at radio and optical wavelengths. Here, we detail the digital data recording system deployed for Breakthrough Listen observations at the 64-m aperture CSIRO Parkes Telescope in New South Wales, Australia. The recording system currently implements two recording modes: a dual-polarization, 1.125 GHz bandwidth mode for single beam observations, and a 26-input, 308-MHz bandwidth mode for the 21-cm multibeam receiver. The system is also designed to support a 3 GHz single-beam mode for the forthcoming Parkes ultra-wideband feed. In this paper, we present details of the system architecture, provide an overview of hardware and software, and present initial performance results.

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 use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.

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 present extensive observations of the radio emission from the remnant of SN 1987A made with the Australia Telescope Compact Array (ATCA), since the first detection of the remnant in 1990. The radio emission has evolved in time providing unique information on the interaction of the supernova shock with the circumstellar medium. We particularly focus on the monitoring observations at 1.4, 2.4, 4.8 and 8.6 GHz, which have been made at intervals of 4-6 weeks. The flux density data show that the remnant brightness is now increasing exponentially, while the radio spectrum is flattening. The current spectral index value of -0.68 represents an 18+/-3% increase over the last 8 years. The exponential trend in the flux is also found in the ATCA imaging observations at 9 GHz, which have been made since 1992, approximately twice a year, as well as in the 843 MHz data set from the Molonglo Observatory Synthesis Telescope from 1987 to March 2007. Comparisons with data at different wavelengths (X-ray, H\alpha) are made. The rich data set that has been assembled in the last 22 years forms a basis for a better understanding of the evolution of the supernova remnant.

The aftermath of supernova (SN) 1987A continues to provide spectacular insights into the interaction between a SN blastwave and its circumstellar en- vironment. We here present 36 GHz observations from the Australia Telescope Compact Array of the radio remnant of SN 1987A. These new images, taken in 2008 Apr and 2008 Oct, substantially extend the frequency range of an ongo- ing monitoring and imaging program conducted between 1.4 and 20 GHz. Our 36.2 GHz images have a diffraction-limited angular resolution of 0.3-0.4 arcseconds, which covers the gap between high resolution, low dynamic range VLBI images of the remnant and low resolution, high dynamic range images at frequencies between 1 and 20 GHz. The radio morphology of the remnant at 36 GHz is an elliptical ring with enhanced emission on the eastern and western sides, similar to that seen previously at lower frequencies. Model fits to the data in the Fourier domain show that the emitting region is consistent with a thick inclined torus of mean radius 0.85 arcsec, and a 2008 Oct flux density of 27 +/- 6 mJy at 36.2 GHz. The spectral index for the remnant at this epoch, determined between 1.4 GHz and 36.2 GHz, is -0.83. There is tentative evidence for an unresolved central source with flatter spectral index.

We present detailed Fourier modeling of the radio remnant of Supernova 1987A, using high-resolution 9 GHz and 18 GHz data taken with the Australia Telescope Compact Array over the period 1992 to 2008. We develop a parameterized three-dimensional torus model for the expanding radio shell, in which the emission is confined to an inclined equatorial belt; our model also incorporates both a correction for light travel-time effects and an overall east-west gradient in the radio emissivity. By deriving an analytic expression for the two-dimensional Fourier transform of the projected three-dimensional brightness distribution, we can fit our spatial model directly to the interferometric visibility data. This provides robust estimates to the radio morphology at each epoch. The best-fit results suggest a constant remnant expansion at 4000 +/- 400 km/s over the 16-year period covered by the observations. The model fits also indicate substantial mid-latitude emission, extending to 40 degree on either side of the equatorial plane. This likely corresponds to the extra-planar structure seen in H$\alpha$ and Ly$\alpha$ emission from the supernova reverse shock, and broadly supports hydrodynamic models in which the complex circumstellar environment was produced by a progression of interacting winds from the progenitor. Our model quantifies the clear asymmetry seen in the radio images: we find that the eastern half of the radio remnant is consistently ~40 brighter than the western half at all epochs, which may result from an asymmetry in the ejecta distribution between these two hemispheres.

Supernova 1987A in the Large Magellanic Cloud provides a spectacularly detailed view of the aftermath of a core-collapse explosion. The supernova ejecta initially coasted outward at more than 10% of the speed of light, but in 1990 were observed to decelerate rapidly as they began to encounter dense circumstellar material expelled by the progenitor star. The resulting shock has subsequently produced steadily brightening radio synchrotron emission, which is resolved by the Australia Telescope Compact Array (ATCA) into an expanding limb-brightened shell. Here we present 15 years of ATCA imaging of Supernova 1987A, at an effective angular resolution of 0.4 arcsec. We find that the radio remnant has accelerated in its expansion over this period, from approx 3600 km/s in 1992 to approx 5200 km/s at the end of 2006. The published diameters of the evolving X-ray shell have been ~15% smaller than the corresponding radio values, but a simultaneous Fourier analysis of both radio and X-ray data eliminates this discrepancy, and yields a current diameter for the shell in both wave-bands of approx 1.7 arcsec. An asymmetric brightness distribution is seen in radio images at all ATCA epochs: the eastern and western rims have higher fluxes than the northern and southern regions, indicating that most of the radio emission comes from the equatorial plane of the system, where the progenitor star's circumstellar wind is thought to be densest. The eastern lobe is brighter than and further from the supernova site than the western lobe, suggesting an additional asymmetry in the initial distribution of supernova ejecta.

The high-energy gamma-ray emission discovered using the H.E.S.S. telescopes from the binary system PSR B1259 -63, is modelled using an extension of the approach that successfully predicted it. We find that the simultaneous INTEGRAL and H.E.S.S. data permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. Additional, simultaneous, X-ray and TeV data sets are needed to lift this degeneracy.

The high-energy gamma-ray light curve of the binary system PSR B1259 -63, is computed using the approach that successfully predicted the spectrum at periastron. The simultaneous INTEGRAL and H.E.S.S. spectra taken 16 days after periastron currently permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. In this paper we shown how the long-term light curve may help to lift this degeneracy.

We report here on extensive radio observations of the binary system containing PSR B1259-63 and the Be star SS 2883, made around the time of the 2004 periastron. This is the fourth periastron to have been observed in detail. As in previous observations, changes in the pulsar's dispersion and rotation measures are detected over a period spanning 200 days. We show that the eclipse of the pulsed emission lasts from 16 days prior to periastron to 15 days after periastron and is consistent from one periastron to the next. We demonstrate that the timing solution proposed by Wang et al. (2004) provides a good fit through the 2004 periastron. The light curve of the transient unpulsed radio emission is broadly similar from one periastron to the next. For this periastron, however, the light curve is strongly peaked post-periastron with rather little enhancement prior to periastron, in contrast to the 2000 periastron where the peak flux densities were more equal. These observations remain consistent with the interpretation that the pulsar passes through the dense circumstellar disk of the Be star just before and just after periastron. The observed differences from one periastron to the next can be ascribed to variations in the local disk density and magnetic field structure at the time the pulsar enters the disk.

We present calculations of the radio images and light curves from supernovae, based on high-resolution numerical simulations of the hydrodynamics and radiation transfer in a spherically symmetric medium. As a specific example we model the emission from SN1993J. This supernova does not appear to be expanding in a self-similar fashion, and cannot be adequately fitted with the often-used analytic mini-shell model. We present a good fit to the radio evolution at a single frequency. Both free-free absorption and synchrotron self-absorption are needed to fit the light curve at early times, and a circumstellar density profile of $\rho \sim r ^{-1.7}$ provides the best fit to the later data. Comparisons of VLBI images of SN1993J with synthetic model images suggest that internal free-free absorption completely obscures emission at 8.4 GHz passing through the center of the supernova for the first few tens of years after explosion.

We present calculations of the radio emission from supernovae based on high-resolution simulations of the hydrodynamics and radiation transfer, using simple energy density relations which link the properties of the radiating electrons and the magnetic field to the hydrodynamics. As a specific example we model the emission from SN1993J, which cannot be adequately fitted with the often-used analytic mini-shell model, and present a good fit to the radio evolution at a single frequency. Both free-free absorption and syncrotron self- absorption are needed to fit the light curve at early times and a circumstellar density profile of rho~r^-1.7 provides the best fit to the later data. We show that the interaction of density structures in the ejecta with the reverse supernova shock may produce features in the radio light curves such as have been observed. We discuss the use of high-resolution radio images of supernovae to distinguish between different absorption mechanisms and determine the origin of specific light curve features. Comparisons of VLBI images of SN1993J with synthetic model images suggest that internal free-free absorption completely obscures emission at 8.4 GHz passing through the center of the supernova for the first few tens of years after explosion. We predict that at 8.4 GHz the internal free-free absorption is currently declining, and that over the next 40 years the surface brightness of the center of the source should increase relative to the bright ring of emission seen in VLBI images. Similar absorption in a nearby supernova would make the detection of a radio pulsar at 1 GHz impossible for ~150 years after explosion.

PSR B1259-63 is a 47ms radio pulsar with a high spin-down luminosity which is in a close, highly eccentric 3.5-year orbit about a bright stellar companion. The binary system may be a detectable source of hard gamma-rays produced by inverse Compton scattering of photons from the B2e star SS2883 by electrons and positrons in the pulsar wind. The star provides an enormous density of optical photons in the vicinity of the pulsar, particularly at epochs near periastron. We calculate the emission from the unshocked region of the pulsar wind, assuming that it terminates at a shock where it attains pressure balance with the companion's wind. The spectra and light curves for the inverse Compton emission from the shock-terminated wind are compared with those for an unterminated wind. If the pulsar's wind is weaker than that from the companion star, the termination of the wind decreases the inverse Compton flux, particularly near periastron. The termination shock geometry has the effect of decreasing the asymmetry of the gamma-ray light curve around periastron, which arises because of the asymmetrical variation of the scattering angle.

We report the flux densities of the evolving radio source SN1987A at 843 MHz measured from observations made with the Molonglo Observatory Synthesis Telescope between 1994 September and 2000 May. The radio light curve shows that the rate of increase of the flux density jumped markedly around days 2800-3000 (i.e. in the first half of 1995), and that since then the radio evolution has been remarkably well fitted by a simple linear increase of (62.7 +- 0.5) microJansky/day. We discuss in detail the relationship between the radio light curve and the recent brightening of the system at optical wavelengths.

We summarise the proceedings of the SRCfTA workshop on ``Supernova Remnants, Pulsars and the Interstellar Medium'' that was held at the University of Sydney on Mar 18 and 19, 1999.

We investigate the effects of inverse Compton scattering by electrons and positrons in the unshocked winds of rotationally-powered binary pulsars. This process can scatter low energy target photons to produce gamma rays with energies from MeV to TeV. The binary radio pulsars PSR B1259-63 and PSR J0045-73 are both in close eccentric orbits around bright main sequence stars which provide a huge density of low energy target photons. The inverse Compton scattering process transfers momentum from the pulsar wind to the scattered photons, and therefore provides a drag which tends to decelerate the pulsar wind. We present detailed calculations of the dynamics of a pulsar wind which is undergoing inverse Compton scattering, showing that the deceleration of the wind of PSR B1259-63 due to `inverse Compton drag' is small, but that this process may confine the wind of PSR J0045-73 before it attains pressure balance with the outflow of its companion star. We calculate the spectra and light curves of the resulting inverse Compton emission from PSR B1259-63 and show that if the size of the pulsar wind nebula is comparable to the binary separation, then the gamma-ray emission from the unshocked wind may be detectable by atmospheric Cerenkov detectors or by the new generation of satellite-borne gamma-ray detectors such as INTEGRAL and GLAST. This mechanism may therefore provide a direct probe of the freely-expanding regions of pulsar winds, previously thought to be invisible.

We discuss the interpretation of transient, unpulsed radio emission detected from the unique pulsar/Be-star binary system PSR B1259-63. Extensive monitoring of the 1994 and 1997 periastron passages has shown that the source flares over a 100-day interval around periastron, varying on time-scales as short as a day and peaking at 60 mJy (~100 times the apastron flux density) at 1.4 GHz. Interpreting the emission as synchrotron radiation, we show that (i) the observed variations in flux density are too large to be caused by the shock interaction between the pulsar wind and an isotropic, radiatively driven, Be-star wind, and (ii) the radio emitting electrons do not originate from the pulsar wind. We argue instead that the radio electrons originate from the circumstellar disk of the Be star and are accelerated at two epochs, one before and one after periastron, when the pulsar passes through the disk. A simple model incorporating two epochs of impulsive acceleration followed by synchrotron cooling reproduces the essential features of the radio light curve and spectrum and is consistent with the system geometry inferred from pulsed radio data.

We derive light curves for the hard gamma-ray emission, at energies up to several TeV, expected from the unique pulsar/Be-star binary system PSR B1259-63. This is the only known system in our galaxy in which a radio pulsar is orbiting a main sequence star. We show that inverse Compton emission from the electrons and positrons in the shocked pulsar wind, scattering target photons from the Be star, produces a flux of hard gamma-rays that should be above the sensitivity threshold of present day atmospheric Cerenkov detectors. Furthermore, we predict that the flux of hard gamma-rays produced via this mechanism has a characteristic variation with orbital phase that should be observable, and which is not expected from any other mechanism.

We propose a model for radio supernovae (RSN) based on the synchrotron emission from relativistic electrons which are diffusively accelerated at the expanding supernova shock. This model was originally developed for application to the optically thin emission observed from SN1987A. Here we generalise it by including the effects of free-free absorption from both an external screen and from material internal to the source, and by relaxing the restriction to an azimuthal B-field. We find a good fit to the entire set of radio data for the best observed highly-luminous RSN -- SN1986J -- with a reduced \chi^2 of 3.85. Applying the new model to SN1988Z, another intrinsically bright RSN, also yields a good fit (\chi^2_{\rm red}\approx 2) but this is less significant, because of the limited data on this distant (z=0.02) source. These fits suggest that the shock expands at constant speed, that the magnetic field within the source decreases with time according to t^{-2}, and that the compression ratio of the shock front is close to the value expected of a strong shock in an ideal gas of adiabatic index 5/3 -- indicating a relatively low value of the cosmic ray pressure compared with SN1987A. In the case of SN1986J we derive an explosion date in August/September 1982, a magnetic field at the position of the shock 1000 days after explosion of 4nT and a spatial diffusion coefficient of the electrons of 4 10^{19} cm^2/s, four orders of magnitude greater than the Bohm value. In addition, we obtain the optical depths to external and internal absorption, and derive an estimate of the mass-loss rate.

The observations of radio emission from SNR~1987A can be accounted for on the basis of diffusive shock acceleration of electrons by the supernova blast wave. However, with this interpretation the observed spectral index implies that the compression ratio of the gas subshock is roughly $2.7$ rather than the value of $4$ expected of a strong shock front. We propose that in SNR~1987A, ions also undergo diffusive acceleration at the shock, a process that is likely to be rapid. Unlike the electron population, the accelerated ions can have an important effect on the gas dynamics. We calculate this coupled gas and energetic particle dynamics on the basis of the two-fluid model, in which the accelerated ions provide an additional component to the total pressure acting on the fluid. By accelerating and possibly heating the upstream plasma, the initially strong shock is modified and a weaker subshock with an upstream precursor results. The electrons behave as test particles. They are accelerated at the evolving subshock, escape downstream, and emit synchrotron radiation in the swept up magnetic field. Two models are considered for the surroundings of the progenitor: that of a freely expanding wind of number density $n\propto r^{-2}$, and that of a wind confined by a shell of denser material, creating a stagnation zone of roughly constant density beyond the standing shock which terminates the free wind. We model the observed radio light curves and the relatively steep spectrum of SNR~1987A using similar values for the ion acceleration parameters to those used in models of cosmic ray acceleration in older SNRs which can also contain high Mach number shocks, and find a good fit for the case in which the termination shock is located at about $2\times 10^{15}\,$m from the progenitor.

When applied to the blast wave formed by the explosion of a massive star as a supernova (SN), the theory of diffusive particle acceleration at shock fronts predicts a very high energy density in cosmic rays. Almost immediately after particles begin to be injected into the process, the cosmic ray pressure rises until comparable to the ram-pressure encountered by the shock front. Those supernovae which are observed in the radio band i.e., radio supernovae (RSNe), provide direct evidence of particle acceleration in the form of synchrotron emitting electrons. Furthermore, these objects are particularly interesting, since they are usually surrounded by a relatively dense confining medium. The acceleration of cosmic rays can then lead to the production of very high energy (VHE) gamma-rays which arise from collisions between energetic particles and target nuclei. We estimate the cosmic ray energy density assuming a fraction $\phi\lesim1$ of the energy available at the shock front is converted into cosmic rays. Combining this with the parameters describing the environment of the SN progenitor, as deduced from observations, and from more detailed modelling, we compute the flux at Earth $\tevflux$ of photons of energy above $1\,$TeV. For the relatively weak but nearby supernova SN1987A we predict $\tevflux=2\times10^{-13}\,\fluxunits$ before the shock front encounters the ring of dense matter seen by the Hubble Space Telescope. Subsequently, the flux is expected to rise further. The medium around SN1993J in M81 is thought to have a density