No Arabic abstract
The digital correlator is a crucial element in a modern radio telescope. In this paper we describe a scalable design of the correlator system for the Tianlai pathfinder array, which is an experiment dedicated to test the key technologies for conducting 21cm intensity mapping survey. The correlator is of the FX design, which firstly performs Fast Fourier Transform (FFT) including Polyphase Filter Bank (PFB) computation using a Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) Reconfigurable Open Architecture Computing Hardware-2 (ROACH2) board, then computes cross-correlations using Graphical Processing Units (GPUs). The design has been tested both in laboratory and in actual observation.
The Tianlai Dish Pathfinder Array is a radio interferometer designed to test techniques for 21~cm intensity mapping in the post-reionization universe as a means for measuring large-scale cosmic structure. It performs drift scans of the sky at constant declination. We describe the design, calibration, noise level, and stability of this instrument based on the analysis of about $sim 5 %$ of 6,200 hours of on-sky observations through October, 2019. Beam pattern determinations using drones and the transit of bright sources are in good agreement, and compatible with electromagnetic simulations. Combining all the baselines, we make maps around bright sources and show that the array behaves as expected. A few hundred hours of observations at different declinations have been used to study the array geometry and pointing imperfections, as well as the instrument noise behaviour. We show that the system temperature is below 80~K for most feed antennas, and that noise fluctuations decrease as expected with integration time, at least up to a few hundred seconds. Analysis of long integrations, from 10 nights of observations of the North Celestial Pole, yielded visibilities with amplitudes of 20-30~mK, consistent with the expected signal from the NCP radio sky with $<10,$mK precision for $1 ~mathrm{MHz} times 1~ mathrm{min}$ binning. Hi-pass filtering the spectra to remove smooth spectrum signal yields a residual consistent with zero signal at the $0.5,$mK level.
In 21~cm intensity mapping, the spectral smoothness of the foreground is exploited to separate it from the much weaker 21~cm signal. However, the non-smooth frequency response of the instrument complicates this process. Reflections and standing waves generate modulations on the frequency response. Here we report the analysis of the standing waves in the bandpass of the signal channels of the Tianlai Cylinder Array. By Fourier transforming the bandpass into the delay time domain, we find various standing waves generated on the telescope. A standing wave with time delay at about 142 ns is most clearly identified which is produced in the 15 meter feed cable. We also find a strong peak at a shorter delay of $tau < 50 s$, which may be a mix of the standing wave between the reflector and feed, and the standing wave on the 4 m intermediate frequency (IF) cable. We also show that a smoother frequency response could be partially recovered by removing the reflection-inducted modulations. However, the standing wave on the antenna is direction-dependent, which poses a more difficult challenge for high precision calibration.
In this paper, we apply our sky map reconstruction method for transit type interferometers to the Tianlai cylinder array. The method is based on the spherical harmonic decomposition, and can be applied to cylindrical array as well as dish arrays and we can compute the instrument response, synthesised beam, transfer function and the noise power spectrum. We consider cylinder arrays with feed spacing larger than half wavelength, and as expected, we find that the arrays with regular spacing have grating lobes which produce spurious images in the reconstructed maps. We show that this problem can be overcome, using arrays with different feed spacing on each cylinder. We present the reconstructed maps, and study the performance in terms of noise power spectrum, transfer function and beams for both regular and irregular feed spacing configurations.
Precision measurement of the beam pattern of an antenna is very important for many applications. While traditionally such measurement is often made in a microwave anechoic chamber or at a test range, measurement using an unmanned aerial vehicle offers a number of advantages: the measurement can be made for the assembled antenna on site, thus reflecting the actual characteristics of the antenna of interest, and more importantly, it can be performed for larger antennas which cannot be steered or easily measured using the anechoic chamber and test range. Here we report our beam measurement experiment with UAV for a 6 meter dish used in the Tianlai array, which is a radio astronomy experiment. Due to the dishs small collecting area, calibration with an astronomical source only allows for determining the antenna beam pattern over a very limited angular range. We describe in detail the setup of the experiment, the components of the signal transmitting system, the design of the flight path and the procedure for data processing. We find the UAV measurement of the beam pattern agrees very well with the astronomical source measurement in the main lobe, but the UAV measurement can be extended to the fourth side lobe. The measured position and width of each lobe also shows good agreement with electromagnetic field simulation. This UAV-based approach of beam pattern measurement is flexible and inexpensive, and the technique may also be applied to other experiments.
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.