We present a multi-wavelength study of the active galactic nucleus in the nearby ($D=14.1$ Mpc) low mass galaxy IC 750, which has circumnuclear 22 GHz water maser emission. The masers trace a nearly edge-on, warped disk $\sim$0.2 pc in diameter, coincident with the compact nuclear X-ray source which lies at the base of the $\sim$kpc-scale extended X-ray emission. The position-velocity structure of the maser emission indicates the central black hole (BH) has a mass less than $1.4 \times 10^5~M_\odot$. Keplerian rotation curves fitted to these data yield enclosed masses between $4.1 \times 10^4~M_\odot$ and $1.4 \times 10^5~M_\odot$, with a mode of $7.2 \times 10^4~M_\odot$. Fitting the optical spectrum, we measure a nuclear stellar velocity dispersion $\sigma_* = 110.7^{+12.1}_{-13.4}$~{\rm km~s}$^{-1}.$ From near-infrared photometry, we fit a bulge mass of $(7.3 \pm 2.7) \times 10^8~M_\odot$ and a stellar mass of $1.4 \times 10^{10}~M_\odot$. The mass upper limit of the intermediate mass black hole in IC 750 falls roughly two orders of magnitude below the $M_{\rm BH}-\sigma_*$ relation and roughly one order of magnitude below the $M_{\rm BH}-M_{\rm Bulge}$ and $M_{\rm BH}-M_*$ relations -- larger than the relations' intrinsic scatters of (0.58 $\pm$ 0.09) dex, 0.69 dex, and (0.65 $\pm$ 0.09) dex, respectively. These offsets could be due to larger scatter at the low mass end of these relations. Alternatively, black hole growth is intrinsically inefficient in galaxies with low bulge and/or stellar masses, which causes the black holes to be under-massive relative to their hosts, as predicted by some galaxy evolution simulations.

A successful ground array Radio Frequency (RF)-only self-trigger on 10 high-energy cosmic ray events is demonstrated with 256 dual-polarization antennas of the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). This RF-only capability is predicated on novel techniques for Radio Frequency Interference (RFI) identification and mitigation with an analysis efficiency of 45\% for shower-driven events with a Signal-to-noise ratio $\gtrsim$ 5 against the galactic background noise power of individual antennas. This technique enables more efficient detection of cosmic rays over a wider range of zenith angles than possible via triggers from in-situ particle detectors and can be easily adapted by neutrino experiments relying on RF-only detection. This paper discusses the system design, RFI characterization and mitigation techniques, and initial results from 10 cosmic ray events identified within a 40-hour observing window. A design for a future optimized commensal cosmic-ray detector for the OVRO-LWA is presented, as well as recommendations for developing a similar capability for other experiments -- these designs either reduce data-rate or increase sensitivity by an order of magnitude for many configurations of radio instruments.

The 21\,cm transition of neutral hydrogen is opening an observational window into the cosmic dawn of the universe---the epoch of first star formation. We use 28\,hr of data from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) to place upper limits on the spatial power spectrum of 21\,cm emission at $z \approx 18.4$ ($\Delta_{21} \lesssim 10^4\,\text{mK}$), and within the absorption feature reported by the EDGES experiment (Bowman et al. 2018). In the process we demonstrate the first application of the double Karhunen-Lo\`{e}ve transform for foreground filtering, and diagnose the systematic errors that are currently limiting the measurement. We also provide an updated model for the angular power spectrum of low-frequency foreground emission measured from the northern hemisphere, which can be used to refine sensitivity forecasts for next-generation experiments.

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.

The Large-Aperture Experiment to Detect the Dark Age (LEDA) was designed to detect the predicted O(100)mK sky-averaged absorption of the Cosmic Microwave Background by Hydrogen in the neutral pre- and intergalactic medium just after the cosmological Dark Age. The spectral signature would be associated with emergence of a diffuse Ly$\alpha$ background from starlight during 'Cosmic Dawn'. Recently, Bowman et al. (2018) have reported detection of this predicted absorption feature, with an unexpectedly large amplitude of 530 mK, centered at 78 MHz. Verification of this result by an independent experiment, such as LEDA, is pressing. In this paper, we detail design and characterization of the LEDA radiometer systems, and a first-generation pipeline that instantiates a signal path model. Sited at the Owens Valley Radio Observatory Long Wavelength Array, LEDA systems include the station correlator, five well-separated redundant dual polarization radiometers and backend electronics. The radiometers deliver a 30-85MHz band (16<z<34) and operate as part of the larger interferometric array, for purposes ultimately of in situ calibration. Here, we report on the LEDA system design, calibration approach, and progress in characterization as of January 2016. The LEDA systems are currently being modified to improve performance near 78 MHz in order to verify the purported absorption feature.

We have conducted the most sensitive low frequency (below 100 MHz) search to date for prompt, low-frequency radio emission associated with short-duration gamma-ray bursts (GRBs), using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). The OVRO-LWA's nearly full-hemisphere field-of-view ($\sim20$,$000$ square degrees) allows us to search for low-frequency (sub-$100$ MHz) counterparts for a large sample of the subset of GRB events for which prompt radio emission has been predicted. Following the detection of short GRB 170112A by Swift, we used all-sky OVRO-LWA images spanning one hour prior to and two hours following the GRB event to search for a transient source coincident with the position of GRB 170112A. We detect no transient source, with our most constraining $1\sigma$ flux density limit of $650~\text{mJy}$ for frequencies spanning $27~\text{MHz}-84~\text{MHz}$. We place constraints on a number of models predicting prompt, low-frequency radio emission accompanying short GRBs and their potential binary neutron star merger progenitors, and place an upper limit of $L_\text{radio}/L_\gamma \lesssim 7\times10^{-16}$ on the fraction of energy released in the prompt radio emission. These observations serve as a pilot effort for a program targeting a wider sample of both short and long GRBs with the OVRO-LWA, including bursts with confirmed redshift measurements which are critical to placing the most constraining limits on prompt radio emission models, as well as a program for the follow-up of gravitational wave compact binary coalescence events detected by advanced LIGO and Virgo.

A host of new low-frequency radio telescopes seek to measure the 21-cm transition of neutral hydrogen from the early universe. These telescopes have the potential to directly probe star and galaxy formation at redshifts $20 \gtrsim z \gtrsim 7$, but are limited by the dynamic range they can achieve against foreground sources of low-frequency radio emission. Consequently, there is a growing demand for modern, high-fidelity maps of the sky at frequencies below 200 MHz for use in foreground modeling and removal. We describe a new widefield imaging technique for drift-scanning interferometers, Tikhonov-regularized $m$-mode analysis imaging. This technique constructs images of the entire sky in a single synthesis imaging step with exact treatment of widefield effects. We describe how the CLEAN algorithm can be adapted to deconvolve maps generated by $m$-mode analysis imaging. We demonstrate Tikhonov-regularized $m$-mode analysis imaging using the Owens Valley Long Wavelength Array (OVRO-LWA) by generating 8 new maps of the sky north of $\delta=-30^\circ$ with 15 arcmin angular resolution, at frequencies evenly spaced between 36.528 MHz and 73.152 MHz, and $\sim$800 mJy/beam thermal noise. These maps are a 10-fold improvement in angular resolution over existing full-sky maps at comparable frequencies, which have angular resolutions $\ge 2^\circ$. Each map is constructed exclusively from interferometric observations and does not represent the globally averaged sky brightness. Future improvements will incorporate total power radiometry, improved thermal noise, and improved angular resolution -- due to the planned expansion of the OVRO-LWA to 2.6 km baselines. These maps serve as a first step on the path to the use of more sophisticated foreground filters in 21-cm cosmology incorporating the measured angular and frequency structure of all foreground contaminants.

Radio astronomy observatories with high throughput back end instruments require real-time data processing. While computing hardware continues to advance rapidly, development of real-time processing pipelines remains difficult and time-consuming, which can limit scientific productivity. Motivated by this, we have developed Bifrost: an open-source software framework for rapid pipeline development. Bifrost combines a high-level Python interface with highly efficient reconfigurable data transport and a library of computing blocks for CPU and GPU processing. The framework is generalizable, but initially it emphasizes the needs of high-throughput radio astronomy pipelines, such as the ability to process data buffers as if they were continuous streams, the capacity to partition processing into distinct data sequences (e.g., separate observations), and the ability to extract specific intervals from buffered data. Computing blocks in the library are designed for applications such as interferometry, pulsar dedispersion and timing, and transient search pipelines. We describe the design and implementation of the Bifrost framework and demonstrate its use as the backbone in the correlation and beamforming back end of the Long Wavelength Array station in the Sevilleta National Wildlife Refuge, NM.

We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-shifted emission in the north and west arms. Comparisons with previous images of the 28SiO v=1,2, J=1-0 transitions at 7mm (observed in 2001-2002) show that the bulk of the J=2-1 transition emission follows the streamlines of the J=1-0 emission and exhibits an overall velocity gradient consistent with the gradient at 7mm. While there is spatial overlap between the 3mm and 7mm transitions, the 3mm emission, on average, lies at larger projected distances from Source I (~44 AU compared with ~35 AU for 7mm). The spatial overlap between the v=1, J=1-0 and J=2-1 transitions is suggestive of a range of temperatures and densities where physical conditions are favorable for both transitions of a same vibrational state. However, the observed spatial offset between the bulk of emission at 3mm and 7mm possibly indicates different ranges of temperatures and densities for optimal excitation of the masers. We discuss different maser pumping models that may explain the observed offset. We interpret the 3mm and 7mm masers as being part of a single wide-angle outflow arising from the surface of an edge-on disk rotating about a northeast-southwest axis, with a continuous velocity gradient indicative of differential rotation consistent with a Keplerian profile in a high-mass proto-binary.

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

Advances in astronomy are intimately linked to advances in digital signal processing (DSP). This special issue is focused upon advances in DSP within radio astronomy. The trend within that community is to use off-the-shelf digital hardware where possible and leverage advances in high performance computing. In particular, graphics processing units (GPUs) and field programmable gate arrays (FPGAs) are being used in place of application-specific circuits (ASICs); high-speed Ethernet and Infiniband are being used for interconnect in place of custom backplanes. Further, to lower hurdles in digital engineering, communities have designed and released general-purpose FPGA-based DSP systems, such as the CASPER ROACH board, ASTRON Uniboard and CSIRO Redback board. In this introductory article, we give a brief historical overview, a summary of recent trends, and provide an outlook on future directions.

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.

The Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) has been working for a decade to reduce the time and cost of designing, building and deploying new digital radio-astronomy instruments. Today, CASPER open-source technology powers over 45 scientific instruments worldwide, and is used by scientists and engineers at dozens of academic institutions. In this paper we catalog the current offerings of the CASPER collaboration, and instruments past and present built by CASPER users and developers. We describe the ongoing state of software development, as CASPER looks to support a broader range of programming environments and hardware and ensure compatibility with the latest vendor tools.

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.

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.

The birth of the first luminous sources and the ensuing epoch of reionization are best studied via the redshifted 21-cm emission line, the signature of the first two imprinting the last. In this work we present a fully-Bayesian method, \textsc{hibayes}, for extracting the faint, global (sky-averaged) 21-cm signal from the much brighter foreground emission. We show that a simplified (but plausible), Gaussian model of the 21-cm emission from the Cosmic Dawn epoch ($15 \lesssim z \lesssim 30$), parameterized by an amplitude $A_{\rm HI}$, a frequency peak $\nu_{\rm HI}$ and a width $\sigma_{\rm HI}$, can be extracted even in the presence of a structured foreground frequency spectrum (parameterized as a $7^{\rm th}$-order polynomial), provided sufficient signal-to-noise (400~hours of observation with a single dipole). We apply our method to an early, 19-minute long observation from the Large aperture Experiment to detect the Dark Ages, constraining the 21-cm signal amplitude and width to be $-890 < A_{\rm HI} < 0$ mK and $\sigma_{\rm HI} > 6.5$ MHz (corresponding to $\Delta z > 1.9$ at redshift $z \simeq 20$) respectively at the 95-per-cent confidence level in the range $13.2 < z < 27.4$ ($100 > \nu > 50$ MHz).

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.

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.

The magnitude of the real-time digital signal processing challenge attached to large radio astronomical antenna arrays motivates use of high performance computing (HPC) systems. The need for high power efficiency (performance per watt) at remote observatory sites parallels that in HPC broadly, where efficiency is an emerging critical metric. We investigate how the performance per watt of graphics processing units (GPUs) is affected by temperature, core clock frequency and voltage. Our results highlight how the underlying physical processes that govern transistor operation affect power efficiency. In particular, we show experimentally that GPU power consumption grows non-linearly with both temperature and supply voltage, as predicted by physical transistor models. We show lowering GPU supply voltage and increasing clock frequency while maintaining a low die temperature increases the power efficiency of an NVIDIA K20 GPU by up to 37-48% over default settings when running xGPU, a compute-bound code used in radio astronomy. We discuss how temperature-aware power models could be used to reduce power consumption for future HPC installations. Automatic temperature-aware and application-dependent voltage and frequency scaling (T-DVFS and A-DVFS) may provide a mechanism to achieve better power efficiency for a wider range of codes running on GPUs

The FITS (Flexible Image Transport System) data format has been the de facto data format for astronomy-related data products since its inception in the late 1970s. While the FITS file format is widely supported, it lacks many of the features of more modern data serialization, such as the Hierarchical Data Format (HDF5). The HDF5 file format offers considerable advantages over FITS, such as improved I/O speed and compression, but has yet to gain widespread adoption within astronomy. One of the major holdbacks is that HDF5 is not well supported by data reduction software packages and image viewers. Here, we present a comparison of FITS and HDF5 as a format for storage of astronomy datasets. We show that the underlying data model of FITS can be ported to HDF5 in a straightforward manner, and that by doing so the advantages of the HDF5 file format can be leveraged immediately. In addition, we present a software tool, fits2hdf, for converting between FITS and a new `HDFITS' format, where data are stored in HDF5 in a FITS-like manner. We show that HDFITS allows faster reading of data (up to 100x of FITS in some use cases), and improved compression (higher compression ratios and higher throughput). Finally, we show that by only changing the import lines in Python-based FITS utilities, HDFITS formatted data can be presented transparently as an in-memory FITS equivalent.

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.

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

A "large-N" correlator that makes use of Field Programmable Gate Arrays and Graphics Processing Units has been deployed as the digital signal processing system for the Long Wavelength Array station at Owens Valley Radio Observatory (LWA-OV), to enable the Large Aperture Experiment to Detect the Dark Ages (LEDA). The system samples a ~100MHz baseband and processes signals from 512 antennas (256 dual polarization) over a ~58MHz instantaneous sub-band, achieving 16.8Tops/s and 0.236 Tbit/s throughput in a 9kW envelope and single rack footprint. The output data rate is 260MB/s for 9 second time averaging of cross-power and 1 second averaging of total-power data. At deployment, the LWA-OV correlator was the largest in production in terms of N and is the third largest in terms of complex multiply accumulations, after the Very Large Array and Atacama Large Millimeter Array. The correlator's comparatively fast development time and low cost establish a practical foundation for the scalability of a modular, heterogeneous, computing architecture.

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.

The most promising near-term observable of the cosmic dark age prior to widespread reionization (z~15-200) is the sky-averaged \lambda 21 cm background arising from hydrogen in the intergalactic medium. Though an individual antenna could in principle detect the line signature, data analysis must separate foregrounds that are orders of magnitude brighter than the \lambda 21 cm background (but that are anticipated to vary monotonically and gradually with frequency). Using more physically motivated models for foregrounds than in previous studies, we show that the intrinsic "spectral smoothness" of the foregrounds is likely not a concern, and that data analysis for an ideal antenna should be able to detect the \lambda 21 cm signal after deprojecting a ~5th order polynomial in log(\nu). However, we find that the foreground signal is corrupted by the frequency-dependent response of a real antenna. The frequency dependence complicates modeling of foregrounds commonly based on the assumption of spectral smoothness. Much of our study focuses on the Large-aperture Experiment to detect the Dark Age (LEDA), which combines both radiometric and interferometric measurements. We show that statistical uncertainty remaining after fitting antenna gain patterns to interferometric measurements does not compromise extraction of the \lambda 21 cm signal for a range of cosmological models after fitting a 7th order polynomial to radiometric data. Our results generalize to most efforts to measure the sky-averaged spectrum.

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.

The FITS (Flexible Image Transport System) data format was developed in the late 1970s for storage and exchange of astronomy-related image data. Since then, it has become a standard file format not only for images, but also for radio interferometer data (e.g. UVFITS, FITS-IDI). But is FITS the right format for next-generation telescopes to adopt? The newer Hierarchical Data Format (HDF5) file format offers considerable advantages over FITS, but has yet to gain widespread adoption within radio astronomy. One of the major holdbacks is that HDF5 is not well supported by data reduction software packages. Here, we present a comparison of FITS, HDF5, and the MeasurementSet (MS) format for storage of interferometric data. In addition, we present a tool for converting between formats. We show that the underlying data model of FITS can be ported to HDF5, a first step toward achieving wider HDF5 support.

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.

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.

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.

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.

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.

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.

Radio astronomical imaging arrays comprising large numbers of antennas, O(10^2-10^3) have posed a signal processing challenge because of the required O(N^2) cross correlation of signals from each antenna and requisite signal routing. This motivated the implementation of a Packetized Correlator architecture that applies Field Programmable Gate Arrays (FPGAs) to the O(N) "F-stage" transforming time domain to frequency domain data, and Graphics Processing Units (GPUs) to the O(N^2) "X-stage" performing an outer product among spectra for each antenna. The design is readily scalable to at least O(10^3) antennas. Fringes, visibility amplitudes and sky image results obtained during field testing are presented.

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.

We report a new geometric maser distance estimate to the active galaxy NGC 4258. The data for the new model are maser line-of-sight velocities and sky positions from 18 epochs of Very Long Baseline Interferometry observations, and line-of-sight accelerations measured from a 10-year monitoring program of the 22 GHz maser emission of NGC 4258. The new model includes both disk warping and confocal elliptical maser orbits with differential precession. The distance to NGC 4258 is 7.60 +/- 0.17 +/- 0.15 Mpc, a 3% uncertainty including formal fitting and systematic terms. The resulting Hubble Constant, based on the use of the Cepheid Variables in NGC 4258 to recalibrate the Cepheid distance scale (Riess et al. 2011), is H_0 = 72.0 +/- 3.0 km/s/Mpc.

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

In Orion BN/KL, proper motions of 7 mm vibrationally-excited SiO masers trace rotation of a nearly edge-on disk and a bipolar wide-angle outflow 10-100 AU from radio Source I, a binary young stellar object (YSO) of ~20 Msun. Here we map ground-state 7 mm SiO emission with the Very Large Array and track proper motions over 9 years. The innermost and strongest emission lies in two extended arcs bracketing Source I. The proper motions trace a northeast-southwest bipolar outflow 100-1000 AU from Source I with a median 3D motion of ~18 km/s. An overlying distribution of 1.3 cm H2O masers betrays similar flow characteristics. Gas dynamics and emission morphology traced by the masers suggest the presence of a magnetocentrifugal disk-wind. Reinforcing evidence lies in the colinearity of the flow, apparent rotation across the flow parallel to the disk rotation, and recollimation that narrows the flow opening angle ~120 AU downstream. The arcs of ground-state SiO emission may mark the transition point to a shocked super-Alfvenic outflow.

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 report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

There are few data available with which to constrain the thermal history of the intergalactic medium (IGM) following global recombination. Thus far, most constraints flow from analyses of the Cosmic Microwave Background and optical spectroscopy along a few lines of sight. However, direct study of the IGM in emission or absorption against the CMB via the 1S hyperfine transition of Hydrogen would enable broad characterization thermal history and source populations. New generations of radio arrays are in development to measure this line signature. Bright foreground emission and the complexity of instrument calibration models are significant hurdles. How to optimize these is uncertain, resulting in a diversity in approaches. We discuss recent limits on line brightness, array efforts including the new Large Aperture Experiment to Detect the Dark Ages (LEDA), and the next generation Hydrogen Reionization Array (HERA) concept.

We study the optical properties of the host galaxies of nuclear 22 GHz (lambda=1.35 cm) water masers. To do so, we cross-match the galaxy sample surveyed for water maser emission (123 detections and 3806 non-detections) with the SDSS low-redshift galaxy sample (z<0.05). Out of 1636 galaxies with SDSS photometry, we identify 48 detections; out of the 1063 galaxies that also have SDSS spectroscopy, we identify 33 detections. We find that maser detection rate is higher at higher optical luminosity (M_B), larger velocity dispersion ($\sigma$), and higher [OIII]5007 luminosity, with [OIII]5007 being the dominant factor. These detection rates are essentially the result of the correlations of isotropic maser luminosity with all three of these variables. These correlations are natural if maser strength increases with central black hole mass and the level of AGN activity. We also find that the detection rate is higher in galaxies with higher extinction. Based on these results, we propose that maser surveys seeking to efficiently find masers should rank AGN targets by extinction-corrected [OIII]5007 flux when available. This prioritization would improve maser detection efficiency, from an overall ~ 3% without pre-selection to ~ 16% for the strongest intrinsic [OIII]5007 emitters, by a factor of ~ 5.

We present a highly parallel implementation of the cross-correlation of time-series data using graphics processing units (GPUs), which is scalable to hundreds of independent inputs and suitable for the processing of signals from "Large-N" arrays of many radio antennas. The computational part of the algorithm, the X-engine, is implementated efficiently on Nvidia's Fermi architecture, sustaining up to 79% of the peak single precision floating-point throughput. We compare performance obtained for hardware- and software-managed caches, observing significantly better performance for the latter. The high performance reported involves use of a multi-level data tiling strategy in memory and use of a pipelined algorithm with simultaneous computation and transfer of data from host to device memory. The speed of code development, flexibility, and low cost of the GPU implementations compared to ASIC and FPGA implementations have the potential to greatly shorten the cycle of correlator development and deployment, for cases where some power consumption penalty can be tolerated.

Using the Expanded Very Large Array, we have mapped the vicinity of the Orion BN/KL Hot Core with sub-arcsecond angular resolution in seven metastable inversion transitions of ammonia: (J,K)=(6,6) to (12,12). This emission comes from levels up to 1500 K above the ground state, enabling identification of source(s) responsible for heating the region. We used this multi-transition dataset to produce images of the rotational/kinetic temperature and the column density of ammonia for ortho and para species separately and on a position-by-position basis. We find rotational temperature and column density in the range 160-490 K and (1-4)x10^17 cm^-2, respectively. Our spatially-resolved images show that the highest (column) density and hottest gas is found in a northeast-southwest elongated ridge to the southeast of Source I. We have also measured the ortho-para ratio of ammonia, estimated to vary in the range 0.9-1.6. Enhancement of ortho with respect to para and the offset of hot ammonia emission peaks from known (proto)stellar sources provide evidence that the ammonia molecules have been released from dust grains into the gas-phase through the passage of shocks and not by stellar radiation. We propose that the combined effect of Source I's proper motion and its low-velocity outflow impinging on a pre-existing dense medium is responsible for the excitation of ammonia and the Orion Hot Core. Finally, we found for the first time evidence of a slow (5 km/s) and compact (1000 AU) outflow towards IRc7.

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.

We present new 7mm continuum observations of Orion BN/KL with the VLA. We resolve the emission from the protostar radio Source I and BN at several epochs. Source I is highly elongated NW-SE, and remarkably stable in flux density, position angle, and overall morphology over nearly a decade. This favors the extended emission component arising from an ionized disk rather than a jet. We have measured the proper motions of Source I and BN for the first time at 43 GHz. We confirm that both sources are moving at high speed (12 and 26 km/s, respectively) approximately in opposite directions, as previously inferred from measurements at lower frequencies. We discuss dynamical scenarios that can explain the large motions of both BN and Source I and the presence of disks around both. Our new measurements support the hypothesis that a close (~50 AU) dynamical interaction occurred around 500 years ago between Source I and BN as proposed by Gomez et al. From the dynamics of encounter we argue that Source I today is likely to be a binary with a total mass on the order of 20 Msun, and that it probably existed as a softer binary before the close encounter. This enables preservation of the original accretion disk, though truncated to its present radius of ~50 AU. N-body numerical simulations show that the dynamical interaction between a binary of 20 Msun total mass (I) and a single star of 10 Msun mass (BN) may lead to the ejection of both and binary hardening. The gravitational energy released in the process would be large enough to power the wide-angle flow traced by H2 and CO emission in the BN/KL nebula. Assuming the proposed dynamical history is correct, the smaller mass for Source I recently estimated from SiO maser dynamics (>7 Msun) by Matthews et al., suggests that non-gravitational forces (e.g. magnetic) must play an important role in the circumstellar gas dynamics.

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.

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.

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.

Correlations between X-ray and water maser emission in AGN have been recently reported. However, the lack of systematic studies affects the confidence level of these results. In the following, we introduce a project aimed at studying all the water maser sources believed to be associated with AGN activity through X-ray data obtained with the XRT and BAT instruments on-board the Swift satellite. Preliminary results of this work indicate a promising rate of XRT detections allowing us to refine follow-up observing strategies focused on investigating the nuclei of individual galaxies and deriving, on statistical basis, the main characteristics of water maser hosts. In addition, a cross-correlation between our sample and the BAT 22-months all-sky survey provides an exceptionally high detection rate at hard X-ray energies when compared to other AGN-related catalogs.

We present multi-epoch VLBA imaging of the 28SiO v=1 & v=2, J=1-0 maser emission toward the massive YSO Orion Source I. Both SiO transitions were observed simultaneously with an angular resolution of ~0.5 mas (~0.2 AU for d=414 pc). Here we explore the global properties and kinematics of the emission through two 19-epoch animated movies spanning 21 months (2001 March 19 to 2002 December 10). These movies provide the most detailed view to date of the dynamics and temporal evolution of molecular material within ~20-100 AU of a massive (~>8M_sun) YSO. The bulk of the SiO masers surrounding Source I lie in an X-shaped locus; emission in the South/East arms is predominantly blueshifted and in the North and West is predominantly redshifted. In addition, bridges of intermediate-velocity emission connect the red and blue sides of the emission distribution. We have measured proper motions of >1000 maser features and find a combination of radially outward migrations along the four arms and motions tangent to the bridges. We interpret the SiO masers as arising from a wide-angle bipolar wind emanating from a rotating, edge-on disk. The detection of maser features along extended, curved filaments suggests that magnetic fields may play a role in launching and/or shaping the wind. Our observations appear to support a picture in which stars with M ~>8 M_sun form via disk-mediated accretion. However, we cannot rule out that the Source I disk may have been formed or altered following a close encounter. (Abridged).

Non-zero neutrino mass would affect the evolution of the Universe in observable ways, and a strong constraint on the mass can be achieved using combinations of cosmological data sets. We focus on the power spectrum of cosmic microwave background (CMB) anisotropies, the Hubble constant H_0, and the length scale for baryon acoustic oscillations (BAO) to investigate the constraint on the neutrino mass, m_nu. We analyze data from multiple existing CMB studies (WMAP5, ACBAR, CBI, BOOMERANG, and QUAD), recent measurement of H_0 (SHOES), with about two times lower uncertainty (5%) than previous estimates, and recent treatments of BAO from the Sloan Digital Sky Survey (SDSS). We obtained an upper limit of m_nu < 0.2eV (95% C.L.), for a flat LambdaCDM model. This is a 40% reduction in the limit derived from previous H_0 estimates and one-third lower than can be achieved with extant CMB and BAO data. We also analyze the impact of smaller uncertainty on measurements of H_0 as may be anticipated in the near term, in combination with CMB data from the Planck mission, and BAO data from the SDSS/BOSS program. We demonstrate the possibility of a 5 sigma detection for a fiducial neutrino mass of 0.1eV or a 95% upper limit of 0.04eV for a fiducial of m_nu = 0eV. These constraints are about 50% better than those achieved without external constraint. We further investigate the impact on modeling where the dark-energy equation of state is constant but not necessarily -1, or where a non-flat universe is allowed. In these cases, the next-generation accuracies of Planck, BOSS, and 1% measurement of H_0 would all be required to obtain the limit m_nu < 0.05 - 0.06eV (95% C.L.) for the fiducial of m_nu = 0eV. The independence of systematics argues for pursuit of both BAO and H_0 measurements.

Based on spectroscopic signatures, about one-third of known H$_2$O maser sources in active galactic nuclei (AGN) are believed to arise in highly inclined accretion disks around central engines. These "disk maser candidates" are of interest primarily because angular structure and rotation curves can be resolved with interferometers, enabling dynamical study. We identify five new disk maser candidates in studies with the Green Bank Telescope, bringing the total number published to 30. We discovered two (NGC1320, NGC17) in a survey of 40 inclined active galaxies (v_{sys}< 20000 kms^{-1}). The remaining three disk maser candidates were identified in monitoring of known sources: NGC449, NGC2979, NGC3735. We also confirm a previously marginal case in UGC4203. For the disk maser candidates reported here, inferred rotation speeds are 130-500 kms^{-1}. Monitoring of three more rapidly rotating candidate disks (CG211, NGC6264, VV340A) has enabled measurement of likely orbital centripetal acceleration, and estimation of central masses (2-7x10^7 M_\odot) and mean disk radii (0.2-0.4pc). Accelerations may ultimately permit estimation of distances when combined with interferometer data. This is notable because the three AGN are relatively distant (10000<v_{sys}<15000 kms^{-1}). As signposts of highly inclined geometries at galactocentric radii of \sim0.1-1pc, disk masers also provide robust orientation references that allow analysis of (mis)alignment between AGN and surrounding galactic stellar disks, even without interferometric mapping. We find no preference among published disk maser candidates to lie in high-inclination galaxies, providing independent support for conclusions that central engines and galactic plane orientations are not correlated. (ABRIDGED)

Modern graphics processing units (GPUs) are inexpensive commodity hardware that offer Tflop/s theoretical computing capacity. GPUs are well suited to many compute-intensive tasks including digital signal processing. We describe the implementation and performance of a GPU-based digital correlator for radio astronomy. The correlator is implemented using the NVIDIA CUDA development environment. We evaluate three design options on two generations of NVIDIA hardware. The different designs utilize the internal registers, shared memory and multiprocessors in different ways. We find that optimal performance is achieved with the design that minimizes global memory reads on recent generations of hardware. The GPU-based correlator outperforms a single-threaded CPU equivalent by a factor of 60 for a 32 antenna array, and runs on commodity PC hardware. The extra compute capability provided by the GPU maximises the correlation capability of a PC while retaining the fast development time associated with using standard hardware, networking and programming languages. In this way, a GPU-based correlation system represents a middle ground in design space between high performance, custom built hardware and pure CPU-based software correlation. The correlator was deployed at the Murchison Widefield Array 32 antenna prototype system where it ran in real-time for extended periods. We briefly describe the data capture, streaming and correlation system for the prototype array.

This is the first of two papers reporting measurements from a program to determine the Hubble constant to 5% precision from a refurbished distance ladder. We present new observations of 110 Cepheid variables in the host galaxies of two recent Type Ia supernovae (SNe Ia), NGC 1309 and NGC 3021, using the Advanced Camera for Surveys on the Hubble Space Telescope (HST). We also present new observations of the hosts previously observed with HST whose SNe Ia provide the most precise luminosity calibrations: SN 1994ae in NGC 3370, SN 1998aq in NGC 3982, SN 1990N in NGC 4639, and SN 1981B in NGC 4536, as well as the maser host, NGC 4258. Increasing the interval between observations enabled the discovery of new, longer-period Cepheids, including 57 with P>60 days, which extend these period-luminosity (PL) relations. We present 93 measurements of the metallicity parameter, 12 + log[O/H], measured from HII regions in the vicinity of the Cepheids and show these are consistent with solar metallicity. We find the slope of the seven dereddened PL relations to be consistent with that of the Large Magellanic Cloud Cepheids and with parallax measurements of Galactic Cepheids, and we address the implications for the Hubble constant. We also present multi-band light curves of SN 2002fk (in NGC 1309) and SN 1995al (in NGC 3021) which may be used to calibrate their luminosities. In the second paper we present observations of the Cepheids in the H-band obtained with the Near Infrared Camera and Multi-Object Spectrometer on HST, further mitigating systematic errors along the distance ladder resulting from dust and chemical variations. The quality and homogeneity of these SN and Cepheid data provide the basis for a more precise determination of the Hubble constant.

We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions (J=1-0) from three SiO isotopologues towards Orion BN/KL: 28SiO v=0,1,2; and 29SiO and 30SiO v=0. For the first time, we have mapped the 29SiO and 30SiO J=1-0 emission, established the maser nature of the emission, and confirmed association with the deeply embedded high-mass young stellar object commonly denoted radio Source I. The 28SiO v=0 maser emission shows a bipolar structure that extends over ~700 AU along a northeast-southwest axis, and we propose that it traces a bipolar outflow driven by Source I. The high-brightness isotopic SiO maser emission imaged with a ~0.2 arcsec resolution has a more compact distribution, generally similar to that of the 28SiO v=1,2 emission, and it probably traces bulk gas flows in a region of diameter <100 AU centered on Source I. On small scales of <10 AU, however, compact 29SiO/30SiO v=0 and 28SiO v=1,2 emission features may be offset from one another in position and line-of-sight velocity. From a radiative transfer analysis based on a large velocity gradient (LVG) pumping model, we derive similar temperatures and densities for the optimum excitation of both 29SiO/30SiO v=0 and 28SiO v=1,2 masers, significantly higher than required for 28SiO v=0 maser excitation. In order to account for the small-scale differences among the isotopologues (v=0) and the main species (v=1,2), follow-up radiative transfer modeling that incorporates non-local line overlap among transitions of all SiO isotopic species may be required.

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.

Joint analysis of Cosmic Microwave Background, Baryon Acoustic Oscillation, and supernova data has enabled precision estimation of cosmological parameters. New programs will push to 1% uncertainty in the dark energy equation of state and tightened constraint on curvature, requiring close attention to systematics. Direct 1% measurement of the Hubble constant (H0) would provide a new constraint. It can be obtained without overlapping systematics directly from recessional velocities and geometric distance estimates for galaxies via the mapping of water maser emission that traces the accretion disks of nuclear black holes. We identify redshifts 0.02<z<0.06 as best for small samples, e.g., 10 widely distributed galaxies, each with 3% distance uncertainty. Knowledge of peculiar radial motion is also required. Mapping requires very long baseline interferometry (VLBI) with the finest angular resolution, sensitivity to individual lines of a few mJy-km/s, and baselines that can detect a complex of ~10 mJy lines (peak) in < 1 min. For 2010-2020, large ground apertures (50-100m diameter) augmenting the VLBA are critical, such as EVLA, GBT, Effelsberg, and the Large Millimeter Telescope, for which we propose a 22 GHz receiver and VLBI instrumentation. A space-VLBI aperture may be required, thus motivating US participation in the Japanese VSOP-2 mission (launch c.2013). This will provide 3-4x longer baselines and ~5x improvement in distance uncertainty. There are now 5 good targets at z>0.02, out of ~100 known masers. A single-dish discovery survey of >10,000 nuclei (>2500 hours on the GBT) would build a sample of tens of potential distance anchors. Beyond 2020, a high-frequency SKA could provide larger maser samples, enabling estimation of H0 from individually less accurate distances, and possibly without the need for peculiar motion corrections.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

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 Megamaser Cosmology Project (MCP) seeks to measure the Hubble Constant (Ho) in order to improve the extragalactic distance scale and constrain the nature of dark energy. We are searching for sources of water maser emission from AGN with sub-pc accretion disks, as in NGC 4258, and following up these discoveries with Very Long Baseline Interferometric (VLBI) imaging and spectral monitoring. Here we present a VLBI map of the water masers toward UGC 3789, a galaxy well into the Hubble Flow. We have observed masers moving at rotational speeds up to 800 km/s at radii as small as 0.08 pc. Our map reveals masers in a nearly edge-on disk in Keplerian rotation about a 10^7 Msun supermassive black hole. When combined with centripetal accelerations, obtained by observing spectral drifts of maser features (to be presented in Paper II), the UGC 3789 masers may provide an accurate determination of Ho, independent of luminosities and metallicity and extinction corrections.

We have surveyed molecular line emission from Orion BN/KL from 42.3 to 43.6 GHz with the Green Bank Telescope. Sixty-seven lines were identified and ascribed to 13 different molecular species. The spectrum at 7 mm is dominated by SiO, SO2, CH3OCH3, and C2H5CN. Five transitions have been detected from the SiO isotopologues 28SiO, 29SiO, and 30SiO. We report here for the first time the spectra of the 29SiO and 30SiO v=0 J=1-0 emission in Orion BN/KL, and we show that they have double-peaked profiles with velocity extents similar to the main isotopologue. The main motivation for the survey was the search of high-velocity (100-1000 km/s) outflows in the BN/KL region as traced by SiO Doppler components. Some of the unidentified lines in principle could be high-velocity SiO features, but without imaging data their location cannot be established. Wings of emission are present in the v=0 28SiO, 29SiO and 30SiO profiles, and we suggest that the v=0 emission from the three isotopologues might trace a moderately high-velocity (~30-50 km/s) component of the flows around the high-mass protostar Source I in the Orion BN/KL region. We also confirm the 7 mm detection of a complex oxygen-bearing species, acetone (CH3COCH3), which has been recently observed towards the hot core at 3 mm, and we have found further indications of the presence of long cyanopolyynes (HC5N and HC7N) in the quiescent cold gas of the extended ridge.

We report Suzaku observations of the active, Compton-thick Circinus galaxy. Observations were obtained with both the X-ray Imaging spectrometer (XIS) and the Hard X-ray Detector (HXD). Below 10 keV, the nuclear spectrum is dominated by radiation reflected from cold dense gas of high column density, while above 13 keV the radiation is directly transmitted nuclear emission seen through a column density of $~ 4 x 10^{24} cm^-2. In the 0.2--10 keV band, the XIS spectrum is contaminated at 5% level by the brightest off-nuclear source in Circinus (CG X-1), but drops to 1% in the 5-10 keV and is negligible at higher energies. We find no significant evidence for variability in the hard (>12 keV) emission. The Circinus is marginally detected with the HXD/GSO in the 50--100 keV band at 2.5\sigma level. We model the 3-70 keV band XIS+PIN spectra with a four components: the Compton transmitted nuclear emission, the reflected nuclear emission, a soft power law (representing a combination of scattered nuclear emission, extended emission and contamination by sources in the galaxy below a few keV). The hard nuclear power-law is found to have a photon index Gamma_h ~= 1.6, very similar to the soft power-law. The high energy cut-off is E_C ~= 49 keV. These results agree with those from BeppoSax. An extrapolation of this model up to the GSO band shows good agreement with the GSO spectrum and supports our detection of the Circinus up to ~= 100keV.

Of 104 AGN known to exhibit H$_2$O maser emission, X-ray data that enable estimation of column densities, or lower limits, are available for 42. Contributing to this, we report analysis of new and archival X-ray data for 8 galaxies and collation of values for three more. Maser emission is indicative of large columns of cold gas, and in five of the eight new cases, maser spectra point toward origins in accretion disks viewed close-to edge-on (a.k.a. "disk-maser" systems). In these, we detect hard continuum and Fe K$\alpha$ emission with equivalent widths on the order of 1 keV, which is consistent with Compton reflection, fluorescence by cold material, and obscuring columns $\ga 10^{24}$ cm$^{-2}$. Reviewing the full sample of 42, 95% exhibit N$_{\rm H} >10^{23}$ cm$^{-2}$ and 60% exhibit N$_{\rm H} >10^{24}$ cm$^{-2}$. Half of these are now recognized to be disk masers (up from 13); in this sub-sample, which is likely to be more homogeneous vis-\'a-vis the origin of maser emission, 76% exhibit N$_{\rm H} >10^{24}$ cm$^{-2}$. The probability of a common parent distribution of columns for disk-masers and other AGN masers is <3%. Because ground-based surveys of AGN to detect new disk masers are relatively unbiased with respect to X-ray brightness and comparatively inexpensive, they may also be efficient guides for the sensitive pointed X-ray observations required to identify Compton-thick objects outside of shallow surveys.

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.

We present a Very Long Baseline Interferometry image of the water maser emission in the nuclear region of NGC3393. The maser emission has a linear distribution oriented at a position angle of $\sim -34\degr$, perpendicular to both the kpc-scale radio jet and the axis of the narrow line region. The position-velocity diagram displays a red-blue asymmetry about the systemic velocity and the estimated dynamical center, and is thus consistent with rotation. Assuming Keplerian rotation in an edge-on disk, we obtain an enclosed mass of $(3.1\pm 0.2) \times 10^7 M_{\sun}$ within $0.36\pm 0.02$ pc ($1.48\pm 0.06$ mas), which corresponds to a mean mass density of $\sim10^{8.2} M_{\sun}$ pc$^{-3}$. We also report the measurement with the Green Bank Telescope of a velocity drift, a manifestation of centripetal acceleration within the disk, of $5\pm 1$ km s $^{-1}$ yr$^{-1}$ in the $\sim3880$ km s$^{-1}$ maser feature, which is most likely located along the line of sight to the dynamical center of the system. From the acceleration of this feature, we estimate a disk radius of $0.17\pm 0.02$ pc, which is smaller than the inner disk radius ($0.36\pm 0.02$ pc) of emission that occurs along the midline (i.e., the line of nodes). The emission along the line of sight to the dynamical center evidently occurs much closer to the center than the emission from the disk midline, contrary to the situation in the archetypal maser systems NGC4258 and NGC1068. The outer radius of the disk as traced by the masers along the midline is about 1.5 pc.

We present a new XMM-Newton spectrum of the Seyfert 2 nucleus of IC 2560, which hosts H$_{2}$O maser emission from an inclined Keplerian accretion disk. The X-ray spectrum shows soft excess due to multi-temperature ionized plasma, a hard continuum and strong emission features, from Mg, Si, S, Ca, Fe and Ni, mainly due to fluorescence. It is consistent with reflection of the continuum from a mostly neutral medium and obscuration due to a high column density, $>$ 10$^{24}$ cm$^{-2}$. The amplitude of the reflected component may exceed 10% of the central unobscured luminosity. This is higher than the reflected fraction, of a few percent, observed in other Seyfert 2 sources like NGC 4945. We observe an emission line at 6.7 keV, possibly due to FeXXV, undetected in previous Chandra observations. The absorption column density associated with this line is less than 10$^{23}$ cm$^{-2}$, lower than the obscuration of the central source. We hypothesize that this highly ionized Fe line emission originates in warm gas, also responsible for a scattered component of continuum emission that may dominate the spectrum between 1 and 3 keV. We compare X-ray and maser emission characteristics of IC 2560 and other AGN that exhibit water maser emission originating in disk structures around central engines. The temperature for the region of the disk associated with maser action is consistent with the expected 400-1000K range. The clumpiness of disk structures (inferred from the maser distribution) may depend on the unobscured luminosities of the central engines.

We report measurements of centripetal accelerations of maser spectral components of NGC 4258 for 51 epochs spanning 1994 to 2004. This is the second paper of a series, in which the goal is determination of a new geometric maser distance to NGC 4258 accurate to possibly ~3%. We measure accelerations using a formal analysis method that involves simultaneous decomposition of maser spectra for all epochs into multiple, Gaussian components. Components are coupled between epochs by linear drifts (accelerations) from their centroid velocities at a reference epoch. For high-velocity emission, accelerations lie in the range -0.7 to +0.7 km/s/yr indicating an origin within 13 degrees of the disk midline (the perpendicular to the line-of-sight to the black hole). Comparison of high-velocity emission projected positions in VLBI images, with those derived from acceleration data, provides evidence that masers trace real gas dynamics. High-velocity emission accelerations do not support a model of trailing shocks associated with spiral arms in the disk. However, we find strengthened evidence for spatial periodicity in high-velocity emission, of wavelength 0.75 mas. This supports suggestions of spiral structure due to density waves in the nuclear accretion disk of an active galaxy. Accelerations of low-velocity (systemic) emission lie in the range 7.7 to 8.9 km/s/yr, consistent with emission originating from a concavity where the thin, warped disk is tangent to the line-of-sight. A trend in accelerations of low-velocity emission as a function of Doppler velocity may be associated with disk geometry and orientation, or with the presence of spiral structure.

Galactic nuclei are well known sources of OH and H2O maser emission. It appears that intense star formation in ultra-luminous infrared galaxies drives most OH sources. In contrast, nuclear activity appears to drive most H2O sources. When H2O emission originates in accretion disk structures, constrained geometry and dynamics enable robust interpretation of spectroscopic and imaging data. The principal science includes study of AGN geometry at parsec and sub-parsec radii and measurement of geometric distances in the Hubble Flow. New high accuracy estimates of the Hubble constant, "Ho," obtained from maser distances may enable new substantively improved constraints on fundamental cosmological parameters (e.g., dark energy).

A wealth of new information about the structure of the maser disk in NGC 4258 has been obtained from a series of 18 VLBA observations spanning three years, as well as from 32 additional epochs of spectral monitoring data from 1994 to the present, acquired with the VLA, Effelsberg, and GBT. The warp of the disk has been defined precisely. The thickness of the maser disk has been measured to be 12 microarcseconds (FWHM), which is slightly smaller than previously quoted upper limits. Under the assumption that the masers trace the true vertical distribution of material in the disk, from the condition of hydrostatic equilibrium the sound speed is 1.5 km/s, corresponding to a thermal temperature of 600K. The accelerations of the high velocity maser components have been accurately measured for many features on both the blue and red side of the spectrum. The azimuthal offsets of these masers from the midline (the line through the disk in the plane of the sky) and derived projected offsets from the midline based on the warp model correspond well with the measured offsets. This result suggests that the masers are well described as discrete clumps of masing gas, which accurately trace the Keplerian motion of the disk. However, we have continued to search for evidence of apparent motions caused by ``phase effects.'' This work provides the foundation for refining the estimate of the distance to NGC 4258 through measurements of feature acceleration and proper motion. The refined estimate of this distance is expected to be announced in the near future.

The relatively nearby (distance=24.1 Mpc) elliptical galaxy NGC821 hosts an extreme example of a quiescent central massive black hole, for which deep Chandra observations revealed a nuclear source for the first time (with L(2-10 keV)/L_Edd ~ 3.6X10^{-8}). We present here a multiwavelength study of this nucleus, including VLA observations that detect a radio counterpart to the Chandra nuclear source at 1.4 GHz, with a flux density of 127 $\mu$m and possibly a flat spectral shape; we also consider new Spitzer IRAC observations and archival HST images. With these data we discuss possible scenarios for the accretion modalities of the sole material that is available for fuelling, i.e., the stellar mass losses steadily replenishing the circumnuclear region. The final stages of accretion could be radiatively inefficient and coupled to a compact nuclear jet/outfow. The stellar mass losses could instead end up in a standard disc only if a Compton-thick AGN is present. Two extended sources detected by Chandra close to the nucleus could be due to several unresolved knots in a jet. If a jet is present, though, its kinetic energy would be only a very small fraction of the energy associated with the rest mass of the material being accreted. Starformation close to the nucleus is not shown by the available data. Deeper NICMOS, radio and far-IR observations should further constrain the accretion process.

A comprehensive picture of high-mass star formation has remained elusive, in part because examples of high-mass YSOs tend to be relatively distant, deeply embedded, and confused with other emission sources. These factors have impeded dynamical investigations within tens of AU of high-mass YSOs--scales that are critical for probing the interfaces where outflows from accretion disks are launched and collimated. Using observations of SiO masers obtained with the VLA and the VLBA, the KaLYPSO project is overcoming these limitations by mapping the structure and dynamical/temporal evolution of the material 10-1000 AU from the nearest high-mass YSO: Radio Source I in the Orion BN/KL region. Our data include ~40 epochs of VLBA observations over a several-year period, allowing us to track the proper motions of individual SiO maser spots and to monitor changes in the physical conditions of the emitting material with time. Ultimately these data will provide 3-D maps of the outflow structure over approximately 30% of the outflow crossing time. Here we summarize recent results from the KaLYPSO project, including evidence that high-mass star formation is occurring via disk-mediated accretion.

We have imaged the enigmatic radio source-I (Orion-I) in the Orion-KL nebula with the VLA at 43 GHz with 34 mas angular resolution. The continuum emission is highly elongated and is consistent with that expected from a nearly edge-on disk. The high brightness and lack of strong molecular lines from Orion-I can be used to argue against emission from dust. Collisional ionization and H-minus free-free opacity, as in Mira variables, require a central star with >10^5 Lsun, which is greater than infrared observations allow. However, if significant local heating associated with accretion occurs, lower total luminosities are possible. Alternatively, photo-ionization from an early B-type star and p+/e- bremsstrahlung can explain our observations, and Orion-I may be an example of ionized accretion disk surrounding a forming massive star. Such accretion disks may not be able to form planets efficiently.

We report a three year, 18 epoch, VLBI monitoring study of H2O masers in the sub-parsec, warped, accretion disk within the NGC4258 AGN. Our immediate goals are to trace the geometry of the underlying disk, track rotation via measurement of proper motion, and ascertain the radii of masers for which centripetal acceleration may be measured separately. The monitoring includes ~ 4 times as many epochs, ~ 3 times denser sampling, and tighter control over sources of systematic error than earlier VLBI investigations. Coverage of a ~ 2400 km/s bandwidth has also enabled mapping of molecular material ~ 30% closer to the black hole than accomplished previously, which will strengthen geometric and dynamical disk models. Through repeated observation we have also measured for the first time a 5 microarsecond (1 sigma) thickness of the maser medium. Assuming this corresponds to the thickness of the accretion disk, hydrostatic equilibrium requires a disk plane temperature of ~ 600 K. Our long-term goal is a geometric distance to NGC4258 that is accurate to ~ 3%, a ~ 2 times improvement over the current best estimate. A geometric estimate of distance can be compared to distances obtained from analysis of Cepheid light curves, with the intent to recalibrate the extragalactic distance scale with reduced systematic uncertainties. This is the first paper in a series. We present here VLBI observations, data reduction, and temporal and spatial characteristics of the maser emission. Later papers will report estimation of orbital acceleration and proper motion, modeling of disk 3-D geometry and dynamics, and estimation of a "maser distance." Estimation of a "Cepheid distance" is presented in a parallel paper series.

We report the discovery of water maser emission in five active galactic nuclei (AGN) with the 100-m Green Bank Telescope (GBT). The positions of the newly discovered masers, measured with the VLA, are consistent with the optical positions of the host nuclei to within 1 sigma (0.3 arcsec radio and 1.3 arcsec optical) and most likely mark the locations of the embedded central engines. The spectra of three sources, 2MASX J08362280+3327383, NGC 6264, and UGC 09618 NED02, display the characteristic spectral signature of emission from an edge-on accretion disk with maximum orbital velocity of ~700, ~800, and ~1300 km s^-1, respectively. We also present a GBT spectrum of a previously known source MRK 0034 and interpret the narrow Doppler components reported here as indirect evidence that the emission originates in an edge-on accretion disk with orbital velocity of ~500 km s^-1. We obtained a detection rate of 12 percent (5 out of 41) among Seyfert 2 and LINER systems with 10000 km s^-1 < v_sys < 15000 km s^-1. For the 30 nuclear water masers with available hard X-ray data, we report a possible relationship between unabsorbed X-ray luminosity (2-10 keV) and total isotropic water maser luminosity, L_{2-10} proportional to L_{H2O}^{0.5+-0.1}, consistent with the model proposed by Neufeld and Maloney in which X-ray irradiation and heating of molecular accretion disk gas by the central engine excites the maser emission.

We present initial results from a time-series BVI survey of two fields in NGC 4258 using the Advanced Camera for Surveys onboard the Hubble Space Telescope. This galaxy was selected because of its accurate maser-based distance, which is anticipated to have a total uncertainty of ~3%. The goal of the HST observations is to provide an absolute calibration of the Cepheid Distance Scale and to measure its dependence on chemical abundance (the so-called "metallicity effect"). We carried out observations of two fields at different galactocentric distances with a mean abundance difference of 0.5 dex. We discovered a total of 281 Cepheids with periods ranging from 4 to 45 days (the duration of our observing window). We determine a Cepheid distance modulus for NGC 4258 (relative to the LMC) of 10.88 +- 0.04 (random) +- 0.05 (systematic) mag. Given the published maser distance to the galaxy, this implies \mu (LMC)=18.41 +- 0.10 (r) +- 0.13 (s) mag or D(LMC)= 48.1 +- 2.3 (r) +- 2.9 (s) kpc. We measure a metallicity effect of \gamma=-0.29 +- 0.09 (r) +- 0.05 (s) mag/dex. We see no evidence for a variation in the slope of the Period-Luminosity relation as a function of abundance. We estimate a Hubble Constant of H_0= 74 +- 3 (r) +- 6 (s) km/s Mpc using a recent sample of 4 well-observed type Ia SNe and our new calibration of the Cepheid Distance Scale. It may soon be possible to measure the value of H_0 with a total uncertainty of 5%, with consequent improvement in the determination of the equation of state of dark energy.

We compare the number of detected 22 GHz H2O masers in the Local Group galaxies M31, M33, NGC6822, IC10, IC1613, DDO187, GR8, NGC185, and the Magellanic Clouds with the water maser population of the Milky Way. To accomplish this we searched for water maser emission in the two Local Group galaxies M33 and NGC6822 using the Very Large Array (VLA) and incorporated results from previous studies. We observed 62 HII regions in M33 and 36 regions with H-alpha emission in NGC6822. Detection limits are 0.0015 and 0.0008 Lsun for M33 and NGC6822, respectively (corresponding to 47 and 50 mJy in three channels with 0.7 km/s width). M33 hosts three water masers above our detection limit, while in NGC6822 no maser source was detected. We find that the water maser detection rates in the Local Group galaxies M31, M33, NGC6822, IC1613, DDO187, GR8, NGC185, and the Magellanic Clouds are consistent with expectations from the Galactic water masers if one considers the different star formation rates of the galaxies. However, the galaxy IC10 exhibits an overabundance of masers, which may result from a compact central starburst.

Having conducted a search for the 22 GHz water vapor line towards galaxies with nuclear activity, large nuclear column densities or high infrared luminosities, we present H_2O spectra for NGC2273, UGC5101 and NGC3393 with isotropic luminosities of 7, 1500, and 400 L_sun. The H_2O maser in UGC5101 is by far the most luminous yet found in an ultraluminous infrared galaxy. NGC3393 reveals the classic spectrum of a `disk maser', represented by three distinct groups of Doppler components. As in all other known cases except NGC4258, the rotation velocity of the putative masing disk is well below 1000 km/s. Based on the literature and archive data, X-ray absorbing column densities are compiled for the 64 galaxies with reported maser sources beyond the Magellanic Clouds. For NGC2782 and NGC5728, we present Chandra archive data that indicate the presence of an active galactic nucleus in both galaxies. The correlation between absorbing column and H_2O emission is analyzed. There is a striking difference between kilo- and megamasers with megamasers being associated with higher column densities. All kilomasers (L_H_2O < 10 L_sun) except NGC2273 and NGC5194 are Compton-thin, i.e. their absorbing columns are < 10^24 cm^-2. Among the H_2O megamasers, 50% arise from Compton-thick and 85% from heavily obscured (> 10^23 cm^-2) active galactic nuclei. These values are not larger but consistent with those from samples of Seyfert 2 galaxies not selected on the basis of maser emission. The similarity in column densities can be explained by small deviations in position between maser spots and nuclear X-ray source and a high degree of clumpiness in the circumnuclear interstellar medium.

We report the first detection of an extragalactic millimeter wavelength H2O maser at 183 GHz towards NGC 3079 using the Submillimeter Array (SMA), and a tentative submillimeter wave detection of the 439 GHz maser towards the same source using the James Clerk Maxwell Telescope (JCMT). These H2O transitions are known to exhibit maser emission in star-forming regions and evolved stars. NGC 3079 is a well-studied nuclear H2O maser source at 22 GHz with a time-variable peak flux density in the range 3 -- 12 Jy. The 183 GHz H2O maser emission, with peak flux density $\sim$0.5 Jy (7$\sigma$ detection), also originates from the nuclear region of NGC 3079 and is spatially coincident with the dust continuum peak at 193 GHz (53 mJy integrated). Peak emission at both 183 and 439 GHz occurs in the same range of velocity as that covered by the 22 GHz spectrum. We estimate the gas to dust ratio of the nucleus of NGC 3079 to be $\approx$150, comparable to the Galactic value of 160. Discovery of maser emission in an active galactic nucleus beyond the long-known 22 GHz transition opens the possibility of future position-resolved radiative transfer modeling of accretion disks and outflows $<1$ pc from massive black holes.

We report the discovery of water maser emission in eight active galactic nuclei (AGN) with the 70-m NASA Deep Space Network (DSN) antennas at Tidbinbilla, Australia and Robledo, Spain. The positions of the newly discovered masers, measured with the VLA, are consistent with the optical positions of the host nuclei to within 1 sigma (0.3 arcsec radio and 1.3 arcsec optical) and most likely mark the locations of the embedded central engines. The spectra of two sources, NGC 3393 and NGC 5495, display the characteristic spectral signature of emission from an edge-on accretion disk, with orbital velocities of ~600 and ~400 km s^-1, respectively. In a survey with DSN facilities of 630 AGN selected from the NASA Extragalactic Database, we have discovered a total of 15 water maser sources. The resulting incidence rate of maser emission among nearby (v_sys < 7000 km s^-1) Seyfert 1.8-2.0 and LINER systems is ~10 percent for a typical rms noise level of ~14 mJy over 1.3 km s^-1 spectral channels. As a result of this work, the number of nearby AGN (v_sys < 7000 km s^-1) observed with <20 mJy rms noise has increased from 130 to 449.

Observation of the H2O megamaser galaxy IC 2560 with the Chandra Observatory reveals a complex spectrum composed of soft X-ray emission due to multi-temperature thermal plasma, and a hard continuum with strong emission lines. The continuum is most likely a Compton reflection (reprocessing) of primary emission that is completely absorbed at least up to 7 keV. The lines can be identified with fluorescence from Si, S and Fe in the lowest ionization stages. The equivalent widths of the Si and S lines are broadly compatible with those anticipated for reprocessing by optically thick cold plasma of Solar abundances, while the large equivalent width of the Fe line requires some overabundance of iron. A contribution to the line from a transmitted component cannot be ruled out, but the limits on the strength of the Compton shoulder make it less likely. From the bolometric luminosity of the nuclear region, we infer that the source radiates at 1 - 10% of its Eddington luminosity, for an adopted central mass of 3 million Solar masses. The overall spectrum is consistent with the hypotheses that the central engines powering the detected megamsers in accretion disks are obscured from direct view by the associated accretion disk material itself, and that there is a correlation between the occurrence of megamaser emission and Compton-thick absorption columns. For the 11 known galaxies with both column density measurements and maser emission believed to arise from accretion disks, eight AGN are Compton thick.

We present the first submm (865mum) imaging spectral line survey at one arcsecond resolution conducted with the Submillimeter Array toward Orion-KL. Within the 2x2 GHz bandpasses (lower and upper sidebands, 337.2-339.2GHz and 347.2-349.2GHz), we find about 145 spectral lines from 13 species, 6 isotopologues, and 5 vibrational excited states. Most nitrogen-bearing molecules are strong toward the hot core, whereas the oxygen-bearing molecules peak toward the south-west in the so-called compact ridge. Imaging of spectral lines is shown to be an additional tool to improve the identifications of molecular lines. Arcsecond spatial resolution allows us to distinguish the molecular line emission of the sources I and n from that of the hot core. The only molecular species detected strongly toward source I is SiO, delineating mainly the collimated north-east south-west low-velocity outflow. The two positions close to source I, which have previously been reported to show maser emission in the v=0 28SiO(1-0) and (2-1) lines, show no detectable maser emission in the v=0 28SiO(8-7) line at our spatial resolution. SiO is weak toward source n, and thus source n may not currently be driving a molecular outflow. CH$_3$OH is the molecule with the highest number of identified lines (46) in this spectral window. This ``line forest'' allows us to estimate temperatures in the region, and we find temperatures between 50 and 350K, with the peak temperatures occurring toward the hot core. The detection of strong vibrational excited line emission from the submm continuum peak SMA1 supports the interpretation that the source SMA1 is likely of protostellar nature.

A maximum likelihood analysis of the NGC 4258 maser positions and velocities reveals a ~2 sigma deviation from Keplerian motion in the projected rotation curve of the high-velocity features, corresponding to a ~9 km/s, or 0.8%, flattening of the LOS velocities with respect to Keplerian motion over the range of the high-velocity masers. While there are a number of potential explanations for this flattening, we argue for pure Keplerian rotation in an inclination-warped disk based on the ability of this model to explain a number of otherwise puzzling features of the system. A radial gradient in the disk inclination of 0.034 mas^{-1} is not only consistent with the observed rotation curve, but it generates a bowl along the near edge of the disk that naturally explains the otherwise puzzling narrow spread in the declinations of the systemic masers. It also explains the existence and location of an apparently recurring flare amongst the systemic masers. There is no significant evidence for non-Keplerian rotation in the inclination-warped disk. An additional implication of the inclination warp is that the disk rises in front of and obscures the central engine at a disk radius of about 8.3 mas, or 0.29 pc. By comparing the observed X-ray column to conditions in the disk at this radius, we argue that the molecular-to-atomic transition occurs just beyond the outermost maser at 0.28 pc and from this we infer an accretion rate of \~10^{-4}alpha M_sun/yr, consistent with the jet-dominated accretion models of Yuan et al. (2002). Our model suggests that most of the observed X-ray column arises in the warped accretion disk at 0.29 pc, and that the maser emission is truncated at large radii predominantly as a result of the molecular-to-atomic phase transition originally proposed by Neufeld & Maloney (1995).

We measured the angular rotation and proper motion of the Triangulum Galaxy (M33) with the Very Long Baseline Array by observing two H2O masers on opposite sides of the galaxy. By comparing the angular rotation rate with the inclination and rotation speed, we obtained a distance of 730 +/- 168 kiloparsecs. This distance is consistent with the most recent Cepheid distance measurement. M33 is moving with a velocity of 190 +/- 59 km/s relative to the Milky Way. These measurements promise a new method to determine dynamical models for the Local Group and the mass and dark matter halos of M31, M33 and the Milky Way.

We present an analysis of polarimetric observations at 22 GHz of the water vapor masers in NGC 4258 obtained with the VLA and the GBT. We do not detect any circular polarization in the spectrum indicative of Zeeman-induced splitting of the maser lines of water, a non-paramagnetic molecule. We have improved the 1-sigma upper limit estimate on the toroidal component of the magnetic field in the circumnuclear disk of NGC 4258 at a radius of 0.2 pc from 300 mG to 90 mG. We have developed a new method for the analysis of spectra with blended features and derive a 1-sigma upper limit of 30 mG on the radial component of the magnetic field at a radius of 0.14 pc. Assuming thermal and magnetic pressure balance, we estimate an upper limit on the mass accretion rate of ~10^(-3.7) M_sun/yr for a total magnetic field of less than 130 mG. We discuss the ramifications of our results on current maser models proposed to explain the observed maser emission structure and the consequences for current accretion theories. We find from our magnetic field limits that the thin-disk model and the jet-disk model are better candidates for accounting for the extremely low-luminosity nature of NGC 4258, than models that include advection-dominated accretion flows.

We report monitoring of the 0.3-10 keV spectrum of NGC4258 with XMM over 1.5 years.We als o report reprocessing of an overlapping series of archival Chandra observations. By including earlier ASCA and SAX observations, we present a new, nine-year time series of models fit to the X-ray spectrum of NGC4258. Over the nine years, the photoelectric absorbing column (~10^23 cm^-2) did not vary detectably, except for a ~40% drop between two ASCA epochs separated by 3 years and a ~60% rise between two XMM epochs separated by just 5 months. In contrast, factor of 2-3 changes are seen in absorbed flux on the timescale of years. These are uncorrelated with changes in absorbing column and indicative of central engine variability. The most rapid change in luminosity (5-10 keV) that we detect is ~30% over 19 days. The warped disk, a known source of H2O maser emission in NGC4258, is believed to cross the line of sight to the central engine. We propose that the variations in absorbing column arise from inhomogeneities sweeping across the line of sight in the rotating disk at the radius where the disk crosses the line of sight. We estimate that the inhomogeneities are ~10^15 cm in size at the crossing radius of 0.29 pc, slightly smaller than the expected scale height of the disk. This result thus provides strong evidence that the warped accretion disk is the absorber. This is the first direct confirmation that obscuration in type-2 AGN may, in some cases, arise in thin, warped accretion disks, rather than in geometrically thick tori. We do not detect Fe Kalpha line emission in any of our XMM spectra. We do not observe evidence of absorption lines in the XMM or reprocessed Chandra data.

We have mapped, for the first time, the full velocity extent of the water maser emission in NGC 3079. The largely north-south distribution of emission, aligned with a kpc-scale molecular disk, and the segregation of blue- and red-shifted emission on the sky are suggestive of a nearly edge-on molecular disk on pc-scales. Positions and line-of-sight velocities of blue- and red-shifted maser emission are consistent with a central mass of ~2 x 10^6 Msun enclosed within a radius of ~0.4 pc. The corresponding mean mass density of 10^6.8 Msun pc^-3 is suggestive of a central black hole, which is consistent with the detection of hard X-ray excess (20-100 keV) and an Fe Kalpha line from the nucleus. Because the rotation curve traced by the maser emission is flat, the mass of the pc-scale disk is significant with respect to the central mass. Since the velocity dispersion of the maser features does not decrease with radius and constitutes a large fraction of the orbital velocity, the disk is probably thick and flared. The rotation curve and the physical conditions necessary to support maser emission imply a Toomre Q-parameter that is << 1. Thus, the disk is most likely clumpy, and we argue that it is probably forming stars. Overall, the accretion disk in NGC 3079 stands in contrast to the compact, thin, warped, differentially rotating disk in the archetypal maser galaxy NGC 4258 (abridged).

Using the Green Bank Telescope, we have conducted a survey for 1.3 cm water maser emission toward the nuclei of nearby active galaxies, the most sensitive large survey for H2O masers to date. Among 145 galaxies observed, maser emission was newly detected in eleven sources and confirmed in one other. Our survey targeted nearby (v < 12,000 km/s, mainly type 2 AGNs north of declination -20 degrees, and includes a few additional sources as well. We find that more than a third of Seyfert 2 galaxies have strong maser emission, though the detection rate declines beyond v about 5000 km/s due to sensitivity limits. Two of the masers discovered during this survey are found in unexpected hosts: NGC 4151 (Seyfert 1.5) and NGC 2782 (starburst). We discuss the possible relations between the large X-ray column to NGC 4151 and a possible hidden AGN in NGC 2782 to the detected masers. Four of the masers discovered here, NGC 591, NGC 4388, NGC 5728 and NGC 6323, have high-velocity lines symmetrically spaced about the systemic velocity, a likely signature of molecular gas in a nuclear accretion disk. The maser source in NGC 6323, in particular, reveals the classic spectrum of a "disk maser" represented by three distinct groups of Doppler components. Future single-dish and VLBI observations of these four galaxies could provide a measurement of the distance to each galaxy, and of the Hubble constant, independent of standard candle calibrations.

We present observations of the 22 GHz water vapor megamasers in the Circinus galaxy made with the Tidbinbilla 70m telescope. These observations confirm the rapid variability seen earlier by Greenhill et al (1997). We show that this rapid variability can be explained by interstellar scintillation, based on what is now known of the interstellar scintillation seen in a significant number of flat spectrum AGN. The observed variability cannot be fully described by a simple model of either weak or diffractive scintillation.

One of the most important problems in the study of active galaxies is understanding the detailed geometry, physics, and evolution of the central engines and their environments. The leading models involve an accretion disk and torus structure around a central dense object, thought to be a supermassive black hole. Gas found in the environment of AGN is associated with different structures: molecular accretion disks, larger scale atomic tori, ionized and neutral "cocoons" in which the nuclear regions can be embedded. All of them can be studied at radio wavelengths by various means. Here, we summarize the work that has been done to date in the radio band to characterize these structures. Much has been learned about the central few parsecs of AGN in the last few decades with contemporary instruments but the picture remains incomplete. In order to be able to define a more accurate model of this region, significant advances in sensitivity, spectral and angular resolution, and bandpass stability are required. The necessary advances will only be provided by the Square Kilometer Array and we discuss the possibilities that these dramatic improvements will open for the study of the gas in the central region of AGN.

The present-day Universe is seemingly dominated by dark energy and dark matter, but mapping the normal (baryonic) content remains vital for both astrophysics - understanding how galaxies form - and astro-particle physics - inferring properties of the dark components. The Square Kilometre Array (SKA) will provide the only means of studying the cosmic evolution of neutral Hydrogen (HI) which, alongside information on star formation from the radio continuum, is needed to understand how stars formed from gas within dark-matter over-densities and the roles of gas accretion and galaxy merging. `All hemisphere' HI redshift surveys to redshift 1.5 are feasible with wide-field-of-view realizations of the SKA and, by measuring the galaxy power spectrum in exquisite detail, will allow the first precise studies of the equation-of-state of dark energy. The SKA will be capable of other uniquely powerful cosmological studies including the measurement of the dark-matter power spectrum using weak gravitational lensing, and the precise measurement of H0 using extragalactic water masers. The SKA is likely to become the premier dark-energy-measuring machine, bringing breakthroughs in cosmology beyond those likely to be made possible by combining CMB (e.g. Planck), optical (e.g. LSST, SNAP) and other early-21st-century datasets.

High precision estimation of the equation of state of dark energy depends on constraints external to analyses of Cosmic Microwave Background fluctuations. A geometric estimation of the local expansion rate, H_o, would provide the most direct and robust constraint. Traditional techniques to estimate H_o have depended on observations of standard candles for which systematic effects can be 10% or more. Observations of water maser sources in the accretion disks that feed the central engines of active galaxies enable simplified, robust, and largely geometric analyses. Many thousand maser sources will be discovered in studies with the SKA, owing to its great sensitivity. Spectroscopic monitoring and interferometric mapping - with intercontinental baselines - will allow estimation of H_o to 1% and possibly better.

We report on the results of a 14 year long VLBI study of proper motions in the IC133 H2O maser source in the galaxy M33. The method of Ordered Motion Parallax was used to model the 3-dimensional structure and dynamics of IC133 and obtain a distance estimate, 800 +/- 180 kpc. Our technique for determining the distance to M33 is independent of calibrations common to other distance indicators, such as Cepheid Period-Luminosity relations, and therefore provides an important check for previous distance determinations.

The 25 years following the serendipitous discovery of megamasers have seen tremendous progress in the study of luminous extragalactic H$_2$O emission. Single-dish monitoring and high resolution interferometry have been used to identify sites of massive star formation, to study the interaction of nuclear jets with dense molecular gas and to investigate the circumnuclear environment of active galactic nuclei (AGN). Accretion disks with radii of 0.1--3 pc were mapped and masses of nuclear engines of order 10$^{6}$--10$^{8}$ M$_{\odot}$ were determined. So far, $\sim$50 extragalactic H$_2$O maser sources have been detected, but few have been studied in detail.

We present the first 865 mu continuum image with sub-arcsecond resolution obtained with the Submillimeter Array. These data resolve the Orion-KL region into the hot core, the nearby radio source I, the sub-mm counterpart to the infrared source n (radio source L), and new sub-mm continuum sources. The radio to submillimeter emission from source I may be modeled as either the result of proton-electron free-free emission that is optically thick to ~100 GHz plus dust emission that accounts for the majority of the submillimeter flux, or H- free-free emission that gives rise to a power-law spectrum with power-law index of ~1.6. The latter model would indicate similar physical conditions as found in the inner circumstellar environment of Mira variable stars. Future sub-arcsecond observations at shorter sub-mm wavelengths should easily discriminate between these two possibilities. The sub-mm continuum emission toward source n can be interpreted in the framework of emission from an accretion disk.

We report the discovery of a radio-loud flat-spectrum QSO at z=5.47 with properties similar to those of the EGRET gamma-ray blazars. This source is the brightest radio QSO at z>5, with a pc-scale radio jet and a black hole mass estimate >10^{10}M_\odot. It appears to be the most distant blazar discovered to date. High energy observations of this source can provide powerful probes of the background radiation in the early universe.

We present Keck LWS images of the Orion BN/KL star forming region obtained in the first multi-wavelength study to have 0.3-0.5" resolution from 4.7 to 22 microns. The young stellar objects designated infrared source-n and radio source-I are believed to dominate the BN/KL region. We have detected extended emission from a probable accretion disk around source-n but infer a stellar luminosity on the order of only 2000 Lsun. Although source-I is believed to be more luminous, we do not detect an infrared counterpart even at the longest wavelengths. However, we resolve the closeby infrared source, IRc2, into an arc of knots ~1000 AU long at all wavelengths. Although the physical relation of source-I to IRc2 remains ambiguous, we suggest these sources mark a high density core (10^7-10^8 pc^-3 over 1000 AU) within the larger BN/KL star forming cluster. The high density may be a consequence of the core being young and heavily embedded. We suggest the energetics of the BN/KL region may be dominated by this cluster core rather than one or two individual sources.

We report the discovery of two binary millisecond pulsars in the core-collapsed globular cluster M30 using the Green Bank Telescope (GBT) at 20 cm. PSR J2140-2310A (M30A) is an eclipsing 11-ms pulsar in a 4-hr circular orbit and PSR J2140-23B (M30B) is a 13-ms pulsar in an as yet undetermined but most likely highly eccentric (e>0.5) and relativistic orbit. Timing observations of M30A with a 20-month baseline have provided precise determinations of the pulsar's position (within 4" of the optical centroid of the cluster), and spin and orbital parameters, which constrain the mass of the companion star to be m_2 >~ 0.1Msun. The position of M30A is coincident with a possible thermal X-ray point source found in archival Chandra data which is most likely due to emission from hot polar caps on the neutron star. In addition, there is a faint (V_555 ~ 23.8) star visible in archival HST F555W data that may be the companion to the pulsar. Eclipses of the pulsed radio emission from M30A by the ionized wind from the compact companion star show a frequency dependent duration (\propto\nu^{-\alpha} with \alpha ~ 0.4-0.5) and delay the pulse arrival times near eclipse ingress and egress by up to 2-3 ms. Future observations of M30 may allow both the measurement of post-Keplerian orbital parameters from M30B and the detection of new pulsars due to the effects of strong diffractive scintillation.

High-mass star formation is not well understood chiefly because examples are deeply embedded, relatively distant, and crowded with sources of emission. Using VLA and VLBA observations of water and SiO maser emission, we have mapped in detail the structure and proper motion of material 20-500 AU from the closest high-mass YSO, radio source-I in the Orion KL region. We observe streams of material driven in a rotating, wide angle, bipolar wind from the surface of an edge-on accretion disk. The example of source-I provides strong evidence that high-mass star formation proceeds via accretion

We report on our ongoing, high-resolution study of H2O masers in the innermost parsec of NGC 4258. Over thirty epochs of VLBA and VLA data, taken over six years, are being used to monitor the velocities, accelerations, positions and proper motions of water masers rotating in a warped, Keplerian disk about a supermassive central object. Our extensive monitoring results in an improved accuracy in the distance determination to this galaxy (i) via a reduction in experimental random errors due to a longer time-baseline than in previous work; and (ii) via a better modeling of sources of systematic error, such as disk eccentricity. These data can therefore yield an extremely accurate geometric distance to NGC 4258.

We present the first VLBI maps of H2O maser emission (lambda 1.3cm) in the nucleus of the Circinus Galaxy, constructed from data obtained with the Australia Telescope Long Baseline Array. The maser emission traces a warped, edge-on accretion disk between radii of 0.11+/-0.02 and ~0.40 pc, as well as a wide-angle outflow that extends up to ~1 pc from the estimated disk center. The disk rotation is close to Keplerian (v varies as 1/sqrt(r)), the maximum detected rotation speed is 260 km/s, and the inferred central mass is 1.7+/-0.3 x 10^6 solar masses. The outflowing masers are irregularly distributed above and below the disk, with relative outflow velocities up to ~+/-160 km/s, projected along the line of sight. The flow probably originates closer than 0.1 pc to the central engine, possibly in an inward extension of the accretion disk, though there is only weak evidence of rotation in the outward moving material. We observe that the warp of the disk appears to collimate the outflow and to fix the extent of the ionization cone observed on larger angular scales. This study provides the first direct evidence (i.e., through imaging) of dusty, high-density, molecular material in a nuclear outflow <1 pc from the central engine of a Seyfert galaxy, as well as the first graphic evidence that warped accretion disks can channel outflows and illumination patterns in AGN. We speculate that the same arrangement, which in some ways obviates the need for a geometrically thick, dusty torus, may apply to other type-2 AGN.

We analyze recent simultaneous Chandra/RXTE observations of the Seyfert 2 galaxy NGC 4945. The unprecedented spatial resolution of Chandra means we are able to separate the spectra of the nucleus, starburst and superwind regions, while the RXTE data extend the spectrum to higher energies. The extreme absorbing column of $N_H\sim 4\times 10^{24}$ cm$^{-2}$ means that the nucleus is only seen directly above 8--10 keV, while the lower energy spectrum from the nuclear region in Chandra is dominated by reflection. By contrast, the superwind is dominated by emission from hot plasma, but the starburst region contains both hot plasma and reflection signatures. To form a reflected spectrum requires that the starburst region contains clumps of cool, optically thick material, perhaps star forming cores, which are irradiated by 7--10 keV photons from the nucleus. Since photons of this energy are obscured along the line sight then this confirms the result of Madejski et al. (2000) that the extreme absorbtion material is disk-like rather than a torus. However, the IR/optical limits on the lack of high excitation emission lines show that by contrast the lower energy photons from the nucleus are obscured in all directions. We discuss the complex absorption structure revealed by these observations, and propose an an overall source geometry in which the nucleus is completely embedded in material with $N_H\sim 10^{23}$ cm$^{-2}$

We report the discovery of H2O maser emission at 1.35 cm wavelength in seven active galactic nuclei (at distances up to <80 Mpc) during a survey conducted at the 70-m diameter antenna of the NASA Deep Space Network near Canberra, Australia. The detection rate was (approx.) 4%. Two of the maser sources are particularly interesting because they display satellite high-velocity emission lines, which are a signature of emission from the accretion disks of supermassive black holes when seen edge on. Three of the masers are coincident, to within uncertainties of 0.''2, with continuum emission sources we observed at about (lamda)1.3 cm. We also report the discovery of new spectral features in the Circinus galaxy H2O maser that broaden the known velocity range of emission therein by a factor of (approx.) 1.7. If the new spectral features originate in the Circinus accretion disk, then molecular material must survive at radii (approx.) 3 times smaller than had been believed previously [(approx.) 0.03 pc or (approx.) 2 x 105 Schwarzschild radii].

Water vapor emission at 22 GHz is reported from the nucleus of the LINER galaxy Mrk 1419 (NGC 2960). Single-dish spectra of the maser source show properties that are similar to those seen in NGC 4258, namely (1) a cluster of systemic H2O features, (2) two additional H2O clusters, one red- and one blue-shifted by about 475 km/s, (3) a likely acceleration of the systemic features, and (4) no detectable velocity drifts in the red- and blue-shifted features. Interpreting the data in terms of the paradigm established for NGC 4258, i.e. assuming the presence of an edge-on Keplerian circumnuclear annulus with the systemic emission arising from the near side of its inner edge, the following parameters are derived: Rotational velocity: 330-600 km/s; radius: 0.13-0.43 pc; binding mass: about 10 million solar masses. With the galaxy being approximately ten times farther away than NGC 4258, a comparison of linear and angular scales (the latter via Very Long Baseline Interferometry) may provide an accurate geometric distance to Mrk 1419 that could be used to calibrate the cosmic distance scale.

We report the discovery of seven new millisecond pulsars during 20 cm searches of 19 globular clusters using the recently upgraded Arecibo Telescope and the new Green Bank Telescope (GBT). Five of the pulsars are in compact binaries, three of which are eclipsing. One of the systems is in an eccentric orbit which is almost certainly highly relativistic. Additional searches and timing observations are underway.

We have mapped the water maser emission associated with the infrared centers IRS1 and IRS3 of the NGC 2071IR star forming region at four epochs over ~4 months with the Very Long Baseline Array (VLBA). We detected 269 maser features with ~1 km/sec linewidths and measured 30 proper motions. In each infrared center, the water maser emission appears to trace parts of a protostellar disk and collimated outflow. The disk components are ~9 and ~17 AU long, in IRS3 and IRS1 respectively, and ~2 AU wide. They are identified as disks by their compact size, elongation parallel to the direction of known IR polarization, central location in the maser maps, small internal proper motions, and proximity to 1.3 cm continuum emission. The outflows have axes perpendicular to the disks and exhibit proper motions of up to ~42 km/sec. They are outlined by maser emission up to ~260 AU from the protostars. The IRS3 outflow appears to be conical on one side, while the IRS1 outflow comprises a narrowly collimated bipolar flow surrounded by outward-facing, funnel-shaped cavities. The detection of water maser emission tracing such compact disk components and specifically conical or funnel-shaped structures is unusual. The fact that the distributions are similar in IRS3 and IRS1 may indicate the two infrared centers are roughly coeval. NGC 2071IR provides a rare opportunity to resolve the structures and dynamics of disks and outflows together, and to do so for two protostars that are only ~2000 AU apart (in projection) in a deeply embedded star forming region of intermediate luminosity.

Key and still largely missing parameters for measuring the mass content and distribution of the Local Group are the proper motion vectors of its member galaxies. The problem when trying to derive the gravitational potential of the Local Group is that usually only radial velocities are known, and hence statistical approaches have to be used. The expected proper motions for galaxies within the Local Group, ranging from 20 to 100 $\mu$as/yr, are detectable with VLBI using the phase-referencing technique. We present phase-referencing observations of bright masers in IC~10 and M33 with respect to background quasars. We observed the H$_2$O masers in IC10 three times over a period of two months to check the accuracy of the relative positions. The relative positions were obtained by modeling the interferometer phase data for the maser sources referenced to the background quasars. The model allowed for a relative position shift for the source and a single vertical atmospheric delay error in the correlator model for each antenna. The rms of the relative positions for the three observations is only 0.01 mas, which is approximately the expected position error due to thermal noise. Also, we present a method to measure the geometric distance to M33. This will allow re-calibration of the extragalactic distance scale based on Cepheids. The method is to measure the relative proper motions of two H$_2$O maser sources on opposite sides of M33. The measured angular rotation rate, coupled with other measurements of the inclination and rotation speed of the galaxy, yields a direct distance measurement.

We present the results of the first high angular resolution observations of SiO maser emission towards the star forming region W51-IRS2 made with the Very Large Array (VLA) and Very Long Baseline Array (VLBA). Our images of the water maser emission in W51-IRS2 reveal two maser complexes bracketing the SiO maser source. One of these water maser complexes appears to trace a bow shock whose opening angle is consistent with the opening angle observed in the distribution of SiO maser emission. A comparison of our water maser image with an image constructed from data acquired 19 years earlier clearly shows the persistence and motion of this bow shock. The proper motions correspond to an outflow velocity of 80 km/s, which is consistent with the data of 19 years ago (that spanned 2 years). We have discovered a two-armed linear structure in the SiO maser emission on scales of ~25 AU, and we find a velocity gradient on the order of 0.1 km/s/AU along the arms. We propose that the SiO maser source traces the limbs of an accelerating bipolar outflow close to an obscured protostar. We estimate that the outflow makes an angle of <20 degrees with respect to the plane of the sky. Our measurement of the acceleration is consistent with a reported drift in the line-of-sight velocity of the W51 SiO maser source.

We report the cumulative results of five surveys for water maser emission at 1.35 cm wavelength in 131 active galactic nuclei (AGNs) and star-forming galaxies, conducted at the Parkes Observatory between 1993 and 1998. We detected one new maser, in the edge-on galaxy IRASF01063-8034, which exhibits a single, about 0.1 Jy spectral feature at 4282+/-6 km/s (heliocentric) with an unusually large 54+/-16 km/s half-power full width. The centroid velocity of the emission increased to 4319.6+/-0.6 km/s (38+/-2 km/s width) over the 13 days between discovery and confirmation of the detection. A similarly broad linewidth and large change in velocity has been noted for the maser in NGC1052, wherein jet activity excites the emission. Neither optical spectroscopy, radio-infrared correlations, nor infrared colors provide compelling evidence of unusual activity in the nucleus of IRASF01063-8034. Since the galaxy appears to be outwardly normal at optical and infrared wavelengths, detection of a water maser therein is unique. The maser emission is evidence that the galaxy harbors an AGN that is probably obscured by the edge-on galactic disk. The detection highlights the possibility that undetected AGNs could be hidden in other relatively nearby galaxies. No other maser emission features have been identified at velocities between 3084 km/s and 6181 km/s.

A small number of dusty Wolf-Rayet stars have been resolved into pinwheel nebulae, defined by their ``rotating'' spiral dust shells observed in the infrared. This morphology is naturally explained by dust formation associated with colliding winds in a binary system. In order to confirm and further explore this hypothesis, we have observed the known pinwheel nebulae (WR 104 and WR 98a) as well as the suspected binary WR 112 at multiple radio wavelengths with the Very Large Array to search for non-thermal radio emission from colliding winds. The spectrum of each target is nearly flat between 5 and 22 GHz, consistent with the presence of non-thermal emission that is reduced at low frequencies by free-free absorption. This emission must lie outside the radio ``photosphere,'' leading us to estimate a lower limit to the physical size of the non-thermal emitting region that is larger than expected from current theory. Based on a radio and infrared comparison to WR 104 and 98a, we conclude that WR 112 is a likely candidate pinwheel nebula, but its temporal variability indicates an eccentric binary orbit or a pinwheel viewed nearly edge-on. A sensitive radio survey of IR-bright WRs would stringently test the hypothesis that colliding winds lie at the heart of_all_ dusty WR systems. We also discuss the effects of dust obscuration in the ultra-violet and how radio-determined mass-loss rates of pinwheel nebulae (and dusty WR stars in general) may be underestimated due to shadowing effects.

Study of extragalactic H2O masers has progressed significantly in the 25 years since their discovery. Existing in star forming regions and in the accretion disks supermassive black holes, they are familiar and unfamiliar at the same time. A review of how our understanding has grown, up to the present day, is followed by comments on future prospects.

This paper discusses our ongoing efforts to characterize dust-enshrouded Wolf-Rayet (WR) stars in the radio and infrared. We have used the Very Large Array to measure the broadband radio spectrum of WR stars in suspected binary systems and discovered non-thermal emission, which is usually attributed to colliding winds. In addition, infrared imaging using aperture masking interferometry on the Keck-I telescope has resolved the dust shells around a number of WR stars with K-magnitudes brighter than ~6. Although this admittedly small study suffers from selection bias, we note that all the dust-enshrouded WR stars with radio detections show evidence for colliding winds, supporting the theory that wind compression in a binary system is necessary for efficient dust production. A consequence of this hypothesis is that virtually all WC8-10 stars must be in binaries, since most are dusty. Single-star and binary stellar evolution models will have to be modified to accommodate this observational result if confirmed.

The first VLBI images of water maser emission in the Circinus Galaxy AGN show both a warped, edge-on accretion disk and an outflow 0.1 to 1 pc from the central engine. The inferred central mass is 1.3 million suns, while the disk mass may be on the order of 0.1 million suns, based on a nearly Keplerian rotation curve. The bipolar, wide-angle outflow appears to contain ``bullets'' ejected from within <0.1 pc of the central mass. The positions of filaments and bullets observed in the AGN ionization cone on kpc-scales suggest that the disk channels the flow to a radius of about 0.4 pc, at which the flow appears to disrupt the disk.

The first VLBI images of water maser emission in the Circinus Galaxy AGN show both an accretion disk and outflow 0.1 to 1 pc from the central engine. The maser traces a warped, edge-on accretion disk between radii of about 0.1 and 0.4 pc that is bound by a 1.3x10^6 Msun central mass. The rotation curve is somewhat flatter than a Keplerian rotation law and is consistent with a disk mass on the order of 10^5 Msun. Clumpy substructure may display a kinematic signature consistent with spiral arms. Away from the disk, a second population of water masers traces broadly a bipolar, wide-angle outflow that contains (bullet-like) clumps ejected from a region centered on the central engine and <0.1 pc in radius. Out to a radius of about 0.4 pc, the warp of the accretion disk appears to channel the outflow, the orientation of which on the sky coincides with the orientation of the known kiloparsec-scale ionization cone. Beyond this radius, the flow crosses the disk and truncates it. This suggests that the current accretion event and associated disk-outflow geometry have a lifetime on the order of 10^7 years.

We report the detection of a new 3.06 ms binary pulsar in the globular cluster NGC6544 using a Fourier-domain ``acceleration'' search. With an implied companion mass of ~0.01 solar masses and an orbital period of only P_b~1.7 hours, it displays very similar orbital properties to many pulsars which are eclipsed by their companion winds. The orbital period is the second shortest of known binary pulsars after 47 Tuc R. The measured flux density of 1.3 +/- 0.4 mJy at 1332 MHz indicates that the pulsar is almost certainly the known steep-spectrum point source near the core of NGC6544.

We present VLBA observations of water maser emission associated with the star-forming region IRAS 21391+5802, which is embedded in a bright rimmed cometary globule in IC1396. The angular resolution of the maps is about 0.8 mas, corresponding to a spatial resolution of about 0.6 AU, at an estimated distance of 750 pc. Proper motions are derived for 10 maser features identified consistently over three epochs, which were separated by intervals of about one month. The masers appear in four groups, which are aligned linearly on the sky, roughly along a northeast-southwest direction, with a total separation of about 520 AU (about 0.7 arcseconds). The 3-D velocities of the masers have a maximum value of about 42 km/s (about 9 AU/yr). The average error on the derived proper motions is about 4 km/s. The overall pattern of proper motions is indicative of a bipolar outflow. Proper motions of the masers in a central cluster, with a projected extent of about 20 AU, show systematic deviations from a radial outflow. However, we find no evidence of Keplerian rotation, as has been claimed elsewhere. A nearly circular loop of masers lies near the middle of the cluster. The radius of this loop is 1 AU and the line-of-sight velocities of the masers in the loop are within 2 km/s of the systemic velocity of the region. These masers presumably exist at the radial distance where significant dust condensation occurs in the outflow emanating from the star.

The water masers in NGC4258 delineate the structure and dynamics of a sub-parsec-diameter accretion disk around a supermassive black hole. We report acceleration measurements for 12 systemic-velocity spectral features and 19 high-velocity spectral features, using a total of 25 epochs of observations from Effelsberg (5 epochs), the VLA (15 epochs), and the VLBA (5 epochs) spanning the years 1994 to 1997. The measured accelerations of the systemic-velocity features are between 7.5 and 10.4 km/s/yr, and there is no evidence for a dip in the spectrum at the systemic velocity. Such a dip has been attributed in the past to an absorbing layer of non-inverted H2O (Watson & Wallin 1994; Maoz & McKee 1998). The accelerations of the high-velocity features, measured here for the first time, range from -0.77 to 0.38 km/s/yr. From the line-of-sight accelerations and velocities, we infer the positions of these high-velocity masers with a simple edge-on disk model. The resulting positions fall between -13.6 and 9.3 degrees in azimuth (measured from the midline, the diameter perpendicular to the line of sight). A model that suggests a spiral shock origin of the masers (Maoz & McKee 1998), in which changes in maser velocity are due to the outward motion of the shock wave, predicts apparent accelerations of -0.05 times theta_p/2.5 degrees km/s/yr, where theta_p is the pitch angle of the spiral arms. Our data are not consistent with these predictions. We also discuss the physical properties of the high-velocity masers. Most notably, the strongest high-velocity masers lie near the midline, where the velocity gradient is smallest, thereby providing the longest amplification path lengths. (Abridged)

Naturally occurring water vapor maser emission at 1.35 cm wavelength provides an accurate probe for the study of accretion disks around highly compact objects, thought to be black holes, in the centers of active galaxies. Because of the exceptionally fine angular resolution, 200 microarcseconds, obtainable with very long baseline interferometry, accompanied by high spectral resolution, < 0.1 km/s, the dynamics and structures of these disks can be probed with exceptional clarity. The data on the galaxy NGC4258 are discussed here in detail. The mass of the black hole binding the accretion disk is 3.9 times 10^7 solar masses. Although the accretion disk has a rotational period of about 800 years, the physical motions of the masers have been directly measured with VLBI over a period of a few years. These measurements also allow the distance from the earth to the black hole to be estimated to an accuracy of 4 percent. The status of the search for other maser/black hole candidates is also discussed.

The water maser in the mildly active nucleus in the nearby galaxy NGC4258 traces a thin, nearly edge-on, subparsec-scale Keplerian disk. Using the technique of very long baseline interferometry, we have detected the proper motions of these masers as they sweep in front of the central black hole at an orbital velocity of about 1100 km/s. The average maser proper motion of 31.5 microarcseconds per year is used in conjunction with the observed acceleration of the masers to derive a purely geometric distance to the galaxy of 7.2 +- 0.3 Mpc. This is the most precise extragalactic distance measured to date, and, being independent of all other distance indicators, is likely to play an important role in calibrating the extragalactic distance scale.

The BN/KL region in the Orion molecular cloud is an archetype in the study of the formation of stars much more massive than the Sun. This region contains luminous young stars and protostars, but it is difficult to study because of overlying dust and gas. Our basic expectations are shaped to some extent by the present theoretical picture of star formation, the cornerstone of which is that protostars acrete gas from rotating equatorial disks, and shed angular momentum by ejecting gas in bipolar outflows. The main source of the outflow in the BN/KL region may be an object known as radio source I, which is commonly believed to be surrounded by a rotating disk of molecular material. Here we report high-resolution observations of silicon monoxide (SiO) and water maser emission from the gas surrounding source I; we show that within 60 AU (about the size of the Solar System), the region is dominated by a conical bipolar outflow, rather than the expected disk. A slower outflow, close to the equatorial plane of the protostellar system, extends to radii of 1,000 AU.

We report VLA 1 sigma upper limits of 1.5% and 3% on the intrinsic circular and linear fractional polarizations, respectively, of the water vapor maser emission 0.2 pc from the central engine of NGC 4258. A corresponding 0.5% upper limit on any Zeeman-splitting-induced circular polarization translates to a 1 sigma upper limit on the parallel, or toroidal, component of the magnetic field of 300 mG. Assuming magnetic and thermal pressure balance in the disk, this magnetic field upper limit corresponds to a model-dependent estimate of the accretion rate through the molecular disk of 10^-1.9 alpha solar masses per year for the case where the magnetic field lies along the line of sight.

We report a 3 sigma upper limit of 220 microJanskys on any 22-GHz continuum emission coincident with the central engine in NGC 4258. If NGC 4258 is powered by an advection-dominated accretion flow, this radio upper limit implies that the inner advection-dominated flow cannot extend significantly beyond about 100 Schwarzschild radii.

We have conducted VLBI observations of water maser and radio continuum emission in the nucleus of NGC 3079. The 22 GHz maser emission arises in compact clumps, distributed along an axis that is aligned with the major axis of the galactic disk. The velocities of the masers are consistent with their lying in the inner parsec of a molecular disk rotating in the same sense as the rest of the galaxy. However, the velocity field has a significant non-rotational component, which may indicate supersonic turbulence. The bright maser emission is not coincident with any detected compact 22 GHz continuum source, suggesting the high apparent luminosity of the maser may not due to beamed amplification of continuum emission. We observed two compact continuum sources that have inverted spectra between 5 and 8 GHz, and steep spectra between 8 and 22 GHz. NGC 3079 may be a nearby, low-luminosity example of the class of compact symmetric gigahertz-peaked spectrum (GPS) radio sources. We detected a third continuum component that lies along the same axis as the other two, strongly suggesting that this galaxy possesses a nuclear jet. Faint maser emission was detected near this axis, which may indicate a second population of masers associated with the jet.

We present a parsec-scale image of the OH maser in the nucleus of the active galaxy IIIZw35, made using the Very Long Baseline Array at a wavelength of 18 cm. We detected two distinct components, with a projected separation of 50 pc (for D=110 Mpc) and a separation in Doppler velocity of 70 km/s, which contain 50% of the total maser flux. Velocity gradients within these components could indicate rotation of clouds with binding mass densities of ~7000 solar masses per cubic parsec, or total masses of more than 500,000 solar masses. Emission in the 1665-MHz OH line is roughly coincident in position with that in the 1667-MHz line, although the lines peak at different Doppler velocities. We detected no 18 cm continuum emission; our upper limit implies a peak apparent optical depth greater than 3.4, assuming the maser is an unsaturated amplifier of continuum radiation.

We present the first interferometer map of the water maser emission in the active nucleus of NGC 4945. Although the declination of the galaxy is about -49 degrees, we were able to make the observations with the southernmost antennas of the Very Long Baseline Array. Strong maser emission is present in three velocity ranges, one near the systemic velocity and two shifted roughly symmetrically by +/-(100-150) km/s. This is the first detection of highly blue-shifted water emission in NGC 4945. We determined the position of the maser to be RA(B1950)= 13 02 32.28 +/- 0.02 ; Dec(B1950)= -49 12 01.9 +/- 0.1. The uncertainties in earlier estimates are at least several arcseconds. The maser lies within 2'' (36 pc at a distance of 3.7 Mpc) of the peaks in 1.4 GHz continuum and 1.6 micron emission from the nucleus. The mappable maser emission is distributed roughly linearly over about 40 milliarcseconds (0.7 pc) at a position angle of about 45 degrees, which is close to the 43 +/- 2 degree position angle of the galactic disk. The red and blue-shifted emission symmetrically stradle the systemic emission on the sky, which suggests material in edge-on circular motion around a central object. The position-velocity structure indicates a binding mass of about one million Suns, within a volume of radius about 0.3 pc. This implies that the central engine radiates on the order of 10% of its Eddington luminosity.

The water maser lines in the nucleus of the Circinus galaxy vary on time scales as short as a few minutes, at least two orders of magnitude less than for other Galactic or extragalactic water masers. The amplitude of one line more than doubled in about 10 minutes. The intensity changes cannot be attributed easily to a mechanism of intrinsic fluctuations. The variability may be the result of strong interstellar diffractive scintillation along the line of sight within our Galaxy, which would be the first example of diffractive scintillation for any source at 22 GHz and for any source other than a pulsar. However, only the very shortest timescales for interstellar scintillation, obtained from pulsar observations and scaled to 22 GHz, correspond to the observed maser variability. Alternatively, the intensity changes may be a reaction to fluctuations in compact background or radiative pump sources and thereby may be related to variability of the central engine. The maser spectral features symmetrically bracket the systemic velocity of the galaxy, offset by about +/-(100-200) km/s. The spectrum of the Circinus maser is similar in some respects to that of the maser in NGC 4258, which probably traces a molecular disk rotating around a supermassive object. VLBI observations could reveal whether the maser source in the heart of the Circinus galaxy is part of a similar dynamical system.

We have made the first VLBI synthesis images of the H2O maser emission associated with the central engine of the Seyfert galaxy NGC 1068. Emission extends about +/-300 km/s from the systemic velocity. Images with submilliarcsecond angular resolution show that the red-shifted emission lies along an arc to the northwest of the systemic emission. (The blue-shifted emission has not yet been imaged with VLBI.) Based on the maser velocities and the relative orientation of the known radio jet, we propose that the maser emission arises on the surface of a nearly edge-on torus, where physical conditions are conducive to maser action. The visible part of the torus is axially thick, with comparable height and radius. The velocity field indicates sub-Keplerian differential rotation around a central mass of about 1e7 Msun that lies within a cylindrical radius of about 0.65 pc. The estimated luminosity of the central engine is about 0.5 of the Eddington limit. There is no detectable compact radio continuum emission near the proposed center of the torus (T_B< 5e6 K on size scales of about 0.1 pc), so that the observed flat-spectrum core cannot be direct self-absorbed synchrotron radiation.