We report on the local structure of the Martian subsolar Magnetic Pileup Boundary (MPB) from minimum variance analysis of the magnetic field measured by the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft for six orbits. In particular, we detect a well defined current layer within the MPB's fine structure and provide a local estimate of its current density which results in a sunward Lorentz force. This force accounts for the deflection of the solar wind ions and the acceleration of electrons which carry the interplanetary magnetic field through the MPB into the Magnetic Pileup Region. We find that the thickness of the MPB current layer is of the order of both the upstream (magnetosheath) solar wind proton inertial length and convective gyroradius. This study provides a high resolution view of one of the components of the current system around Mars reported in recent works.

To make a power spectrum (PS) detection of the 21 cm signal from the Epoch of Reionisation (EoR), one must avoid/subtract bright foreground sources. Sources such as Fornax A present a modelling challenge due to spatial structures spanning from arc seconds up to a degree. We compare modelling with multi-scale (MS) CLEAN components to 'shapelets', an alternative set of basis functions. We introduce a new image-based shapelet modelling package, SHAMFI. We also introduce a new CUDA simulation code (WODEN) to generate point source, Gaussian, and shapelet components into visibilities. We test performance by modelling a simulation of Fornax A, peeling the model from simulated visibilities, and producing a residual PS. We find the shapelet method consistently subtracts large-angular-scale emission well, even when the angular-resolution of the data is changed. We find that when increasing the angular-resolution of the data, the MS CLEAN model worsens at large angular-scales. When testing on real MWA data, the expected improvement is not seen in real data because of the other dominating systematics still present. Through further simulation we find the expected differences to be lower than obtainable through current processing pipelines. We conclude shapelets are worthwhile for subtracting extended galaxies, and may prove essential for an EoR detection in the future, once other systematics have been addressed.

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

The Parker Solar Probe was launched on 2018 August 12 and completed its second orbit on 2019 June 19 with perihelion of 35.7 solar radii. During this time, the Energetic particle Instrument-Hi (EPI-Hi, one of the two energetic particle instruments comprising the Integrated Science Investigation of the Sun, ISOIS) measured seven proton intensity increases associated with stream interaction regions (SIRs), two of which appear to be occurring in the same region corotating with the Sun. The events are relatively weak, with observed proton spectra extending to only a few MeV and lasting for a few days. The proton spectra are best characterized by power laws with indices ranging from -4.3 to -6.5, generally softer than events associated with SIRs observed at 1 au and beyond. Helium spectra were also obtained with similar indices, allowing He/H abundance ratios to be calculated for each event. We find values of 0.016-0.031, which are consistent with ratios obtained previously for corotating interaction region events with fast solar wind < 600 km s-1. Using the observed solar wind data combined with solar wind simulations, we study the solar wind structures associated with these events and identify additional spacecraft near 1 au appropriately positioned to observe the same structures after some corotation. Examination of the energetic particle observations from these spacecraft yields two events that may correspond to the energetic particle increases seen by EPI-Hi earlier.

The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $\Delta^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$.

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

We draw a comparison between a solar energetic particle event associated with the release of a slow coronal mass ejection close to the sun, and the energetic particle population produced in high current density field-aligned current structures associated with auroral phenomena in planetary magnetospheres. We suggest that this process is common in CME development and lift-off in the corona, and may account for the electron populations that generate Type III radio bursts, as well as for the prompt energetic ion and electron populations typically observed in interplanetary space.

A solar energetic particle event was detected by the Integrated Science Investigation of the Sun (ISOIS) instrument suite on Parker Solar Probe (PSP) on 2019 April 4 when the spacecraft was inside of 0.17 au and less than 1 day before its second perihelion, providing an opportunity to study solar particle acceleration and transport unprecedentedly close to the source. The event was very small, with peak 1 MeV proton intensities of ~0.3 particles (cm^2 sr s MeV)^-1, and was undetectable above background levels at energies above 10 MeV or in particle detectors at 1 au. It was strongly anisotropic, with intensities flowing outward from the Sun up to 30 times greater than those flowing inward persisting throughout the event. Temporal association between particle increases and small brightness surges in the extreme-ultraviolet observed by the Solar TErrestrial RElations Observatory, which were also accompanied by type III radio emission seen by the Electromagnetic Fields Investigation on PSP, indicates that the source of this event was an active region nearly 80 degrees east of the nominal PSP magnetic footpoint. This suggests that the field lines expanded over a wide longitudinal range between the active region in the photosphere and the corona.

A series of solar energetic particle (SEP) events were observed at Parker Solar Probe (PSP) by the Integrated Science Investigation of the Sun (\ISOIS) during the period from April 18, 2019 through April 24, 2019. The PSP spacecraft was located near 0.48 au from the Sun on Parker spiral field lines that projected out to 1 au within $\sim 25^\circ$ of near Earth spacecraft. These SEP events, though small compared to historically large SEP events, were amongst the largest observed thus far in the PSP mission and provide critical information about the space environment inside 1 au during SEP events. During this period the Sun released multiple coronal mass ejections (CMEs). One of these CMEs observed was initiated on April 20, 2019 at 01:25 UTC, and the interplanetary CME (ICME) propagated out and passed over the PSP spacecraft. Observations by the Electromagnetic Fields Investigation (FIELDS) show that the magnetic field structure was mostly radial throughout the passage of the compression region and the plasma that followed, indicating that PSP did not directly observe a flux rope internal to the ICME, consistent with the location of PSP on the ICME flank. Analysis using relativistic electrons observed near Earth by the Electron, Proton and Alpha Monitor (EPAM) on the Advanced Composition Explorer (ACE) demonstrates the presence of electron seed populations (40--300 keV) during the events observed. The energy spectrum of the \ISOIS~ observed proton seed population below 1 MeV is close to the limit of possible stationary state plasma distributions out of equilibrium. \ISOIS~ observations reveal the \revise{enhancement} of seed populations during the passage of the ICME, which \revise{likely indicates a key part} of the pre-acceleration process that occurs close to the Sun.

The Parker Solar Probe (PSP) spacecraft has flown into the most dense and previously unexplored region of our solar system's zodiacal cloud. While PSP does not have a dedicated dust detector, multiple instruments onboard are sensitive to the effects of meteoroid bombardment. Here, we discuss measurements taken during PSP's first two orbits and compare them to models of the zodiacal cloud's dust distribution. Comparing the radial impact rate trends and the timing and location of a dust impact to an energetic particle detector, we find the impactor population to be consistent with dust grains on hyperbolic orbits escaping the solar system. Assuming PSP's impact environment is dominated by hyperbolic impactors, the total quantity of dust ejected from our solar system is estimated to be 1-14 tons/s. We expect PSP will encounter an increasingly more intense impactor environment as its perihelion distance and semi-major axis are decreased.

The Integrated Science Investigations of the Sun (IS$\odot$IS) instrument suite includes two Energetic Particle instruments: EPI-Hi, designed to measure ions from ~1-200 MeV/nuc, and EPI-Lo, designed to measure ions from ~20 keV/nuc to ~15 MeV/nuc. We present an analysis of eight energetic proton events observed across the energy range of both instruments during PSP's first two orbits in order to examine their combined energy spectra. Background corrections are applied to help resolve spectral breaks between the two instruments and are shown to be effective. In doing so we demonstrate that, even in the early stages of calibration, IS$\odot$IS is capable of producing reliable spectral observations across broad energy ranges. In addition to making groundbreaking measurements very near the Sun, IS$\odot$IS also characterizes energetic particle populations over a range of heliocentric distances inside 1 au. During the first two orbits, IS$\odot$IS observed energetic particle events from a single corotating interaction region (CIR) at three different distances from the Sun. The events are separated by two Carrington rotations and just 0.11 au in distance, however the relationship shown between proton intensities and proximity of the spacecraft to the source region shows evidence of the importance of transport effects on observations of energetic particles from CIRs. Future IS$\odot$IS observations of similar events over larger distances will help disentangle the effects of CIR-related acceleration and transport. We apply similar spectral analyses to the remaining five events, including four that are likely related to stream interaction regions (SIRs) and one solar energetic particle (SEP) event.

Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $\Delta^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude.

We report "ground truth", 28-3500 keV in-situ ion and 5.2-55 keV remotely sensed ENA measurements from Voyager 2/Low Energy Charged Particle (LECP) detector and Cassini/Ion and Neutral Camera (INCA), respectively, that assess the components of the ion pressure in the heliosheath. In this process, we predict an interstellar neutral hydrogen density of ~0.12 cm-3 and an interstellar magnetic field strength of ~0.5 nT upstream of the heliopause in the direction of V2, i.e. consistent with the measured magnetic field and neutral density measurements at Voyager 1 from August 2012, when the spacecraft entered interstellar space, to date. Further, this analysis results in an estimated heliopause crossing by V2 of ~119 AU, as observed, suggesting that the parameters deduced from the pressure analysis are valid. The shape of the >5.2 keV ion energy spectra play a critical role towards determining the pressure balance and acceleration mechanisms inside the heliosheath.

We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.

The Murchison Widefield Array is a low-frequency Square Kilometre Array precursor located at the Murchison Radio-astronomy Observatory in Western Australia. Primarily designed as an imaging telescope, but with a flexible signal path, the capabilities of this telescope have recently been extended to include off-line incoherent and tied-array beam formation using recorded antenna voltages. This has provided the capability for high-time and frequency resolution observations, including a pulsar science program. This paper describes the algorithms and pipeline that we have developed to form the tied array beam products from the summation of calibrated signals of the antenna elements, and presents example polarimetric profiles for PSRs J0437-4715 and J1900-2600 at 185 MHz.

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

We present observations of linear polarisation in the southern radio lobe of Centaurus A, conducted during commissioning of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. We used 16 antennas to observe a 30 square degree region in a single 12 hour pointing over a 240 MHz band centred on 913 MHz. Our observations achieve an angular resolution of $26\times33$ arcseconds (480 parsecs), a maximum recoverable angular scale of 30 arcminutes, and a full-band sensitivity of 85 $\muup$Jy beam$^{-1}$. The resulting maps of polarisation and Faraday rotation are amongst the most detailed ever made for radio lobes, with of order 10$^5$ resolution elements covering the source. We describe several as-yet unreported observational features of the lobe, including its detailed peak Faraday depth structure, and intricate networks of depolarised filaments. These results demonstrate the exciting capabilities of ASKAP for widefield radio polarimetry.

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

Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky signal. We present a general framework for understanding the signatures of ionospheric behaviour in the two-dimensional (2D) neutral hydrogen power spectrum measured by a low-frequency radio interferometer. Two primary classes of ionospheric behaviour are considered, corresponding to dominant modes observed in Murchison Widefield Array (MWA) EoR data; namely, anisotropic structured wave behaviour, and isotropic turbulence. Analytic predictions for power spectrum bias due to this contamination are computed, and compared with simulations. We then apply the ionospheric metric described in Jordan et al. (2017) to study the impact of ionospheric structure on MWA data, by dividing MWA EoR datasets into classes with good and poor ionospheric conditions, using sets of matched 30-minute observations from 2014 September. The results are compared with the analytic and simulated predictions, demonstrating the observed bias in the power spectrum when the ionosphere is active (displays coherent structures or isotropic turbulence). The analysis demonstrates that unless ionospheric activity can be quantified and corrected, active data should not be included in EoR analysis in order to avoid systematic biases in cosmological power spectra. When data are corrected with a model formed from the calibration information, bias reduces below the expected 21~cm signal level. Data are considered `quiet' when the median measured source position offsets are less than 10-15~arcseconds.

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.

In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes from the young Sun. We found that the photochemical loss of atomic hot oxygen dominates over the total ion loss at the current epoch whilst the atmospheric ion loss is likely much more important at ancient times. We briefly discuss the ensuing implications of high atmospheric ion escape rates in the context of ancient Mars, and exoplanets with similar atmospheric compositions around young solar-type stars and M-dwarfs.

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment.

We detail new techniques for analysing ionospheric activity, using Epoch of Reionisation (EoR) datasets obtained with the Murchison Widefield Array (MWA), calibrated by the `Real-Time System' (RTS). Using the high spatial- and temporal-resolution information of the ionosphere provided by the RTS calibration solutions over 19 nights of observing, we find four distinct types of ionospheric activity, and have developed a metric to provide an `at a glance' value for data quality under differing ionospheric conditions. For each ionospheric type, we analyse variations of this metric as we reduce the number of pierce points, revealing that a modest number of pierce points is required to identify the intensity of ionospheric activity; it is possible to calibrate in real-time, providing continuous information of the phase screen. We also analyse temporal correlations, determine diffractive scales, examine the relative fractions of time occupied by various types of ionospheric activity, and detail a method to reconstruct the total electron content responsible for the ionospheric data we observe. These techniques have been developed to be instrument agnostic, useful for application on LOFAR and SKA-Low.

The current generation of experiments aiming to detect the neutral hydrogen signal from the Epoch of Reionisation (EoR) is likely to be limited by systematic effects associated with removing foreground sources from target fields. In this paper we develop a model for the compact foreground sources in one of the target fields of the MWA's EoR key science experiment: the `EoR1' field. The model is based on both the MWA's GLEAM survey and GMRT 150 MHz data from the TGSS survey, the latter providing higher angular resolution and better astrometric accuracy for compact sources than is available from the MWA alone. The model contains 5049 sources, some of which have complicated morphology in MWA data, Fornax A being the most complex. The higher resolution data show that 13% of sources that appear point-like to the MWA have complicated morphology such as double and quad structure, with a typical separation of 33~arcsec. We derive an analytic expression for the error introduced into the EoR two-dimensional power spectrum due to peeling close double sources as single point sources and show that for the measured source properties, the error in the power spectrum is confined to high $k_\bot$ modes that do not affect the overall result for the large-scale cosmological signal of interest. The brightest ten mis-modelled sources in the field contribute 90% of the power bias in the data, suggesting that it is most critical to improve the models of the brightest sources. With this hybrid model we reprocess data from the EoR1 field and show a maximum of 8% improved calibration accuracy and a factor of two reduction in residual power in $k$-space from peeling these sources. Implications for future EoR experiments including the SKA are discussed in relation to the improvements obtained.

Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts.

We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis).

Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between ~20--200 mJy, providing the best limits to date for hour- and month-long transients.

We present new software to cross-match low-frequency radio catalogues: the Positional Update and Matching Algorithm (PUMA). PUMA combines a positional Bayesian probabilistic approach with spectral matching criteria, allowing for confusing sources in the matching process. We go on to create a radio sky model using PUMA based on the Murchison Widefield Array Commissioning Survey, and are able to automatically cross-match 98.5% of sources. Using the characteristics of this sky model, we create simple simulated mock catalogues on which to test PUMA, and find that PUMA can reliably find the correct spectral indices of sources, along with being able to recover ionospheric offsets. Finally, we use this sky model to calibrate and remove foreground sources from simulated interferometric data, generated using OSKAR (the Oxford University visibility generator). We demonstrate that there is a substantial improvement in foreground source removal when using higher frequency and higher resolution source positions, even when correcting positions by an average of 0.3 given a synthesized beam-width of 2.3.

Context: For over 12 yr, we have carried out a precise radial velocity survey of a sample of 373 G and K giant stars using the Hamilton \'Echelle Spectrograph at Lick Observatory. There are, among others, a number of multiple planetary systems in our sample as well as several planetary candidates in stellar binaries. Aims: We aim at detecting and characterizing substellar+stellar companions to the giant star HD 59686 A (HR 2877, HIP 36616). Methods: We obtained high precision radial velocity (RV) measurements of the star HD 59686 A. By fitting a Keplerian model to the periodic changes in the RVs, we can assess the nature of companions in the system. In order to discriminate between RV variations due to non-radial pulsation or stellar spots we used infrared RVs taken with the CRIRES spectrograph at the Very Large Telescope. Additionally, to further characterize the system, we obtain high-resolution images with LMIRCam at the Large Binocular Telescope. Results: We report the likely discovery of a giant planet with a mass of $m_{p}~\sin i=6.92_{-0.24}^{+0.18}~M_{Jup}$ orbiting at $a_{p}=1.0860_{-0.0007}^{+0.0006}$ au from the giant star HD 59686 A. Besides the planetary signal, we discover an eccentric ($e_{B}=0.729_{-0.003}^{+0.004}$) binary companion with a mass of $m_{B}~\sin i=0.5296_{-0.0008}^{+0.0011}~M_{Sun}$ orbiting at a semi-major axis of just $a_{B}=13.56_{-0.14}^{+0.18}$ au. Conclusions: The existence of the planet HD 59686 Ab in a tight eccentric binary system severely challenges standard giant planet formation theories and requires substantial improvements to such theories in tight binaries. Otherwise, alternative planet formation scenarios such as second generation planets or dynamical interactions in an early phase of the system's lifetime should be seriously considered in order to better understand the origin of this enigmatic planet.

Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307,455 radio sources with 20 separate flux density measurements across 72--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5\sigma$, which itself is typically 50mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low frequency flux density scale of the southern sky to better than 10%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results.

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.

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 $\Delta^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.

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.

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.

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.

Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year.

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.

We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < {\delta} < 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.

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.

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.

The New Horizons spacecraft carried three instruments that measured the space environment near Pluto as it flew by on 14 July 2015. The Solar Wind Around Pluto instrument revealed an interaction region confined sunward of Pluto to within about 6 Pluto radii. The surprisingly small size is consistent with a reduced atmospheric escape rate as well as a particularly high solar wind flux. The Pluto Energetic Particle Spectrometer Science Investigation (PEPSSI) observations suggested ions are accelerated and-or deflected around Pluto. In the wake of the interaction region PEPSSI observed suprathermal particle fluxes about one tenth the flux in the interplanetary medium, increasing with distance downstream. The Student Dust Counter, which measures radius greater than 1.4 um grains, detected 1 candidate impact from 5days before to 5 days after closest approach, indicating an upper limit for the dust density in the Pluto system of 4.6 per cubic km.

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.

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

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.

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.

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.

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.

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$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

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$\sigma$ 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 $\alpha<-1.2$ ($S\propto\nu^{\alpha}$). 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.

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.

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.

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$\sigma$ and 6.5$\sigma$ respectively. We detected 51$\%$ of the MWA giant pulses at the Parkes radio telescope, with spectral indices in the range of $-3.6>\alpha> -4.9$ ($S_{\rm \nu} \propto \nu^\alpha$). 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.

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.

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.

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.

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.

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.

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 $\alpha$ = $-1.17\pm0.06$ and $-1.37\pm0.08$ for the west and east relics, respectively ($S \propto \nu^{\alpha}$). 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.

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.

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40 kHz frequency resolution and 0.5 s time resolution; the second year using 10 kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing and image weighting scheme. At 154 MHz the image resolution is approximately 2.5 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA's primary beam.

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.

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.

The solar system contains solids of all sizes, ranging from km-size bodies to nano-sized particles. Nanograins have been detected in situ in the Earth's atmosphere, near cometary and giant planet environments, and more recently in the solar wind at 1 AU. These latter nano grains are thought to be formed in the inner solar system dust cloud, mainly through collisional break-up of larger grains and are then picked-up and accelerated by the magnetized solar wind because of their large charge-to-mass ratio. In the present paper, we analyze the low frequency bursty noise identified in the Cassini radio and plasma wave data during the spacecraft cruise phase inside Jupiter's orbit. The magnitude, spectral shape and waveform of this broadband noise is consistent with the signature of nano particles impinging at nearby the solar wind speed on the spacecraft surface. Nanoparticles were observed whenever the radio instrument was turned on and able to detect them, at different heliocentric distances between Earth and Jupiter, suggesting their ubiquitous presence in the heliosphere. We analyzed the radial dependence of the nano dust flux with heliospheric distance and found that it is consistent with the dynamics of nano dust originating from the inner heliosphere and picked-up by the solar wind. The contribution of the nano dust produced in asteroid belt appears to be negligible compared to the trapping region in the inner heliosphere. In contrast, further out, nano dust are mainly produced by the volcanism of active moons such as Io and Enceladus.

An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements.

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.

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.

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

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.

Recent observations with the Murchison Widefield Array at 185~MHz have serendipitously unveiled a heretofore unknown giant and relatively nearby ($z = 0.0178$) radio galaxy associated with NGC\,1534. The diffuse emission presented here is the first indication that NGC\,1534 is one of a rare class of objects (along with NGC\,5128 and NGC\,612) in which a galaxy with a prominent dust lane hosts radio emission on scales of $\sim$700\,kpc. We present details of the radio emission along with a detailed comparison with other radio galaxies with disks. NGC1534 is the lowest surface brightness radio galaxy known with an estimated scaled 1.4-GHz surface brightness of just 0.2\,mJy\,arcmin$^{-2}$. The radio lobes have one of the steepest spectral indices yet observed: $\alpha=-2.1\pm0.1$, and the core to lobe luminosity ratio is $<0.1$\%. We estimate the space density of this low brightness (dying) phase of radio galaxy evolution as $7\times10^{-7}$\,Mpc$^{-3}$ and argue that normal AGN cannot spend more than 6\% of their lifetime in this phase if they all go through the same cycle.

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.

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.

After 9 years in the Saturn system, the Cassini spacecraft finally observed Titan in the supersonic solar wind. These unique observations reveal that Titan interaction with the solar wind is in many ways similar to un-magnetized planets Mars and Venus in spite of the differences in the properties of the solar plasma in the outer solar system. In particular, Cassini detected a collisionless, supercritical bow shock and a well-defined induced magnetosphere filled with mass-loaded interplanetary magnetic field lines, which drape around Titan ionosphere. Although the flyby altitude may not allow the detection of an ionopause, Cassini reports enhancements of plasma density compatible with plasma clouds or streamers in the flanks of its induced magnetosphere or due to an expansion of the induced magnetosphere. Because of the upstream conditions, these observations are also relevant for unmagnetized bodies in the outer solar system such as Pluto, where kinetic processes are expected to dominate.

We present the results of an approximately 6,100 square degree 104--196MHz radio sky survey performed with the Murchison Widefield Array during instrument commissioning between 2012 September and 2012 December: the Murchison Widefield Array Commissioning Survey (MWACS). The data were taken as meridian drift scans with two different 32-antenna sub-arrays that were available during the commissioning period. The survey covers approximately 20.5 h < Right Ascension (RA) < 8.5 h, -58 deg < Declination (Dec) < -14 deg over three frequency bands centred on 119, 150 and 180 MHz, with image resolutions of 6--3 arcmin. The catalogue has 3-arcmin angular resolution and a typical noise level of 40 mJy/beam, with reduced sensitivity near the field boundaries and bright sources. We describe the data reduction strategy, based upon mosaiced snapshots, flux density calibration and source-finding method. We present a catalogue of flux density and spectral index measurements for 14,110 sources, extracted from the mosaic, 1,247 of which are sub-components of complexes of sources.

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.

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.

The detection of the Epoch of Reionization (EoR) in the redshifted 21-cm line is a challenging task. Here we formulate the detection of the EoR signal using the drift scan strategy. This method potentially has better instrumental stability as compared to the case where a single patch of sky is tracked. We demonstrate that the correlation time between measured visibilities could extend up to 1-2 hr for an interferometer array such as the Murchison Widefield Array (MWA), which has a wide primary beam. We estimate the EoR power based on cross-correlation of visibilities across time and show that the drift scan strategy is capable of the detection of the EoR signal with comparable/better signal-to-noise as compared to the tracking case. We also estimate the visibility correlation for a set of bright point sources and argue that the statistical inhomogeneity of bright point sources might allow their separation from the EoR signal.

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

Astronomical widefield imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new widefield interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependencies of CASA's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarisation correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2-12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square-Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.

Nanodust grains of a few nanometer in size are produced near the Sun by collisional break-up of larger grains and picked-up by the magnetized solar wind. They have so far been detected at 1 AU by only the two STEREO spacecraft. Here we analyze the spectra measured by the radio and plasma wave instrument onboard Cassini during the cruise phase close to Earth orbit; they exhibit bursty signatures similar to those observed by the same instrument in association to nanodust stream impacts on Cassini near Jupiter. The observed wave level and spectral shape reveal impacts of nanoparticles at about 300 km/s, with an average flux compatible with that observed by the radio and plasma wave instrument onboard STEREO and with the interplanetary flux models.

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.

We present techniques for bridging the gap between idealized inverse covariance weighted quadratic estimation of 21 cm power spectra and the real-world challenges presented universally by interferometric observation. By carefully evaluating various estimators and adapting our techniques for large but incomplete data sets, we develop a robust power spectrum estimation framework that preserves the so-called "EoR window" and keeps track of estimator errors and covariances. We apply our method to observations from the 32-tile prototype of the Murchinson Widefield Array to demonstrate the importance of a judicious analysis technique. Lastly, we apply our method to investigate the dependence of the clean EoR window on frequency--especially the frequency dependence of the so-called "wedge" feature--and establish upper limits on the power spectrum from z = 6.2 to z = 11.7. Our lowest limit is Delta(k) < 0.3 Kelvin at 95% confidence at a comoving scale k = 0.046 Mpc^-1 and z = 9.5.

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.

In this paper we explore for the first time the relative magnitudes of three fundamental sources of uncertainty, namely, foreground contamination, thermal noise and sample variance in detecting the HI power spectrum from the Epoch of Reionization (EoR). We derive limits on the sensitivity of a Fourier synthesis telescope to detect EoR based on its array configuration and a statistical representation of images made by the instrument. We use the Murchison Widefield Array (MWA) configuration for our studies. Using a unified framework for estimating signal and noise components in the HI power spectrum, we derive an expression for and estimate the contamination from extragalactic point-like sources in three-dimensional k-space. Sensitivity for EoR HI power spectrum detection is estimated for different observing modes with MWA. With 1000 hours of observing on a single field using the 128-tile MWA, EoR detection is feasible (S/N > 1 for $k\lesssim 0.8$ Mpc$^{-1}$). Bandpass shaping and refinements to the EoR window are found to be effective in containing foreground contamination, which makes the instrument tolerant to imaging errors. We find that for a given observing time, observing many independent fields of view does not offer an advantage over a single field observation when thermal noise dominates over other uncertainties in the derived power spectrum.

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.

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.

From the 27th to the 28th January 2009, the Cassini spacecraft remotely acquired combined observations of Saturns southern aurorae at radio, ultraviolet and infrared wavelengths, while monitoring ion injections in the middle magnetosphere from energetic neutral atoms. Simultaneous measurements included the sampling of a full planetary rotation, a relevant timescale to investigate auroral emissions driven by processes internal to the magnetosphere. In addition, this interval coincidently matched a powerful substorm-like event in the magnetotail, which induced an overall dawnside intensification of the magnetospheric and auroral activity. We comparatively analyze this unique set of measurements to reach a comprehensive view of kronian auroral processes over the investigated timescale. We identify three source regions in atmospheric aurorae, including a main oval associated with the bulk of Saturn Kilometric Radiation (SKR), together with polar and equatorward emissions. These observations reveal the co-existence of corotational and sub-corototational dynamics of emissions associated with the main auroral oval. Precisely, we show that the atmospheric main oval hosts short-lived sub-corotating isolated features together with a bright, longitudinally extended, corotating region locked at the southern SKR phase. We assign the susbtorm-like event to a regular, internally driven, nightside ion injection possibly triggerred by a plasmoid ejection. We also investigate the total auroral energy budget, from the power input to the atmosphere, characterized by precipitating electrons up to 20 keV, to its dissipation through the various radiating processes. Finally, through simulations, we confirm the search-light nature of the SKR rotational modulation and we show that SKR arcs relate to isolated auroral spots. The resulting findings are discussed in the frame of pending questions.

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 $(\alpha,\delta) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

We aim to detect and characterize substellar companions to K giant stars, to further our knowledge of planet formation and stellar evolution of intermediate-mass stars. For more than a decade we have used Doppler spectroscopy to acquire high precision radial velocity measurements of K giant stars. All data for this survey have been taken at Lick Observatory. Our survey includes 373 G and K giants. Radial velocity data showing periodic variations are fitted with Keplerian orbits using a chi^2 minimization technique. We report the presence of two substellar companions to the K giant stars tau Gem and 91 Aqr. The brown dwarf orbiting tau Gem has an orbital period of 305.5 +/- 0.1 days, a minimum mass of 20.6 M_Jup, and an eccentricity of 0.031 +/- 0.009. The planet orbiting 91 Aqr has an orbital period of 181.4 +/- 0.1 days, a minimum mass of 3.2 M_Jup, and an eccentricity of 0.027 +/- 0.026. Both companions have exceptionally circular orbits for their orbital distance, as compared to all previously discovered planetary companions.

Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives.

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.

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.

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.

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{\sigma} detection of the EoR signal with one full season of observation on two fields (900 and 700 hours).

We introduce the Fast Holographic Deconvolution method for analyzing interferometric radio data. Our new method is an extension of A-projection/software-holography/forward modeling analysis techniques and shares their precision deconvolution and widefield polarimetry, while being significantly faster than current implementations that use full direction-dependent antenna gains. Using data from the MWA 32 antenna prototype, we demonstrate the effectiveness and precision of our new algorithm. Fast Holographic Deconvolution may be particularly important for upcoming 21 cm cosmology observations of the Epoch of Reionization and Dark Energy where foreground subtraction is intimately related to the precision of the data reduction.

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.

The first station of the Long Wavelength Array (LWA1) was completed in April 2011 and is currently performing observations resulting from its first call for proposals in addition to a continuing program of commissioning and characterization observations. The instrument consists of 258 dual-polarization dipoles, which are digitized and combined into beams. Four independently-steerable dual-polarization beams are available, each with two "tunings" of 16 MHz bandwidth that can be independently tuned to any frequency between 10 MHz and 88 MHz. The system equivalent flux density for zenith pointing is ~3 kJy and is approximately independent of frequency; this corresponds to a sensitivity of ~5 Jy/beam (5sigma, 1 s); making it one of the most sensitive meter-wavelength radio telescopes. LWA1 also has two "transient buffer" modes which allow coherent recording from all dipoles simultaneously, providing instantaneous all-sky field of view. LWA1 provides versatile and unique new capabilities for Galactic science, pulsar science, solar and planetary science, space weather, cosmology, and searches for astrophysical transients. Results from LWA1 will detect or tightly constrain the presence of hot Jupiters within 50 parsecs of Earth. LWA1 will provide excellent resolution in frequency and in time to examine phenomena such as solar bursts, and pulsars over a 4:1 frequency range that includes the poorly understood turnover and steep-spectrum regimes. Observations to date have proven LWA1's potential for pulsar observing, and just a few seconds with the completed 256-dipole LWA1 provide the most sensitive images of the sky at 23 MHz obtained yet. We are operating LWA1 as an open skies radio observatory, offering ~2000 beam-hours per year to the general community.

The process of wide-field synthesis imaging is explored, with the aim of understanding the implications of variable, polarised primary beams for forthcoming Epoch of Reionisation experiments. These experiments seek to detect weak signatures from redshifted 21cm emission in deep residual datasets, after suppression and subtraction of foreground emission. Many subtraction algorithms benefit from low side-lobes and polarisation leakage at the outset, and both of these are intimately linked to how the polarised primary beams are handled. Building on previous contributions from a number of authors, in which direction-dependent corrections are incorporated into visibility gridding kernels, we consider the special characteristics of arrays of fixed dipole antennas operating around 100-200 MHz, looking towards instruments such as the Square Kilometre Array (SKA) and the Hydrogen Epoch of Reionization Arrays (HERA). We show that integrating snapshots in the image domain can help to produce compact gridding kernels, and also reduce the need to make complicated polarised leakage corrections during gridding. We also investigate an alternative form for the gridding kernel that can suppress variations in the direction-dependent weighting of gridded visibilities by 10s of dB, while maintaining compact support.

The Murchison Widefield Array (MWA) is a new low-frequency, wide field-of-view radio interferometer under development at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. We have used a 32-element MWA prototype interferometer (MWA-32T) to observe two 50-degree diameter fields in the southern sky in the 110 MHz to 200 MHz band in order to evaluate the performance of the MWA-32T, to develop techniques for epoch of reionization experiments, and to make measurements of astronomical foregrounds. We developed a calibration and imaging pipeline for the MWA-32T, and used it to produce ~15' angular resolution maps of the two fields. We perform a blind source extraction using these confusion-limited images, and detect 655 sources at high significance with an additional 871 lower significance source candidates. We compare these sources with existing low-frequency radio surveys in order to assess the MWA-32T system performance, wide field analysis algorithms, and catalog quality. Our source catalog is found to agree well with existing low-frequency surveys in these regions of the sky and with statistical distributions of point sources derived from Northern Hemisphere surveys; it represents one of the deepest surveys to date of this sky field in the 110 MHz to 200 MHz band.

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

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Shapley 1 (Sp 1), which is known to contain a close-binary central star system. Close-binary central stars have been identified as a likely source of shaping in planetary nebulae, but with little observational support to date. Deep narrowband imaging in the light of [O III] {\lambda}5007A suggests the presence of a large bow-shock to the west of the nebula, indicating that it is undergoing the first stages of an interaction with the interstellar medium. Further narrowband imaging in the light of H{\alpha}+[NII] {\lambda}6584A combined with longslit observations of the H{\alpha} emission have been used to develop a spatio-kinematical model of Sp 1. The model clearly reveals Sp 1 to be a bipolar, axisymmetric structure viewed almost pole-on. The symmetry axis of the model nebula is within a few degrees of perpendicular to the orbital plane of the central binary system - strong evidence that the central close-binary system has played an important role in shaping the nebula. Sp 1 is one of very few nebulae to have this link, between nebular symmetry axis and binary plane, shown observationally.

We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emission features are evident. Our data represent a significant advance in low radio frequency solar imaging, enabling us to follow the spatial, spectral, and temporal evolution of events simultaneously and in unprecedented detail. The rich variety of features seen here reaffirms the coronal diagnostic capability of low radio frequency emission and provides an early glimpse of the nature of radio observations that will become available as the next generation of low frequency radio interferometers come on-line over the next few years.

We present a method for subtracting point sources from interferometric radio images via forward modeling of the instrument response and involving an algebraic nonlinear minimization. The method is applied to simulated maps of the Murchison Wide-field Array but is generally useful in cases where only image data are available. After source subtraction, the residual maps have no statistical difference to the expected thermal noise distribution at all angular scales, indicating high effectiveness in the subtraction. Simulations indicate that the errors in recovering the source parameters decrease with increasing signal-to-noise ratio, which is consistent with the theoretical measurement errors. In applying the technique to simulated snapshot observations with the Murchison Wide-field Array, we found that all 101 sources present in the simulation were recovered with an average position error of 10 arcsec and an average flux density error of 0.15%. This led to a dynamic range increase of approximately 3 orders of magnitude. Since all the sources were deconvolved jointly, the subtraction was not limited by source sidelobes but by thermal noise. This technique is a promising deconvolution method for upcoming radio arrays with a huge number of elements, and a candidate for the difficult task of subtracting foreground sources from observations of the 21 cm neutral Hydrogen signal from the epoch of reionization.

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

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell~41, which is known to contain the well-studied close-binary system MT Ser. This object represents an important test case in the study of the evolution of planetary nebulae with binary central stars as current evolutionary theories predict that the binary plane should be aligned perpendicular to the symmetry axis of the nebula. Longslit observations of the \NII\ emission from Abell~41 were obtained using the Manchester Echelle Spectrometer on the 2.1-m San Pedro M\'artir Telescope. These spectra, combined with deep, narrowband imagery acquired using ACAM on the William Herschel Telescope, were used to develop a spatio-kinematical model of \NII\ emission from Abell~41. The best fitting model reveals Abell~41 to have a waisted, bipolar structure with an expansion velocity of $\sim$40\kms{} at the waist. The symmetry axis of the model nebula is within 5$^\circ$ of perpendicular to the orbital plane of the central binary system. This provides strong evidence that the close-binary system, MT Ser, has directly affected the shaping of its host nebula, Abell~41.

It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.

Bright point sources associated with extragalactic AGN and radio galaxies are an important foreground for low frequency radio experiments aimed at detecting the redshifted 21cm emission from neutral hydrogen during the epoch of reionization. The frequency dependence of the synthesized beam implies that the sidelobes of these sources will move across the field of view as a function of observing frequency, hence frustrating line-of-sight foreground subtraction techniques. We describe a method for subtracting these point sources from dirty maps produced by an instrument such as the MWA. This technique combines matched filters with an iterative centroiding scheme to locate and characterize point sources in the presence of a diffuse background. Simulations show that this technique can improve the dynamic range of EOR maps by 2-3 orders of magnitude.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell 41, which is known to contain the well-studied close-binary system MT Ser. This object represents an important test case in the study of the evolution of planetary nebulae with binary central stars as current evolutionary theories predict that the binary plane should be aligned perpendicular to the symmetry axis of the nebula. Deep narrowband imaging in the light of [NII], [OIII] and [SII], obtained using ACAM on the William Herschel Telescope, has been used to investigate the ionisation structure of Abell 41. Longslit observations of the H-alpha and [NII] emission were obtained using the Manchester Echelle Spectrometer on the 2.1-m San Pedro M\'artir Telescope. These spectra, combined with the narrowband imagery, were used to develop a spatio-kinematical model of [NII] emission from Abell 41. The best fitting model reveals Abell 41 to have a waisted, bipolar structure with an expansion velocity of ~40km\s at the waist. The symmetry axis of the model nebula is within 5$\degr$ of perpendicular to the orbital plane of the central binary system. This provides strong evidence that the close-binary system, MT Ser, has directly affected the shaping of its nebula, Abell 41. Although the theoretical link between bipolar planetary nebulae and binary central stars is long established, this nebula is only the second to have this link, between nebular symmetry axis and binary plane, proved observationally.

The Murchison Widefield Array (MWA) is a next-generation radio telescope currently under construction in the remote Western Australia Outback. Raw data will be generated continuously at 5GiB/s, grouped into 8s cadences. This high throughput motivates the development of on-site, real time processing and reduction in preference to archiving, transport and off-line processing. Each batch of 8s data must be completely reduced before the next batch arrives. Maintaining real time operation will require a sustained performance of around 2.5TFLOP/s (including convolutions, FFTs, interpolations and matrix multiplications). We describe a scalable heterogeneous computing pipeline implementation, exploiting both the high computing density and FLOP-per-Watt ratio of modern GPUs. The architecture is highly parallel within and across nodes, with all major processing elements performed by GPUs. Necessary scatter-gather operations along the pipeline are loosely synchronized between the nodes hosting the GPUs. The MWA will be a frontier scientific instrument and a pathfinder for planned peta- and exascale facilities.

In radio astronomy, reference signals from auxiliary antennas that receive only the radio frequency interference (RFI) can be modified to model the RFI environment at the astronomy receivers. The RFI can then be canceled from the astronomy signal paths. However, astronomers typically only require signal statistics. If the RFI statistics are changing slowly, the cancellation can be applied to the signal correlations at a much lower rate than is required for standard adaptive filters. In this paper we describe five canceler setups; precorrelation and postcorrelation cancelers that use one or two reference signals in different ways. The theoretical residual RFI and added noise levels are examined and are demonstrated using microwave television RFI at the Australia Telescope Compact Array. The RFI is attenuated to below the system noise, a reduction of at least 20 dB. While dual-reference cancelers add more reference noise than single-reference cancelers, this noise is zero-mean and only adds to the system noise, decreasing the sensitivity. The residual RFI that remains in the output of single-reference cancelers (but not dual-reference cancelers) sets a nonzero noise floor that does not act like random system noise and may limit the achievable sensitivity. Thus, dual-reference cancelers often result in superior cancellation. Dual-reference precorrelation cancelers require a double-canceler setup to be useful and to give equivalent results to dual-reference postcorrelation cancelers.

We investigate characteristics of radio frequency interference (RFI) signals that can affect the excision potential of some interference mitigation algorithms. The techniques considered are those that modify signals from auxiliary reference antennas to model and cancel interference from an astronomical observation. These techniques can be applied in the time domain, where the RFI voltage is modeled and subtracted from the astronomy signal path (adaptive noise canceling), or they can be applied to the autocorrelated and cross-correlated voltage spectra in the frequency domain (postcorrelation canceling). For ideal receivers and a single, statistically stationary interfering signal, both precorrelation and postcorrelation filters can result in complete cancellation of the interference from the observation. The postcorrelation method has the advantage of being applied on tens or hundreds of millisecond timescales rather than tens or hundreds of nanosecond timescales. However, this can be a disadvantage if the RFI transmitter location is changing, since the cross-correlated power measurements which link the interference power in the astronomy and reference signal paths can decorrelate. If the decorrelation is not too severe, it can be allowed for, at the expense of a noise increase. The time domain adaptive cancelers are allowed to slightly vary their internal coefficients and adapt to changing phases during the integrations, which means that they avoid the decorrelation problem. However, the freedom to adapt also results in a noise increase. In this paper the ability of both types of cancelers to excise interference originating from a moving source is compared. The cancelers perform well on both observed and simulated data, giving complete cancellation.

Because energetic particles populate both planetary magnetospheres and interplanetary space in significant quantities, energetic-ion sensors have been flown since the beginning of the space age. Early sensors were solid-state detector (SSD) telescopes, with conical fields of view, often swept through a circle by virtue of the spin motion of the spacecraft (e.g., IMP 7 and 8, ISEE 1 and 2). In the 1980s and 1990s, foil/microchannel plate (MCP) time-of-flight (TOF) measurements were added to the energy measurement provided by the SSD (eg, AMPTE/CCE MEPA, Geotail EPIC/ICS, Galileo EPD). The resulting energy and velocity uniquely identified ion mass. More recently, we have developed a 2-D fan acceptance angle sensor that includes both energy and TOF. When mounted on a spinning spacecraft, this 160^\circ x 12^\circ FOV sweeps out nearly 4\pi steradians in one spin. This sensor, dubbed the "hockey puck" for its shape, is currently in flight on MESSENGER (EPS) and New Horizons Pluto (PEPPSI).Increasingly, energetic-ion sensors fly on 3-axis stabilized spacecraft (e.g., MESSENGER EPS, New Horizons (Pluto) PEPPSI, Cassini MIMI. While 3-axis stabilization serves imaging science well, it hampers the goal of obtaining 4\pi-steradian ion measurements. We are developing an energetic-ion sensor that measures ion energy and composition, and covers 2\pi steradians on a 3-axis-stabilized spacecraft without an articulation mechanism. Based on its shape, we refer to this design as the "mushroom". We describe the internally funded development of the concept and its status at the start of development funding by NASA under the Planetary Instrument Definition and Development Program (PIDDP).

We present the first detailed spatio-kinematical analysis and modelling of the Southern planetary nebula SuWt 2. This object presents a problem for current theories of planetary nebula formation and evolution, as it is not known to contain a central post-main sequence star. Deep narrowband [NII]6584 images reveal the presence of faint bipolar lobes emanating fromthe edges of the nebular ring. Longslit observations of the H-alpha and [NII]6584 emission lines were obtained using EMMI on the 3.6-m ESO-NTT. The spectra reveal the nebular morphology as a bright torus encircling the waist of an extended bipolar structure. By deprojection, the inclination of the ring is found to be 68$\degr$ $\pm$ 2$\degr$ (c.f. ~90$\degr$ for the double A-type binary believed to lie at the centre of the nebula), and the ring expansion velocity is found to be 28 km/s. Our findings are discussed with relation to possible formation scenarios for SuWt 2. Through comparison of the nebular heliocentric systemic velocity, found here to be -25 $\pm$ 5 km/s, and the heliocentric systemic velocity of the double A-type binary, we conclude that neither component of the binary could have been the nebular progenitor. However, we are unable to rule out the presence of a third component to the system, which would have been the nebula progenitor.

We present a detailed kinematical analysis of the young compact hourglass-shaped planetary nebula Hb 12. We performed optical imaging and longslit spectroscopy of Hb 12 using the Manchester echelle spectrometer with the 2.1m San Pedro Martir telescope. We reveal, for the first time, the presence of end caps (or knots) aligned with the bipolar lobes of the planetary nebula shell in a deep [NII]6584 image of Hb 12. We measured from our spectroscopy radial velocities of 120 km/s for these knots. We have derived the inclination angle of the hourglass shaped nebular shell to be 65 degrees to the line of sight. It has been suggested that Hb 12's central star system is an eclipsing binary (Hsia et al. 2006) which would imply a binary inclination of at least 80 degrees. However, if the central binary has been the major shaping influence on the nebula then both nebula and binary would be expected to share a common inclination angle. Finally, we report the discovery of high-velocity knots with Hubble-type velocities, close to the core of Hb 12, observed in Halpha and oriented in the same direction as the end caps. Very different velocities and kinematical ages were calculated for the outer and inner knots showing that they may originate from different outburst events.

Time-to-science is an important figure of merit for digital instrumentation serving the astronomical community. A digital signal processing (DSP) community is forming that uses shared hardware development, signal processing libraries, and instrument architectures to reduce development time of digital instrumentation and to improve time-to-science for a wide variety of projects. We suggest prioritizing technological development supporting the needs of this nascent DSP community. After outlining several instrument classes that are relying on digital instrumentation development to achieve new science objectives, we identify key areas where technologies pertaining to interoperability and processing flexibility will reduce the time, risk, and cost of developing the digital instrumentation for radio astronomy. These areas represent focus points where support of general-purpose, open-source development for a DSP community should be prioritized in the next decade. Contributors to such technological development may be centers of support for this DSP community, science groups that contribute general-purpose DSP solutions as part of their own instrumentation needs, or engineering groups engaging in research that may be applied to next-generation DSP instrumentation.

The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities.

A soft ultraviolet radiation field, 10.2 eV < E <13.6 eV, that permeates neutral intergalactic gas during the Epoch of Reionization (EoR) excites the 2p (directly) and 2s (indirectly) states of atomic hydrogen. Because the 2s state is metastable, the lifetime of atoms in this level is relatively long, which may cause the 2s state to be overpopulated relative to the 2p state. It has recently been proposed that for this reason, neutral intergalactic atomic hydrogen gas may be detected in absorption in its 3-cm fine-structure line (2s_1/2 -> 2p_3/2) against the Cosmic Microwave Background out to very high redshifts. In particular, the optical depth in the fine-structure line through neutral intergalactic gas surrounding bright quasars during the EoR may reach tau~1e-5. The resulting surface brightness temperature of tens of micro K (in absorption) may be detectable with existing radio telescopes. Motivated by this exciting proposal, we perform a detailed analysis of the transfer of Lyman beta,gamma,delta,... radiation, and re-analyze the detectability of the fine-structure line in neutral intergalactic gas surrounding high-redshift quasars. We find that proper radiative transfer modeling causes the fine-structure absorption signature to be reduced tremendously to tau< 1e-10. We therefore conclude that neutral intergalactic gas during the EoR cannot reveal its presence in the 3-cm fine-structure line to existing radio telescopes.

The MWA is a next-generation radio interferometer under construction in remote Western Australia. The data rate from the correlator makes storing the raw data infeasible, so the data must be processed in real-time. The processing task is of order ~10 TFLOPS. The remote location of the MWA limits the power that can be allocated to computing. We describe the design and implementation of elements of the MWA real-time data processing system which leverage the computing abilities of modern graphics processing units (GPUs). The matrix algebra and texture mapping capabilities of GPUs are well suited to the majority of tasks involved in real-time calibration and imaging. Considerable performance advantages over a conventional CPU-based reference implementation are obtained.

The interferometric technique known as peeling addresses many of the challenges faced when observing with low-frequency radio arrays, and is a promising tool for the associated calibration systems. We investigate a real-time peeling implementation for next-generation radio interferometers such as the Murchison Widefield Array (MWA). The MWA is being built in Australia and will observe the radio sky between 80 and 300 MHz. The data rate produced by the correlator is just over 19 GB/s (a few Peta-Bytes/day). It is impractical to store data generated at this rate, and software is currently being developed to calibrate and form images in real time. The software will run on-site on a high-throughput real-time computing cluster at several tera-flops, and a complete cycle of calibration and imaging will be completed every 8 seconds. Various properties of the implementation are investigated using simulated data. The algorithm is seen to work in the presence of strong galactic emission and with various ionospheric conditions. It is also shown to scale well as the number of antennas increases, which is essential for many upcoming instuments. Lessons from MWA pipeline development and processing of simulated data may be applied to future low-frequency fixed dipole arrays.

Since 1999, we have been conducting a radial velocity survey of 179 K giants using the CAT at UCO/Lick observatory. At present ~20-100 measurements have been collected per star with a precision of 5 to 8 m/s. Of the stars monitored, 145 (80%) show radial velocity (RV) variations at a level >20 m/s, of which 43 exhibit significant periodicities. Our aim is to investigate possible mechanism(s) that cause these observed RV variations. We intend to test whether these variations are intrinsic in nature, or possibly induced by companions, or both. In addition, we aim to characterise the parameters of these companions. A relation between log g and the amplitude of the RV variations is investigated for all stars in the sample. Furthermore, the hypothesis that all periodic RV variations are caused by companions is investigated by comparing their inferred orbital statistics with the statistics of companions around main sequence stars. A strong relation is found between the amplitude of the RV variations and log g in K giant stars, as suggested earlier by Hatzes & Cochran (1998). However, most of the stars exhibiting periodic variations are located above this relation. These RV variations can be split in a periodic component which is not correlated with log g and a random residual part which does correlate with log g. Compared to main-sequence stars, K giants frequently exhibit periodic RV variations. Interpreting these RV variations as being caused by companions, the orbital param eters are different from the companions orbiting dwarfs. Intrinsic mechanisms play an important role in producing RV variations in K giants stars, as suggested by their dependence on log g. However, it appears that periodic RV variations are additional to these intrinsic variations, consistent with them being caused by companions.

Since 1999, a radial velocity survey of 179 red giant stars is ongoing at Lick Observatory with a one month cadence. At present ~20-100 measurements have been collected per star with an accuracy of 5 to 8 m/s. Of the stars monitored, 145 (80%) show radial velocity (RV) variations at a level >20 m/s, of which 43 exhibit significant periodicities. Here, we investigate the mechanism causing the observed radial velocity variations. Firstly, we search for a correlation between the radial velocity amplitude and an intrinsic parameter of the star, in this case surface gravity (log g). Secondly, we investigate line profile variations and compare these with theoretical predictions.

The Pluto Energetic Particle Spectrometer Science Investigation (PEPSSI) comprises the hardware and accompanying science investigation on the New Horizons spacecraft to measure pick-up ions from Pluto's outgassing atmosphere. To the extent that Pluto retains its characteristics similar to those of a "heavy comet" as detected in stellar occultations since the early 1980s, these measurements will characterize the neutral atmosphere of Pluto while providing a consistency check on the atmospheric escape rate at the encounter epoch with that deduced from the atmospheric structure at lower altitudes by the ALICE, REX, and SWAP experiments on New Horizons. In addition, PEPSSI will characterize any extended ionosphere and solar wind interaction while also characterizing the energetic particle environment of Pluto, Charon, and their associated system. First proposed for development for the Pluto Express mission in September 1993, what became the PEPSSI instrument went through a number of development stages to meet the requirements of such an instrument for a mission to Pluto while minimizing the required spacecraft resources. The PEPSSI instrument provides for measurements of ions (with compositional information) and electrons from 10s of keV to ~1 MeV in a 120 deg x 12 deg fan-shaped beam in six sectors for 1.5 kg and ~2.5 W.

Context: The Rossiter-McLaughlin (RM) effect, a rotational effect in eclipsing systems, provides unique insight into the relative orientation of stellar spin axes and orbital axes of eclipsing binary systems. Aims: Our aim is to develop a robust method to analyze the RM effect in an eclipsing system with two nearly equally bright components. This gives access to the orientation of the stellar rotation axes and may shed light on questions of binary formation and evolution. Methods: High-resolution spectra have been obtained both out of eclipse and during the primary and secondary eclipses in the V1143Cyg system, using the high-resolution Hamilton Echelle Spectrograph at the Lick Observatory. The Rossiter-McLaughlin effect is analyzed in two ways: (1) by measuring the shift of the line center of gravity during different phases of the eclipses and (2) by analysis of the line shape change of the rotational broadening function during eclipses. Results: The projected axes of both stars are aligned with the orbital spin within the observational uncertainties, with the angle of the primary rotation axis beta_p=0.3+-1.5 deg, and the angle of the secondary rotation axis beta_s=-1.2+-1.6 deg, thereby showing that the remaining difference between the theoretical and observed apsidal motion for this system is not due to a misalignment of the stellar rotation axes. Both methods utilized in this paper work very well, even at times when the broadening profiles of the two stars overlap.[abridged]

We study the motion of stars in a star cluster which revolves in a circular orbit about its parent galaxy. The star cluster is modelled as an ellipsoid of uniform spatial density. We exhibit two 2-parameter families of self-consistent equilibrium models in which the velocity at each point is confined to a line in velocity space. We exhibit the link between this problem and that of a uniform rotating ellipsoidal galaxy. With minimal adaptation, Freeman's bar models yield a third family.

We present the first detailed kinematical analysis of the planetary nebula Abell 63, which is known to contain the eclipsing close-binary nucleus UU Sge. Abell 63 provides an important test case in investigating the role of close-binary central stars on the evolution of planetary nebulae. Longslit observations were obtained using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. The spectra reveal that the central bright rim of Abell 63 has a tube-like structure. A deep image shows collimated lobes extending from the nebula, which are shown to be high-velocity outflows. The kinematic ages of the nebular rim and the extended lobes are calculated to be 8400+/-500 years and 12900+/-2800 years, respectively, which suggests that the lobes were formed at an earlier stage than the nebular rim. This is consistent with expectations that disk-generated jets form immediately after the common envelope phase. A morphological-kinematical model of the central nebula is presented and the best-fit model is found to have the same inclination as the orbital plane of the central binary system; this is the first proof that a close-binary system directly affects the shaping of its nebula. A Hubble-type flow is well-established in the morphological-kinematical modelling of the observed line profiles and imagery. Two possible formation models for the elongated lobes of Abell 63 are considered (1) a low-density, pressure-driven jet excavates a cavity in the remnant AGB envelope; (2) high-density bullets form the lobes in a single ballistic ejection event.

It has long been speculated that the observed periodic radial velocity pattern for the K giant Pollux might be explained in terms of an orbiting planetary companion. We have collected 80 high-resolution spectra for Pollux at Lick Observatory yielding precise radial velocities with a mean error of 3.8 m/s, providing the most comprehensive and precise data set available for this star. Our data confirm the periodicity previously seen in the radial velocities. We derive a period of 589.7+/-3.5 days and, assuming a primary mass of 1.86 M_Sun, a minimum companion mass of 2.9+/-0.3 M_Jup, consistent with earlier determinations. No evidence for any periodicities is visible in our analysis of the shapes of the spectral lines via the bisector method, so that we conclude that evidence is accumulating and compelling for a planet around Pollux. However, some last doubt remains about this interpretation, because non-radial pulsations which might be present in giant stars could in principle also explain the observed radial velocities, while the accompanying bisector variations might be too small to be detectable with current data.

Context: Future astrometric missions such as SIM PlanetQuest need very stable reference stars. K giants have large luminosities, which place them at large distances and thus the jitter of their photocenters by companions is relatively small. Therefore K giants would be best suited as references. To confirm this observationally a radial velocity survey is performed to quantify the level of intrinsic variability in K giants. Aims: From this radial velocity survey we present 34 K giants with an observed standard deviation of the radial velocity of less than 20 m/s. These stars are considered ``stable'' and can be used as radial velocity standards. Methods: The radial velocity survey contains 179 K giants. All K giants have a declination between -30 and +65 degrees and visual magnitude of 3-6 mag. The Coude Auxiliary Telescope (CAT) at UCO/Lick Observatory is used to obtain radial velocities with an accuracy of 5-8 m/s. The number of epochs for the 34 stable stars ranges from 11 to 28 with a total timespan of the observations between 1800 and a little over 2200 days. Results: The observational results of the 34 ``stable'' stars are shown together with a discussion about their position in the M_{V} vs. B-V diagram and some conclusions concerning the radial velocity variability of K giants. These results are in agreement with the theoretical predictions. K giants in a certain range of the M_{V} vs. B-V diagram are suitable reference stars.

A detailed analysis is presented of the dominant ionised knot in the halo of the planetary nebula NGC 6543. Observations were made at high spectral and spatial resolution of the [OIII] 5007 line using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. A 20-element multislit was stepped across the field to give almost complete spatial coverage of the knot and surrounding halo. The spectra reveal, for the first time, gas flows around the kinematically inert knot. The gas flows are found to have velocities comparable to the sound speed as gas is photo-evaporated off an ionised surface. No evidence is found of fast wind interaction with the knot and we find it likely that the fast wind is still contained in a pressure-driven bubble in the core of the nebula. This rules out the possibility of the knot having its origin in instabilities at the interface of the fast and AGB winds. We suggest that the knot is embedded in the slowly expanding Red Giant wind and that its surfaces are being continually photoionised by the central star.

We present results of our survey for planetary transits in the field of NGC 6940. We think nearly all of our observed stars are field stars. We have obtained high precision (3-10 millimags at the bright end) photometric observations of 50,000 stars spanning 18 nights in an attempt to identify low amplitude and short period transit events. We have used a matched filter analysis to identify 14 stars that show multiple events, and four stars that show single transits. Of these 18 candidates, we have identified two that should be further researched. However, none of the candidates are convincing hot Jupiters.

This paper is the first of a series describing the results of the Australia Telescope Hubble Deep Field South (ATHDFS) radio survey. The survey was conducted at four wavelengths - 20, 11, 6, and 3 cm, over a 4-year period, and achieves an rms sensitivity of about 10 microJy at each wavelength. We describe the observations and data reduction processes, and present data on radio sources close to the centre of the HDF-S. We discuss in detail the properties of a subset of these sources. The sources include both starburst galaxies and galaxies powered by an active galactic nucleus, and range in redshift from 0.1 to 2.2. Some of them are characterised by unusually high radio-to-optical luminosities, presumably caused by dust extinction.

We present results from 30 nights of observations of the open cluster NGC 7789 with the WFC camera on the INT telescope in La Palma. From ~900 epochs, we obtained lightcurves and Sloan r-i colours for ~33000 stars, with ~2400 stars with better than 1% precision. We expected to detect ~2 transiting hot Jupiter planets if 1% of stars host such a companion and that a typical hot Jupiter radius is ~1.2RJ. We find 24 transit candidates, 14 of which we can assign a period. We rule out the transiting planet model for 21 of these candidates using various robust arguments. For 2 candidates we are unable to decide on their nature, although it seems most likely that they are eclipsing binaries as well. We have one candidate exhibiting a single eclipse for which we derive a radius of 1.81+0.09-0.00RJ. Three candidates remain that require follow-up observations in order to determine their nature.

We present results from our survey for planetary transits in the field of the intermediate age (~2.5 Gyr), metal-rich ([Fe/H]~+0.07) open cluster NGC 6819. We have obtained high-precision time-series photometry for over 38,000 stars in this field and have developed an effective matched-filter algorithm to search for photometric transits. This algorithm identified 8 candidate stars showing multiple transit-like events, plus 3 stars with single eclipses. On closer inspection, while most are shown to be low mass stellar binaries, some of these events could be due to brown dwarf companions. The data for one of the single-transit candidates indicates a minimum radius for the companion similar to that of HD 209458b.

Archetypal analysis represents each individual member of a set of data vectors as a mixture (a constrained linear combination) of the pure types or archetypes of the data set. The archetypes are themselves required to be mixtures of the data vectors. Archetypal analysis may be particularly useful in analysing data sets comprising galaxy spectra, since each spectrum is, presumably, a superposition of the emission from the various stellar populations, nebular emissions and nuclear activity making up that galaxy, and each of these emission sources corresponds to a potential archetype of the entire data set. We demonstrate archetypal analysis using sets of composite synthetic galaxy spectra, showing that the method promises to be an effective and efficient way to classify spectra. We show that archetypal analysis is robust in the presence of various types of noise.

We are using the Isaac Newton Telescope Wide Field Camera to survey open cluster fields for transiting hot Jupiter planets. Clusters were selected on the basis of visibility, richness of stars, age and metallicity. Observations of NGC 6819, 6940 and 7789 began in 1999 and continued in 2000. We have developed an effective matched-filter transit-detection algorithm which has proved its ability to identify very low amplitude eclipse events in real data. Here we present our results for NGC 6819. We have identified 7 candidates showing transit-like events. Colour information suggests that most of the companion bodies are likely to be very-low-mass stars or brown dwarfs, intrinsically interesting objects in their own right.