The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

We explore the cosmological multi-tracer synergies between a DESI-like Emission Line Galaxy (ELG) distribution and a Tianlai-like 21cm intensity map. We use simulated maps generated from a particle simulation in the light-cone volume (Horizon Run 4), sky-trimmed and including the effects of foreground contamination, its removal, and instrument noise. We first validate how the foreground residual affects the recovered 21cm signal by putting different levels of foreground contamination into the 21cm maps. We find that the contamination cannot be ignored in the angular auto-power spectra of HI even when it is small, but has no influence on the accuracy of the angular cross-power spectra between HI and galaxies. In the foreground-cleaned map case, as information is lost in the cleaning procedure, there is also a bias in the cross-power spectrum. However, we found that the bias from the cross-power spectrum is scale-independent, which is easily parameterised as part of the model, while the offset in the HI auto-power spectrum is non-linear. In particular, we tested that the cross power also benefits from the cancellation of the bias in the power spectrum measurement that is induced by the instrument noise, which changes the shape of the auto-power spectra but leaves the cross-power unaffected. We then modelled the angular cross-power spectra to fit the BAO feature in broadband shape of the angular cross power spectrum, including contamination from the residual foreground and the effect of instrument noise. We forecast a constraint on the angular diameter distance $D_\mathrm{A}$ for the Tianlai Pathfinder redshift $0.775<z<1.03$, giving a distance measurement with a precision of 3.3% at that redshift.

We present a simple and fast method for incoherent dedispersion and fast radio burst (FRB) detection based on the Hough transform, which is widely used for feature extraction in image analysis. The Hough transform maps a point in the time-frequency data to a straight line in the parameter space, and points on the same dispersed $f^{-2}$ curve to a bundle of lines all crossing at the same point, thus the curve is transformed to a single point in the parameter space, enabling an easier way for the detection of radio burst. By choosing an appropriate truncation threshold, in a reasonably radio quiet environment, i.e. with radio frequency interferences (RFIs) present but not dominant, the computing speed of the method is very fast. Using simulation data of different noise levels, we studied how the detected peak varies with different truncation thresholds. We also tested the method with some real pulsar and FRB data.

The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate HI intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the HI sky for the frequency range 700-800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21cm distribution on large-scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information.

We propose an eigenvector-based formalism for the calibration of radio interferometer arrays. In the presence of a strong dominant point source, the complex gains of the array can be obtained by taking the first eigenvector of the visibility matrix. We use the stable principle component analysis (SPCA) method to help separate outliers and noise from the calibrator signal to improve the performance of the method. This method can be applied with poorly known beam model of the antenna, and is insensitive to outliers or imperfections in the data, and has low computational complexity. It thus is particularly suitable for the initial calibration of the array, which can serve as the initial point for more accurate calibrations. We demonstrate this method by applying it to the cylinder pathfinder of the Tianlai experiment, which aims to measure the dark energy equation of state using the baryon acoustic oscillation (BAO) features in the large scale structure by making intensity mapping observation of the redshifted 21~cm emission of the neutral hydrogen (HI). The complex gain of the array elements and the beam profile in the East-West direction (short axis of the cylinder) are successfully obtained by applying this method to the transit data of bright radio sources.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

Radio astronomical observation below 30 MHz is hampered by the refraction and absorption of the ionosphere, and the radio frequency interference (RFI), so far high angular resolution sky intensity map is not available. An interferometer array on lunar orbit provides a perfect observatory in this frequency band: it is out of the ionosphere and the Moon helps to block the RFIs from the Earth. The satellites can make observations on the far side of the Moon and then send back the data on the near side part of the orbit. However, for such array the traditional imaging algorithm is not applicable: the field of view is very wide (almost whole sky), and for baselines distributed on a plane, there is a mirror symmetry between the two sides of the plane. A further complication is that for each baseline, the Moon blocks part of the sky, but as the satellites orbit the Moon, both the direction of the baseline and the blocked sky change, so even imaging algorithms which can deal with non-coplanar baseline may not work in this case. Here we present an imaging algorithm based on solving the linear mapping equations relating the sky intensity to the visibilities. We show that the mirror symmetry can be broken by the three dimensional baseline distribution generated naturally by the precession of the orbital plane of the satellites. The algorithm is applicable and good maps could be reconstructed, even though for each baseline the sky blocking by the Moon is different. We also investigate how the map-making is affected by inhomogeneous baseline distributions.

The redshifted 21~cm signal from neutral hydrogen (HI) is potentially a very powerful probe for cosmology, but a difficulty in its observation is that it is much weaker than foreground radiation from the Milky Way as well as extragalactic radio sources. The foreground radiation at different frequencies are however coherent along one line of sight, and various methods of foreground subtraction based on this property have been proposed. In this paper, we present a new method based on the Robust Principal Component Analysis (RPCA) to subtract foreground and extract 21~cm signal, which explicitly uses both the low-rank property of the frequency covariance matrix (i.e. frequency coherence) of the foreground and the sparsity of the frequency covariance matrix of the 21~cm signal. The low-rank property of the foregrounds frequency covariance has been exploited in many previous works on foreground subtraction, but to our knowledge the sparsity of the frequency covariance of the 21~cm signal is first explored here. By exploiting both properties in the RPCA method, in principle, the foreground and signal may be separated without the signal loss problem. Our method is applicable to both small patch of sky with the flat-sky approximation, and to large area of sky where the sphericity has to be considered. It is also easy to be extended to deal with more complex conditions such as sky map with defects.

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.