Digital Signal Processing using Stream High Performance Computing: A 512-input Broadband Correlator for Radio Astronomy

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 correlators comparatively fast development time and low cost establish a practical foundation for the scalability of a modular, heterogeneous, computing architecture.

Foreground Model and Antenna Calibration Errors in the Measurement of the Sky-Averaged lambda 21 cm Signal at z~20

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

A Scalable Hybrid FPGA/GPU FX Correlator

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.

A 189 MHz, 2400 square degree polarization survey with the Murchison Widefield Array 32-element prototype

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(alpha,delta) = (4^h,-27^circ.6)$ in terms of three dimensional power spectra

Masers in AGN environments

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