Person:

Greenhill, Lincoln

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 5 GiB s(^{−1}), grouped into 8 s 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 8 s data must be completely reduced before the next batch arrives. Maintaining real time operation will require a sustained performance of around 2.5 TFLOP s(^{−1}) (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.

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 the 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^{−1} yr^{−1}), indicating an origin within 13° of the disk midline (the perpendicular to the line of sight [LOS] to the black hole). Comparison of the projected positions of high-velocity emission 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^{−1} yr^{−1}), consistent with emission originating from a concavity where the thin, warped disk is tangent to the LOS. 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.

We present a Very Long Baseline Interferometry image of the water maser emission in the nuclear region of NGC 3393. The maser emission has a linear distribution oriented at a position angle of ~–34°, perpendicular to both the kiloparsec-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 ± 0.2) × 10^7 M☉) within 0.36 ± 0.02 pc (1.48 ± 0.06 mas), which corresponds to a mean mass density of (\sim10^{8.2} M☉ 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 ± 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 ± 0.02 pc, which is smaller than the inner disk radius (0.36 ± 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 NGC 4258 and NGC 1068. The outer radius of the disk as traced by the masers along the midline is about 1.5 pc.

Based on spectroscopic signatures, about one-third of known H2O maser sources in active galactic nuclei (AGNs) 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 (NGC 1320, NGC 17) in a survey of 40 inclined active galaxies (vsys < 20,000 km s−1). The remaining three disk maser candidates were identified in monitoring of known sources: NGC 449, NGC 2979, and NGC 3735. We also confirm a previously marginal case in UGC 4203. For the disk maser candidates reported here, inferred rotation speeds are 130–500 km s−1. Monitoring of three more rapidly rotating candidate disks (CG 211, NGC 6264, VV 340A) has enabled measurement of likely orbital centripetal acceleration, and estimation of central masses ((2–7) ×107 M) and mean disk radii (0.2–0.4 pc). Accelerations may ultimately permit estimation of distances when combined with interferometer data. This is notable because the three AGNs are relatively distant (10,000 km s−1 < vsys < 15,000 km s−1), and fractional error in a derived Hubble constant, due to peculiar motion of the galaxies, would be small. As signposts of highly inclined geometries at galactocentric radii of ∼0.1–1 pc, disk masers also provide robust orientation references that allow analysis of (mis)alignment between AGNs and surrounding galactic stellar disks, even without extensive interferometric mapping. We find no preference among published disk maser candidates to lie in high-inclination galaxies. This provides independent support for conclusions that in late-type galaxies, central engine accretion disks and galactic plane orientations are not correlated.

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 (LOS) velocities and sky positions from 18 epochs of very long baseline interferometry observations, and LOS accelerations measured from a 10 yr 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 \pm 0.17 \pm 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, is (H_0 = 72.0 \pm 3.0 km s^{–1} Mpc^{–1}).

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∥,k⊥) 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.

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.

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 104 mK on comoving scales k≲0.5hMpc−1. This represents the first upper limits on the 21 cm power spectrum fluctuations at redshifts 12≲z≲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.