Fengquan Wu

We present the first blind search for OH 18 cm absorption with the Five-hundred-meter Aperture Spherical Telescope (FAST), conducted alongside the H i 21 cm absorption search. Our previous FAST blind H i absorption search identified 34 systems. In this work, we extend the search using 2024 and part of the 2025 CRAFTS and FAST All Sky H i Survey data (394.4 hr and 1622.1 deg ² ) together with FATHOMER observations, yielding four known and three new H i absorbers, for a total of 41 H i absorption systems. We search for OH absorption in 19 H i absorption systems whose OH redshifted frequencies fall within the FAST band. The known OH absorber toward PKS 1413+135 was redetected, making our survey the first blind survey to detect OH absorption. No new OH absorbers were identified. We examine the relationship between N _OH and N _HI , applying survival analysis to account for upper limits. The analysis does not provide statistically significant evidence for either an N _OH – N _HI correlation or redshift evolution of N _OH / N _HI . Finally, spectral stacking sets 3 σ OH column density upper limits of 4.93, 1.64, and 1.72 T _ex / c _f,OH × 10 ¹² cm ⁻² K ⁻¹ for associated, intervening, and combined samples, corresponding to [OH]/[H i ] ratios of −8 , −8 , and −8 , assuming T _ex = 10 K for OH and T _s = 100 K for H i . These results place the strongest constraints to date on the OH content in radio-selected H i absorbers and establish a blind-survey benchmark for future studies of molecular gas in H i -selected systems. They also demonstrate that known H i 21 cm absorbers provide an effective parent sample for systematic OH absorption searches, paving the way for future larger surveys.

We present the beam pattern measurement of the Tianlai Cylinder Pathfinder Array. As it is a pure drift-scan instrument, we exploit the North-South motion of the Sun to demonstrate that the primary beam is factorizable. Leveraging this property, we decompose the primary beam into independent East-West (E-W) and North-South (N-S) components. Using the Sun as a calibration source, we obtain the E-W beam profiles at various elevations, applying normalization to eliminate the effects of solar activity. Subsequently, we simulate the observed signals using a sky map model to derive the best-fit N-S beam. The results of this work are consistent with previous expectations.

Ground-based radio astronomical observation at frequencies below 30 MHz is hampered by the Ionosphere and radio frequency interference. The Discovering Sky at the Longest wavelength mission, also known as the Hongmeng mission, employs a linear array of satellites on a circular orbit around the Moon to make interferometric observations in this band. Though vastly different from the usual ground-based arrays, the interferometric visibility data collected by such an array are linearly related to the sky map, and the reconstruction is, in principle, an inversion problem of linear mapping. In this paper, we investigate a number of issues in the algorithm of global map reconstruction, focusing on the impact of subpixel noise induced by the finite pixelization of the sky and errors due to regularization. We find that in the reconstruction process, if one builds up the beam matrix, which relates the sky pixels to the visibilities, by naively evaluating its elements at each of the pixel centers, then the subpixel noise can give rise to a significant aliasing effect. However, this effect can be effectively mitigated by a simple pixel-averaging method. Based on the evaluation of the image quality using the correlation coefficient between the input and reconstructed map, and the signal-to-noise ratio, we discuss the selection strategy of the regularization parameter, and show that the sky can be well reconstructed with a reasonable choice of the regularization parameter.

The Tianlai Cylinder Pathfinder Array consists of three adjacent cylindrical reflectors fixed on the ground, each 40 m long and 15 m wide, with the cylinder axis oriented along the North–South (N–S) direction. Dual linear polarization feeds are distributed along the focus line, parallel to the cylinder axis. Measurement of the primary beam profile of these cylindrical reflectors is difficult, as they are too large to be placed in an anechoic chamber. While the beam profile along the East–West (E–W) direction can be measured with the transit observations of bright astronomical radio sources, the beam profile along the N–S direction remains very uncertain. We present a preliminary measurement of the beam profile of the Tianlai cylindrical antenna along both the N–S direction and E–W direction in the frequency range of 700–800 MHz, using a calibrator source carried by an unmanned aerial vehicle (UAV) flying in the far field. The beam profile of the Tianlai cylindrical antenna is determined from the analysis of the auto-correlation signals from the cylinder array correlator, taking into account the emitter antenna beam profile, itself measured with a dipole antenna on the ground. The accuracy of the UAV-based determination of the cylinder beam profiles is validated by comparing the results with the one derived from bright astronomical source transits, and with simulated beams.

For an interferometric array, an image of the sky can be synthesized from interferometric visibilities, which are the cross-correlations of the received electric voltages of pairs of array elements. However, to search for transient targets such as the fast radio burst (FRB), it is more convenient to use the beam-forming technique, where the real-time voltage outputs of the array elements are used to generate data streams (beams) which are sensitive to a specific direction. This is usually achieved by a weighted sum of the array element voltages, with the complex weight adjusted so that all outputs have the same phase for that direction. Alternatively, beams can also be formed from the weighted sum of the short time averaged correlation (visibility) data. We shall call these two approaches the electric voltage beam forming (EBF) and cross-correlation beam forming (XBF), respectively. All beams formed with the EBF can also be formed by the XBF method, but the latter can also generate beams which can not be generated by the former. We discuss the properties of these two kinds of beams, and the amount of computation required in each case. For an array with large number of elements, the XBF would require much more computation resource, although this is partly compensated by the fact that it allows integration over time. We study the impact of cross-coupling between array elements on the beamforming, first using a toy model, then for the case of the Tianlai Cylinder Pathfinder Array. In both cases, we find that the impact of the cross-coupling on the beam profile is relatively small. The understanding gained in this study is helpful in designing and understanding the beam-forming FRB digital backend for compact arrays such as the Tianlai array.

For an interferometric array, an image of the sky can be synthesized from interferometric visibilities, which are the cross-correlations of the received electric voltages of pairs of array elements. However, to search for transient targets such as the fast radio burst (FRB), it is more convenient to use the beam-forming technique, where the real-time voltage outputs of the array elements are used to generate data streams (beams) which are sensitive to a specific direction. This is usually achieved by a weighted sum of the array element voltages, with the complex weight adjusted so that all outputs have the same phase for that direction. Alternatively, beams can also be formed from the weighted sum of the short time averaged correlation (visibility) data. We shall call these two approaches the electric voltage beam forming (EBF) and cross-correlation beam forming (XBF), respectively. All beams formed with the EBF can also be formed by the XBF method, but the latter can also generate beams which can not be generated by the former. We discuss the properties of these two kinds of beams, and the amount of computation required in each case. For an array with large number of elements, the XBF would require much more computation resource, although this is partly compensated by the fact that it allows integration over time. We study the impact of cross-coupling between array elements on the beamforming, first using a toy model, then for the case of the Tianlai Cylinder Pathfinder Array. In both cases, we find that the impact of the cross-coupling on the beam profile is relatively small. The understanding gained in this study is helpful in designing and understanding the beam-forming FRB digital backend for compact arrays such as the Tianlai array.

Neutral hydrogen (H i) intensity mapping (IM) presents great promise for future cosmological large-scale structure surveys. However, a major challenge for H i IM cosmological studies is to accurately subtract the foreground contamination. An accurate beam model is crucial for improving the quality of foreground subtraction. In this work, we develop a stacking-based beam reconstruction method utilizing the radio continuum point sources within the drift-scan field. Based on the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we employ two sets of drift-scan survey data and merge the measurements to construct the beam patterns of the 19 FAST L-band feeds. To model the beams, we utilize the Zernike polynomial, which effectively captures asymmetric features of the main beam and the different side lobes. Due to the symmetric location of the beams, the main features of the beams are closely related to the distance from the center of the feed array, e.g., as the distance increases, side lobes become more pronounced. This modeling pipeline leverages the stable drift-scan data to extract beam patterns while accounting for and excluding the reflector’s changing effects. It provides a more accurate measurement beam and a more precise model beam for FAST H i IM cosmology surveys.

The cosmic 21 cm signal serves as a crucial probe for studying the evolutionary history of the Universe. However, detecting the 21 cm signal poses significant challenges due to its extremely faint nature. To mitigate the interference from the Earth’s radio frequency interference (RFI), the ground and the ionospheric effects, the Discovering the Sky at the Longest Wavelength (DSL) project will deploy a constellation of satellites in lunar orbit, with its high-frequency daughter satellite tasked with detecting the global 21 cm signal from cosmic dawn and reionization era (CD/EoR). We intend to employ the vari-zeroth-order polynomial (VZOP) for foreground fitting and subtracting. We have studied the effect of thermal noise, thermal radiation from the Moon, the lunar reflection, anisotropic frequency-dependent beam, inaccurate antenna beam pattern, and RFI contamination. We discovered that the RFI contamination can significantly affect the fitting process and thus prevent us from detecting the signal. Therefore, experimenting on the far side of the moon is crucial. We also discovered that using VZOP together with DSL, after 1080 orbits around the Moon, which takes about 103 days, we can successfully detect the CD/EoR 21 cm signal.

We present an extragalactic H i 21 cm absorption lines catalog from a blind search at z ≤ 0.35, using drift-scan data collected in 1325.6 hr by the ongoing Commensal Radio Astronomy Fast Survey and FAST All Sky H i Survey, which spans a sky area of 6072.0 deg² and covers 84,533 radio sources with a flux density greater than 12 mJy. Fourteen previously identified H i absorbers and 20 newly discovered H i absorbers were detected, comprising 15 associated systems, 10 intervening systems, and nine systems with undetermined classifications. Through spectral stacking, the mean peak optical path, mean velocity-integrated optical path, mean FWHM, and mean H i column density are measured to be 0.47 and 0.30; 27.19 and 4.36 km s⁻¹; 42.61 and 9.33 km s⁻¹; 0.49 and 0.08 T _s × 10²⁰ cm⁻² K⁻¹, for the associated and intervening samples, respectively. Statistical analysis also reveals that associated systems tend to be hosted by red (g − r > 0.7) galaxies at lower redshifts, whereas galaxies hosting intervening H i absorption are typically found at higher redshifts and are of a bluer (g − r ≤ 0.7) type. A noticeable difference is observed in the positions of foregrounds, backgrounds of intervening systems, and high-redshift and low-redshift associated systems on the Wide-field Infrared Survey Explorer color–color diagram. All identified foreground sources in our sample have W1 – W2 magnitudes below 0.8, suggesting no active galactic nuclei (AGNs). In contrast, backgrounds of intervening systems tend to have W1 – W2 magnitudes above 0.8, indicating AGN presence. For associated absorption, most low-redshift (z ≤ 0.5) systems show W1 – W2 values below 0.8, while higher-redshift associated absorption (z > 0.5) displays a broader range of W1 − W2 values.

This paper presents the design, calibration, and survey strategy of the Fast Radio Burst (FRB) digital backend and its real-time data processing pipeline employed in the Tianlai Cylinder Pathfinder Array. The array, consisting of three parallel cylindrical reflectors and equipped with 96 dual-polarization feeds, is a radio interferometer array designed for conducting drift scans of the northern celestial semi-sphere. The FRB digital backend enables the formation of 96 digital beams, effectively covering an area of approximately 40 square degrees with the 3 dB beam. Our pipeline demonstrates the capability to conduct an automatic search of FRBs, detecting at quasi-real-time and classifying FRB candidates automatically. The current FRB searching pipeline has an overall recall rate of 88%. During the commissioning phase, we successfully detected signals emitted by four well-known pulsars: PSR B0329+54, B2021+51, B0823+26, and B2020+28. We report the first discovery of an FRB by our array, designated as FRB 20220414A. We also investigate the optimal arrangement for the digitally formed beams to achieve maximum detection rate by numerical simulation.