We present a deep H I survey at L band conducted with the Fivehundredmetre Aperture Spherical radio Telescope FAST over the COSMOS field. This survey is strategically designed to overlap with the MIGHTEE COSMOS field aiming to combine the sensitivity of the FAST and highresolution of the MeerKAT. We observed the field with FAST for approximately 11 hours covering inlineformulatexmath idTM0001 notationLaTeXsim 2texmathinlineformula square degrees and reduced the raw data to H I spectral cubes over the frequency range 13101420 MHz. The FASTH I data reach a median 3inlineformulatexmath idTM0002 notationLaTeXsigmatexmathinlineformula column density of inlineformulatexmath idTM0003 notationLaTeXNrm Hsmall I sim 2times 1017texmathinlineformula cminlineformulatexmath idTM0004 notationLaTeX2texmathinlineformula over a inlineformulatexmath idTM0005 notationLaTeXsim 5 rm km rm s1texmathinlineformula channel width allowing for studies of the distribution of H I gas in various environments such as in galaxies the CircumGalactic Medium CGM and Intergalactic Medium IGM. We visually searched the spectral cubes for H I sources and found a total of 80 H I detections of which 56 have been crossmatched with the MIGHTEEH I catalogue. With the crossmatched sources we compare their H I masses and find that the total H I mass fraction in the IGM and CGM surrounding the galaxy pairs is statistically higher than the H I fraction surrounding the isolated galaxies by a difference of inlineformulatexmath idTM0006 notationLaTeXsim 13pm 4texmathinlineformula per cent indicating that the CGM and IGM associated with interacting systems are richer in neutral hydrogen compared to those around isolated galaxies in the local Universe. We also describe several FASTMeerKAT synergy projects highlighting the full potential of exploiting both singledish and interferometric observations to study the distribution and evolution of the diffuse H I gas.
['Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
Deep extragalactic H I survey of the COSMOS field with FAST
{'Deep extragalactic H /i.pc survey of the COSMOS field with FAST': "Hengxing Pan 1 ★ , Matt J. Jarvis 1 , 2 , Ming Zhu 3 , 4 , Yin-Zhe Ma 5 , Mario G. Santos 2 , 6 , Anastasia A. Ponomareva 1 ,  \nIan Heywood 1 , 7 , 6 , Yingjie Jing 4 , Chen Xu 4 , 8 , Ziming Liu 4 , 8 , Yogesh Chandola 9 and Yipeng Jing 10 , 11  \n- 1 Astrophysics, University of Oxford, Denys Wilkinson Buiding, Keble Road, Oxford OX1 3RH, UK\n- 2 Department of Physics and Astronomy, University of the Western Cape, Cape Town 7535, South Africa\n- 3 Guizhou Radio Astronomical Observatory, Guizhou University, Guiyang 550000, People's Republic of China\n- 4 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, People's Republic of China\n- 5 Department of Physics, Stellenbosch University, Matieland 7602, South Africa\n- 6 South African Radio Astronomy Observatory (SARAO), 2 Fir Street, Observatory, 7925, South Africa\n- 7\n- Centre for Radio Astronomy Techniques and Technologies, Department of Physics and Electronics, Rhodes University, PO Box 94, Makhanda 6140, South Africa\n- 8 School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China\n- 9 Purple Mountain Observatory, Chinese Academy of Sciences, 10, Yuan Hua Road, Qixia District, Nanjing, 210023, China\n- 10 Department of Astronomy, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China\n- 11 Tsung-Dao Lee Institute, and Shanghai Key Laboratory for Particle Physics and Cosmology,\n- Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China  \nAccepted XXX. Received YYY; in original form ZZZ", 'ABSTRACT': 'We present a deep H /i.pc survey at L-band conducted with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) over the COSMOS field. This survey is strategically designed to overlap with the MIGHTEE COSMOS field, aiming to combine the sensitivity of the FAST and high-resolution of the MeerKAT. We observed the field with FAST for approximately 11 hours covering ∼ 2 square degrees, and reduced the raw data to H /i.pc spectral cubes over the frequency range 1310-1420 MHz. The FAST-H/i.pc data reach a median 3 𝜎 column density of 𝑁 HI ∼ 2 × 10 17 cm -2 over a ∼ 5kms -1 channel width, allowing for studies of the distribution of H /i.pc gas in various environments, such as in galaxies, the Circum-Galactic Medium (CGM) and Intergalactic Medium (IGM). We visually searched the spectral cubes for H /i.pc sources, and found a total of 80 H /i.pc detections, of which 56 have been cross-matched with the MIGHTEE-H /i.pc catalogue. With the cross-matched sources, we compare their H /i.pc masses and find that the total H /i.pc mass fraction in the IGM and CGM surrounding the galaxy pairs is statistically higher than the H /i.pc fraction surrounding the isolated galaxies by a difference of ∼ 13 ± 4%, indicating that the CGM and IGM associated with interacting systems are richer in neutral hydrogen compared to those around isolated galaxies in the local Universe. We also describe several FAST-MeerKAT synergy projects, highlighting the full potential of exploiting both single-dish and interferometric observations to study the distribution and evolution of the diffuse H /i.pc gas.  \nKey words: galaxies: intergalactic medium - radio lines: galaxies', '1 INTRODUCTION': "Galaxy growth relies on the continuous gas accretion from the surrounding environments. Cold gas within galaxies in the form of neutral hydrogen (H /i.pc) serves as the reservoir from which molecules and stars subsequently form. The Intergalactic Medium (IGM) connects galaxies and acts as a route to fueling them, and the Circum-Galactic Medium (CGM) traces the flow of matter, energy, and enriched elements between galaxies and the IGM (e.g. Tumlinson et al. 2017; Martin 2019; Saintonge & Catinella 2022; Decataldo et al. 2024). Understanding the distribution of baryons within galaxies, the CGM  \nand IGM is therefore a fundamental problem in astrophysics, with critical implications for the growth mechanisms of galaxies (e.g. Sancisi et al. 2008; Putman 2017; Faucher-Gigu'ere & Oh 2023) and can help solve the long-lasting 'missing' baryon problem (Shull et al. 2012; Ma et al. 2015; Macquart et al. 2020; Li et al. 2024).  \nIn the past few decades, several H /i.pc galaxy surveys with large sky coverage have been undertaken by radio telescopes, such as the H/i.pc Parkes All-Sky (HIPASS) Survey (Barnes et al. 2001) and the Arecibo Legacy Fast ALFA (ALFALFA) survey (Giovanelli et al. 2005), to detect the 21-cm emission line from neutral hydrogen in galaxies. However, these surveys are constrained to the local Universe ( 𝑧 ≲ 0 . 06), due to the limited sensitivity and frequency range of the telescopes. On the other hand, our view of the IGM and CGM is  \nbased largely on the powerful but restricted information from QSO absorption line surveys that probe the Lyman𝛼 forest (e.g. Becker et al. 2011; Keating et al. 2013; Telikova et al. 2019). Unfortunately, emission from intergalactic baryons is difficult to observe because of current telescope sensitivities, which limit studies to column densities 𝑁 HI ≳ 10 19 (atoms) cm -2 , which is the realm of Damped Ly 𝛼 (DLA) systems and sub-DLAs. Below column densities of ∼ 10 19 cm -2 , the neutral fraction of hydrogen decreases rapidly due to the transition from optically-thick to optically-thin gas ionized by the metagalactic ultraviolet radiation. However, below ∼ 10 18 cm -2 the gas is optically thin and the decline in neutral fraction is much more gradual (Popping et al. 2009; Popping et al. 2015). This gas, generally thought to be residing in filamentary structures, serves as the the pristine reservoir capable of fueling future star formation in galaxies, and could provide a direct signature of the predicted smooth cold-mode accretion to dominate gas acquisition in star-forming galaxies today (Kereˇs et al. 2009; de Blok et al. 2016). Hence, exploring the 𝑁 HI < 10 18 cm -2 regime is crucial to deepen our understanding of the reservoir of baryons that shape galaxy formation and evolution.  \nWith the sensitivity of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST; Nan et al. 2011), we now have a facility with the potential to detect these faint H /i.pc emissions within a reasonable integration time, as indicated, for example, by the FAST All Sky HI survey (FASHI; Zhang et al. 2024). However, due to the large beam size ( ∼ 3 arcmin) of FAST, disentangling the IGM/CGM from the H /i.pc gas within individual galaxies is limited to the very local nearby Universe. Fortunately, with the excellent angular resolution ( ∼ 10 arcsec) of the MeerKAT radio telescope (Jonas & MeerKAT Team 2016), we can identify the H /i.pc within these galaxies directly, and then use a combination of data from FAST and MeerKAT to extract the IGM/CGM signal at the 𝑁 HI < 10 18 cm -2 regime thoroughly. To achieve this, we carried out a series of FAST observations in L-band across the COSMOSfield, overlapping with one of the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE; Jarvis et al. 2016) fields.  \nIn this paper, we introduce our deep H /i.pc survey aimed at detecting the H /i.pc gas not only in the galaxies but also in the IGM and CGM with FAST. We describe the data processing for FAST and the ancillary data in Section 2. We then present the FAST-H /i.pc detections, and investigate the 𝑀 HI -𝑀 ★ relation between the FAST-H /i.pc and the MIGHTEE-H/i.pc samples in Section 3 ending with a list of FASTMeerKAT synergy projects, and then conclude in Section 4. We use the standard Λ CDM cosmology with a Hubble constant 𝐻 0 = 67 . 4 kms -1 Mpc -1 , total matter density Ω m = 0 . 315 and dark energy density Ω Λ = 0 . 685 (Planck Collaboration et al. 2020) in our analysis.", '2.1 FAST-H /i.pc': "FAST has a 500-meter diameter dish constructed in a natural depression in Guizhou province of China and is the largest single-dish telescope in the world. We made observations between the 1st and 6th of January 2022 for six nights over ∼ 2 square degrees of the COSMOS field with FAST's MultiBeamOTF (On-The-Fly) mode for a total amount of 11-hours integration time on source. The footprints are shown in Figure 1. The 19 beams are rotated clockwise by a angle of 23.4 · to ensure a uniform coverage of the field in this horizontal scanning mode. The scan gap is fixed at 3.5 arcmins to reduce the frequency of the source changing. The flux calibrator 3C237 was observed on the first day of the observations with FAST's MultiBeamCalibration mode. The observations are listed in Table 1.  \nFigure 1. Footprints of the FAST 19 beams in the COSMOS field. Each beam size is 2.9 ' at 𝑧 = 0. The arrows and the blue line show the direction of the beam 01 (i.e. the centre beam denoted by a red circle) scanning the field with the MultibeamOTF mode.  \n<!-- image -->  \nWe use the HiFAST 1 pipeline (Jing et al. 2024) to reduce the raw spectral data recorded by the 19-beam L-band receiver of FAST. This pipeline is a dedicated, modular, and self-contained calibration and imaging system designed for processing the H /i.pc data from FAST. The modules in the pipeline can be combined as needed to process the data from various FAST observation modes, including tracking, drift scanning, On-The-Fly mapping, and most of their variants. We provide information on a few key modules in the following subsections.", '2.1.1 Calibration': 'Wefirst calibrate the antenna temperature with a standard 10 K noise diode, and the noise source was injected every 32 seconds. However, the gain fluctuations depend on the condition of the telescope and the receiver during observation. Therefore, we observed the calibrator 3C237 during the initial stage of observing the target sources to calibrate the flux density scale. We use the latest continuum measurements of the calibrator 3C237 as the model across the frequency range of 900-1670 MHz from MeerKAT. The flux density as a function of frequency is described by a power law,  \n𝑆 ( 𝑣 ) = 𝑆 1 . 4GHz GLYPH<16> 𝑣 1 . 4 GHz GLYPH<17> 𝛼 , (1)  \nwhere 𝑆 1 . 4GHz = 6 . 533 ± 0 . 147 Jy, and 𝛼 = -0 . 919 ± 0 . 01 (Taylor & Legodi 2021). The uncertainty on the flux density is ∼ 2 percent at the frequency of ∼ 1.4 GHz, which is taken into account when the error of the H /i.pc flux is determined in addition to the effect of the thermal noise. The data processing of the calibrator observation is detailed by Liu et al. (2024).', '2.1.2 Baseline subtraction': 'We then fit the spectral baseline with two approaches: 1) the asymmetrically reweighted penalized least squares smoothing (arPLS)  \nTable 1. Observation log for the COSMOS field and the calibrator 3C237.  \nmethod (Baek et al. 2015) and 2) the minimum of the medians (MinMed) method (Putman et al. 2002). The arPLS method iteratively refines the baseline by adjusting and applying weights to the data points until the weights converge. The MinMed method splits the scanned spectra into several parts and each part comprises several segments along the time axis for each channel, and then the minimum of the medians in each segment is taken as the reference (off-source) spectra which is subtracted in that part afterwards. These two approaches complement each other to account for continuum variations across a broad range of the frequency scales.', '2.1.3 Standing wave removal': "Standing waves arise because the radio signal enters the receiving system by two paths of different length due to the present design of radio telescopes. The signal taking the longer path suffers a delay, and then correlates with the signal in the shorter path as a function of time. To obtain the power spectrum, the correlation function is Fourier-transformed to a sinusoidal variation across the frequency band. The difference of the two path lengths is ∼ 276 meters for FAST, which corresponds to a time delay of ∼ 0.92 microseconds, therefore the resulting standing wave has a period of ∼ 1.09 MHz in frequency space. Hence we fit the standing wave with a Fourier approach (Xu et al. in prep). This approach involves applying a Fourier transform to the 'waterfall data' (i.e. spectral data as a function of time and frequency) in the frequency space, transforming them into the delay space. We then select the peak amplitude and its surrounding modes for performing an inverse Fourier transform at each timestamp. The standing wave is modelled in this way, and then removed from the baseline subtracted spectra.", '2.1.4 RFI flagging': 'The radio frequency interference (RFI) often refers to artificial radio emission from human-made electronics such as satellites and civil aircraft. It corrupts astronomical measurements and cannot be well modelled, therefore has to be masked out. The RFI presented in our data can be loosely classified as three types: the time domain RFI that cover a wide frequency band from geosynchronous satellites, the narrow band single channel RFI, and polarized RFI. We identify the time domain RFI by averaging the 2-dimensional waterfall data along the frequency axis, and select the narrow band single channel RFI by averaging the waterfall data along the time axis. The polarized RFI is flagged as the H /i.pc signal tends to be non-polarized. We also perform a visual checkup to manually mask out the RFI that is not well defined by the above types. Overall, 10-20 percent of the data is masked out across the bandwidth of 1310-1420 MHz, with most of the affected data due to the wide band RFI. The RFI contamination at the lower frequency ( < 1310 MHz) band is complicated (Zhang et al. 2022) and we leave it for a separate paper to describe the flagging procedures.  \nTable 2. Key observational characteristics of the FAST-H /i.pc COSMOS data. Note that the total sky coverage is ∼ 2 deg 2 where the centre sensitive area is ∼ 1.2 deg 2 .', '2.1.5 Gridding and smoothing': 'We correct the coordinate frame for the Doppler effect due to the motion of the earth, and finally grid the individual spectra into a 3-dimensional cube with a weighting scheme that accounts for the distance of each spectrum from the grid point (Mangum et al. 2007). The channel width of raw FAST data is 7.63 kHz, and we smooth the H /i.pc data cube with a Hanning window along the frequency to a resolution of 22.9 kHz which is ∼ 4 . 83kms -1 at 𝑧 = 0. At this stage, the data cube still has a low-level continuum residual, which has been further subtracted by fitting a second-degree polynomial function to the residual per pixel after applying a 3 sigma-clipping. We note that a H /i.pc cube with a lower-velocity resolution of 26 kms -1 is also made for the visual source finding only.', '2.1.6 Data cube': 'We show the RMS against frequency in Figure 2, and show that our final H /i.pc cube has a median channel noise of ∼ 385 𝜇 Jy beam -1 , which roughly meets our expectation considering that the COSMOS field is close to the equator where the aperture efficiency of FAST is reduced by ∼ 10%with a zenith angle of ∼ 30 · (Jiang et al. 2020), and the additional effects of flagging the RFI. The noise level is relatively constant across the whole band between 1310 and 1420 MHz with a few pronounced spikes due to the flagging of the narrow band single channel RFI. However, the values from the median absolute deviation (MAD) is ∼ 20% lower than the standard deviation. This means that the noise does not follow a perfect Gaussian distribution in the spatial domain as we can see from the color-coded noise map in Figure 3, although the areas outside the dashed orange box are excluded. Indeed, our observations suffered from the strong wide band RFI from a Geosynchronous satellite and lost about half of the integration time in the upper spectral coverage of the COSMOS field, likely leading to suboptimal continuum subtractions and a moderate level of nonGaussian noise behaviour. We provide a few key parameters for our FAST data in Table 2.  \nFigure 2. Noise vs frequency for the H /i.pc cube with a channel width of ∼ 23 kHz. The blue and orange lines are the noise levels estimated by using the standard deviation and median absolute deviation (MAD) for each channel, respectively. The pronounced noise spikes are due to the flagging of the single channel RFI.  \n<!-- image -->  \nFigure 3. FAST observation coverage (in color-coded area) overlapped over MIGHTEE Early Science COSMOS field (in black circle). The blue and orange stars are the MIGHTEE-H /i.pc and FAST-H /i.pc detections, respectively. The colour scheme indicates the noise level on a projected 2-dimensional sky. The dashed orange box delineates the perimeter within which the centre beam of the 19 beams has passed over, and the noise level increases gradually towards the boundaries of the whole field due to limited number of beams available. The black circle shows the MeerKAT primary beam at ∼ 10% level.  \n<!-- image -->', '2.2 Ancillary data': 'The COSMOS field is covered by various multi-wavelength photometric and spectroscopic surveys ranging from X-ray to radio bands. In particular, we exploit the MIGHTEE-H /i.pc data which is the H /i.pc emission project within the MIGHTEE survey undertaken by the interferometric MeerKAT telescope as one of eight large survey projects (Jarvis et al. 2016; Maddox et al. 2021). We make use of the H /i.pc galaxy catalogue compiled from the MIGHTEE-H /i.pc Early Science data which are collected with the MeerKAT 4k correlator mode in the L-band with a frequency range from 900 to 1670 MHz (Ponomareva et al. 2023).  \nHowever, we use the latest MIGHTEE Data Release 1 which is observed with the 32k channels (Heywood et al. 2024) for compar-  \ning the H /i.pc correlator mode. We smooth a cubelet surrounding each source in the Early Science catalogue provided by Ponomareva et al. (2023), and clip it at a 3 𝜎 threshold as a mask for removing the noise following (Ponomareva et al. 2021). The MeerKAT 32k correlator mode has 32,768 channels with a spectral resolution of 26.1 kHz (i.e. 5.5 kms -1 at 1420 MHz). The COSMOS field was observed in this mode for a total of 15 × 8h tracks in a tightly-dithered mosaic that spans ∼ 2 deg 2 , and each pointing is imaged with robustness parameters of 0.0 and 0.5 (Briggs 1995) following up on two rounds of RFI flagging, self-calibration, and the visibility-domain continuum subtraction. Then all pointings are brought together for homogenisation, mosaicking, and a further image-domain continuum subtraction (Heywood et al. 2024; Jarvis et al. 2024). We use the data product with the robustness of 0.5 (i.e. angular resolution ∼ 16 arcseconds across 1300-1420 MHz) to have a well balanced combination of the image resolution and sensitivity for detecting the H /i.pc gas within galaxies. We note that a robust catalogue is currently not available for the 32k data, which is why we use the Early Science catalogue to match to the galaxy positions.  \nThe stellar masses of the H /i.pc galaxies are derived using the Spectral Energy Distribution (SED) fitting code LePhare (Ilbert et al. 2006), and the uncertainty on the stellar mass is conservatively assumed to be ∼ 0.1 dex, due to assumptions made in the SED fitting process on star formation history and initial mass function etc. (Adams et al. 2021). The optical spectroscopic data are mainly taken from the DeepExtragalactic VIsible Legacy Survey (DEVILS) survey (Davies et al. 2018; Hashemizadeh et al. 2021) for cross-matching the optical counterparts with our FAST-H /i.pc detections. The ALFALFA H /i.pc data are from Haynes et al. (2018).', '3.1 FAST-H /i.pc source finding': "We find a total of 80 sources after eyeballing the H /i.pc images with the Cube Analysis and Rendering Tool for Astronomy (CARTA; Comrie et al. 2021), and then employ SoFiA 2 (Serra et al. 2015; Westmeier et al. 2021) to refine the source finding and parameterise the H /i.pc source in subcubes centered at the visually-identified positions. The subcubes have an angular size ranging from 3 to 10 beam sizes (i.e. ∼ 9-30 arcmins) based on the source size. SoFiA 2 is a fully automated 3D source finding pipeline for extragalactic H /i.pc surveys.  \n<!-- image -->  \nFigure 4. FAST-H/i.pc mass as a function of redshift (left) and FAST-H /i.pc flux against the velocity width 𝑊 50 (right) for the cross-matched (blue) and non crossmatched (orange) samples between FAST- and MIGHTEE-H /i.pc detections.  \n<!-- image -->  \nIt first convolves the H /i.pc cube or subcube with various smoothing kernels, and selects voxels with absolute values above a predefined threshold, then links the selected voxels together if they are close enough, to form the mask of the source. For our FAST data, the flux detection threshold is 5 𝜎 , where 𝜎 is the local noise level. We chose this relatively high threshold to mitigate the effect of generally less stable baselines of single-dish data and demonstrate that we can recover the intrinsic fluxes with a median difference fraction of ∼ 1 percent with injected fake sources (see Appendix A). The kernel has an angular size of 0, 3 and 6 pixels (i.e. ∼ 2 beam sizes maximum as each pixel size is 1 ' ) in the projected sky direction, and a velocity width of 0, 3, 7, 15, and 31 channels (i.e. 155 kms -1 maximum).  \nFor each source, the integrated flux 𝑆 is measured by summing the flux density values of all pixels contained in the source mask, and multiplying by the spectral channel width, and then dividing by the beam solid angle. The uncertainty of the integrated flux is estimated by assuming Gaussian error propagation of the local noise and correcting for the correlation of the spatial pixels due to the finite beam size. The H /i.pc mass is determined under the optically thin gas assumption by  \n𝑀 HI 𝑀 ⊙ ≃ 49 . 7 GLYPH<18> 𝐷 𝐿 Mpc GLYPH<19> 2 GLYPH<18> 𝑆 JyHz GLYPH<19> , (2)  \nwhere 𝑀 HI is the H /i.pc mass, and 𝐷 𝐿 is the luminosity distance (Meyer et al. 2017).  \nThe line width 𝑊 50 of the integrated spectral profile is measured by moving inwards from both ends of the spectrum until the signal exceeds 50% of the peak flux density in the spectrum, and the error of the 𝑊 50 is typically less than 10 percent for sources with the SNR > 10, and 5 percent is adopted for our data if the SNR is greater than 20 (Westmeier et al. 2021).", '3.2.1 Cross matching': "After cross-matching the sample of 80 FAST-H /i.pc detections with the MIGHTEE-H/i.pc Early Science catalogue (Ponomareva et al. 2023),  \nbut using the measurements from Date Release 1 (Heywood et al. 2024), by setting the maximum angular separation of 3 arcminites and velocity offset of ∼ 200 kms -1 , we find 56 cross-matched sources. A comparison of their spatial distribution is shown in Figure 3. The rate of finding a MIGHTEE-H /i.pc counterpart is clearly higher in the centre region with the highest sensitivity. We plot their H /i.pc masses against redshift in the left panel of Figure 4, and the H /i.pc flux against 𝑊 50 in the right panel to investigate the differences between the FAST-H /i.pc and MIGHTEE-H/i.pc detections. It shows that the non MIGHTEEdetected sources are dominated by the low-flux or narrow-line width H/i.pc galaxies as would be expected given the coarse spectrum resolution of 44 kms -1 for the MIGHTEE-H /i.pc Early Science data and the limited sensitivity. However, there are a few non cross-matched H/i.pc sources that have a moderate amount of H /i.pc gas and a velocity width larger than 150 kms -1 in the FAST data. In other words, some non cross-matched H /i.pc sources appear to be relatively bright galaxies, and are detected in FAST-H /i.pc but not in the MIGHTEE-H /i.pc 4k catalogue. This difference is likely indicating the fundamental difference between MeerKAT and FAST as an interferometer is generally less sensitive to extended/diffuse emission while a single dish telescope is sensitive to both compact and extended emission. The source finding with the MIGHTEE-H /i.pc 32k data will help further examine these sources. There are also some sources that are detected by MeerKAT but not FAST from our initial visual source finding. The presence of systematics in the FAST data complicates baseline subtraction, potentially contributing to the 'missing sources' from the FAST data, especially for the faint ones. By using the positions from the MIGHTEE-H/i.pc detections, we extract the FAST-H /i.pc fluxes for the missing ones and find that these non-detections are indeed very faint with signal to noise ratios less than ∼ 10. Their spectra and masses are shown in Appendix B.  \nThe cross-matched detections allow us to study the H /i.pc galaxies and their surrounding CGM and IGM in detail by comparing the FASTand MIGHTEE-H/i.pc images, and also further help to investigate the Tully-Fisher relation, 𝑀 HI -size relation, 𝑀 HI -𝑀 ★ relation, and the H /i.pc mass function that have been investigated separately by the MIGHTEE-H/i.pc team (e.g. Ponomareva et al. 2021, 2023; Rajohnson et al. 2022; Pan et al. 2023). An example of a FAST-H /i.pc detection  \nFigure 5. Example of the FAST-H /i.pc detection, with the blue halo of the interpolated H /i.pc moment-0 map overlaid on DECaLS (Dey et al. 2019) false color image, where the contour (light blue line) with a column density of 4.6 × 10 17 cm -2 is ∼ 1.5 times the local noise level on the moment-0 map. The red circles denote the MIGHTEE-H /i.pc detections with their names labelled. The hatched grey circle on the upper right corner shows the FAST beam size of 2.9 arcmins.  \n<!-- image -->  \nis shown in Figure 5, where a thin H /i.pc bridge appears to be between the two nearby galaxies at 𝑧 ∼ 0 . 0068 with a column density of 𝑁 HI < 10 18 cm -2 possibly due to the tidal interaction given their close distance of ∼ 100 kpc. However, we note that the lowest contour is ∼ 1.5 times the local noise level on the moment-0 map.  \nThe red circles are the MIGHTEE-H/i.pc counterparts J100153.8+022449 and J100128.0+022025 with log 10 ( 𝑀 HI / 𝑀 ⊙ ) = 7 . 66 ± 0 . 02 and 7.23 ± 0 . 02, respectively. On the other hand, the total logarithmic FAST-H /i.pc mass of both galaxies is 7.97 ± 0 . 03. Compared to the total MIGHTEE-H /i.pc mass, the FAST-H /i.pc excess fraction for this pair of galaxies is 32 ± 6 percent, which is in line with the median value of 30 percent for the excess H /i.pc in compact galaxy groups with the Green Bank Telescope (GBT) and the Very Large Array (VLA) data from Borthakur et al. (2010). This excess of H /i.pc gas is presumably outside the galaxies and may be considered to be in the CGM and IGM. We will return to this in the next section.  \nWe compare the H /i.pc masses between the cross-matched FAST and MIGHTEE-H/i.pc detections and ALFALFA where available in Figure 6. Overall, the H /i.pc masses from FAST detections are in good agreement with those from the ALFALFA survey in orange dots, with median difference of ∼ 3 percent. The FAST-H /i.pc detections tend to be systematically more H /i.pc-massive than the MIGHTEE-H /i.pc detections by a median fraction of ∼ 7 percent for the isolated galaxies, i.e. those with a single MIGHTEE catalogued source within the FAST beam, denoted by blue dots. This tendency is noticeable from the intermediate to high H /i.pc mass end (i.e. ≥ 10 8 M ⊙ ), which is not unexpected as the FAST-H /i.pc data are likely to pick up more faint diffuse H /i.pc gas than the interferometric MIGHTEE-H /i.pc observations. The origin of this limited diffuse gas could be due to the gas accretion from the large scale structure through the circumgalactic medium (Hummels et al. 2019; Sardone et al. 2021), and/or feedback from the galactic winds and fountains given that the galaxies in blue dots  \nFigure 6. 𝑀 HI comparison of FAST-H /i.pc with MIGHTEE-H /i.pc (blue dots) and ALFALFA detections (orange stars). The x-axis is the MIGHTEE-H /i.pc or ALFALFA H/i.pc mass. The MIGHTEE-H /i.pc data with the robustness of 0.5 are adopted for the mass comparison. The diagonal grey line is the one-to-one relation. Note that the black circles are the detections that have at least one close H /i.pc or optical companion within the FAST beam, therefore their masses could be elevated due to multiple close companions.  \n<!-- image -->  \nare relatively isolated (Suresh et al. 2015; Tumlinson et al. 2017). We present the FAST-H /i.pc catalogue for isolated galaxies in Table 3.  \nThe black circles in Figure 6 are the H /i.pc detections that have at least one additional close H /i.pc or optical companion within the FAST beam, therefore their masses could be elevated as their flux measurements cannot be well separated from their companions on the FAST-H /i.pc images with a beam size of ∼ 3 ' . For example, the larger oranges dots (case 1) in Figure 7 are the cases where there are two MIGHTEE-H/i.pc counterparts within a FAST beam, and the blue line connected dots (case 2) are the FAST-H /i.pc galaxy pairs bridged by a thin H /i.pc structure as shown for one of them in Figure 5. The green circles (case 3) symbolise at least two DEVILS spectroscopic counterparts within the FAST beam and relative frequency range corresponding to a velocity width of ∼ 200 kms -1 , and the symbol size correlates with the number of the optical counterparts. In cases 1 and 2, we add up all the MIGHTEE-H /i.pc masses within the FAST beam or all the H /i.pc masses in a galaxy pair to study the total H /i.pc mass difference. We exclude the FAST-H /i.pc sources covered by the green circles (case 3) in Figure 7 from the remaining cases to investigate H/i.pc mass difference between FAST-H /i.pc and MIGHTEE-H /i.pc in the next section. However, they are included for studying the FAST-H /i.pc massstellar mass relation where the DEVILS stellar masses within the FAST beam can be combined together.", '3.2.2 𝑀 HI -𝑀 ★ relation': 'The H/i.pc and stellar mass ( 𝑀 HI -𝑀 ★ ) relation is often used to investigate the processes of gas consumption and star formation in galaxies. In particular, exploring this relation indicates that there is an upper limit for H /i.pc mass as a function of the stellar mass at high masses for H /i.pc-selected samples (Huang et al. 2012; Maddox et al. 2015; Parkash et al. 2018; Pan et al. 2023). We show the 𝑀 HI -𝑀 ★ relation from our FAST-H /i.pc sample in the upper panel of Figure 8. We bin the H /i.pc sources in stellar mass with a bin width of 0.6 dex to  \nTable 3. FAST-H/i.pc COSMOS catalogue of isolated galaxies. Column 1: Source name; Column 2; FAST-H /i.pc redshift with uncertainties less than 0.0001 estimated by the flux-weighted centroid position of source; Column 3: MIGHTEE-H/i.pc redshift with uncertainties of ∼ 0.0005; Column 4: Optical spectroscopic redshift with the index indicating the reference; Column 5: FAST-H/i.pc mass with 1 𝜎 uncertainties; Column 6: Rest-frame velocity width 𝑊 50 in the unit of kms -1 . Note that the indexes of a, b, c, d, e, f, g and h for the optical redshifts are from Polzin et al. (2021), Ann et al. (2015), Ahumada et al. (2020), Davies et al. (2018), Prescott et al. (2006), DESI Collaboration et al. (2024), Lilly et al. (2007) and Sohn et al. (2019).  \nestimate the mean H /i.pc mass (blue dashed line) for isolated galaxies (small blue dots) and galaxy pairs (large blue dots). We find that the 𝑀 HI -𝑀 ★ relation from FAST-H /i.pc is in good agreement with the relation from the ALFALFA-SDSS galaxy sample in Maddox et al. (2015) above the stellar mass of ∼ 10 8 M ⊙ , and detect a slope transition at 𝑀 ★ ∼ 10 9 . 3 M ⊙ . Although our sample size of 35 is small after excluding a few case 3 sources, the FAST-H /i.pc shows a higher H /i.pc mass by a mean (or median) value of ∼ 0.04 dex (ranging up to ∼ 0.32 dex) when compared to the 𝑀 HI -𝑀 ★ relation from the MIGHTEE-H /i.pc COSMOS sample, which is consistent with what we found in the one-to-one mass comparison in Figure 6 as the number of case 1 and 2 galaxies are low. The FAST-H /i.pc mass excess is most noticeable in the intermediate stellar masses which corresponds to intermediate and high H /i.pc masses due to the non-linear relationship between H /i.pc and stellar masses.  \nIn the lower panel of Figure 8, we show the H /i.pc mass difference between the FAST-H /i.pc and MIGHTEE-H /i.pc samples. The small dots are the isolated galaxies, whereas the large dots are case 1 and 2 galaxies which reside in a galaxy pair environment. Here we assume that the MIGHTEE-H /i.pc gas is within the galaxies, hence we consider that the H /i.pc mass difference between MIGHTEE and FAST for  \nFigure 7. Angular distribution of the cross-matched FAST-H /i.pc COSMOS sample with MIGHTEE-H /i.pc and DEVILS surveys. The blue dots are the crossmatched FAST-H /i.pc sources. The orange and green dots are the MIGHTEE-H /i.pc and DEVILS counterparts within the FAST beam, respectively. The symbol size correlates with the number of the counterparts within the FAST beam, and we only show the DEVILS sample with at least two counterparts. The blue line connected dots are the FAST-H /i.pc galaxy pairs bridged by a thin H /i.pc structure. The blue circles are the faint FAST-H /i.pc detections missed from our initial visual source finding but retrieved from the MIGHTEE-H /i.pc positions. The red and black boundaries enclose the main surveyed areas from FAST-H /i.pc and MIGHTEE-H/i.pc respectively.  \n<!-- image -->  \nisolated galaxies indicates the amount of the diffuse H /i.pc gas in the CGM and/or IGM, and the mass difference between the FAST-H /i.pc and MIGHTEE-H/i.pc for galaxy pairs is indicative of the H /i.pc gas in the CGMand IGM together.  \nIt appears that the total H /i.pc mass fraction in the IGM and CGM for the galaxy pairs is statistically higher than the H /i.pc fraction for the isolated galaxies by ∼ 13 percent with the former being ∼ 15 ± 4 percent and the latter being ∼ 2 ± 2 percent. The median fractions of H /i.pc gas surrounding the galaxy pairs and isolated galaxies are 18 percent and 7 percent respectively, and the mass excess is weakly dependent on the stellar mass with slightly larger excess at the intermediate stellar masses based on our limited size of the H /i.pc sample. This fraction difference is likely due to the contribution of the faint diffuse H /i.pc gas in the IGM where the pristine gas is distributed as a reservoir or from a recent tidal interaction event and will possibly be accreted onto galaxies at a later stage (Sancisi et al. 2008; Wolfe et al. 2013; Zhu et al. 2021; Zhou et al. 2023).  \nThegreendotsintheupperpanelofFigure8arelikely small galaxy groups including potential galaxy pairs, and their FAST-H /i.pc massstellar mass relation aligns with the average H /i.pc mass-stellar mass relation of the rest H /i.pc sample. The upcoming catalogue from source finding with the MIGHTEE 32k data should bring more information about the H /i.pc gas distribution in the member galaxies of these groups, and help to determine the H /i.pc gas fraction in the intragroup medium, which will shed light on the galaxy evolutionary stage.  \nFigure 8. Top: 𝑀 HI -𝑀 ★ relation of the cross-matched FAST-H /i.pc COSMOS sample with MIGHTEE-H /i.pc. The blue and orange dots are from FAST-H /i.pc and MIGHTEE-H/i.pc samples respectively, while the green dots are the FAST-H /i.pc galaxies with at least two DEVILS counterparts. The grey circles are the derived 𝑀 HI -𝑀 ★ relation from the ALFALFA galaxies by Maddox et al. (2015). Bottom: logarithmic H /i.pc mass difference between the FAST-H /i.pc and MIGHTEE-H/i.pc samples as a function of the stellar mass. The symbol size correlates with the number of the MIGHTEE-H /i.pc and DEVILS counterparts within a FAST beam or in a galaxy pair in both panels.  \n<!-- image -->', '3.2.3 Sources of uncertainty': "The approach of measuring the diffuse H /i.pc gas in CGM and IGM in this paper is to subtract the MIGHTEE-H /i.pc gas from the total H/i.pc content determined by the FAST-H /i.pc gas, which relies on an assumption that the MIGHTEE-H /i.pc measurement can well describe the H /i.pc gas within the galaxies. To examine this assumption, we investigate the impact of two sources of uncertainty: 1) definition of H/i.pc galaxy boundary, 2) missing flux due to the lack of zero spacing. We assume that the lack of additional field sources in the sky model for the primary calibrator used for the MIGHTEE-H /i.pc observations does not significantly bias the resulting H /i.pc mass measurements (see e.g. Heywood et al. 2020).  \nThe H /i.pc galaxy boundary is usually defined by exceeding a certain H/i.pc column density, ranging from 10 19 to ∼ 10 20 cm -2 . For example, Wang et al. (2016) and Rajohnson et al. (2022) use a column density of 1 . 25 × 10 20 cm -2 that corresponds to a surface mass density of 1 M ⊙ pc -2 to determine the H /i.pc disk diameter. However, there could be fainter or smaller H /i.pc structures below that threshold but still associated with the galaxies, such as extraplanar H /i.pc clouds and H /i.pc tails (Dedes & Kalberla 2010; Xu et al. 2021). Pingel et al. (2018) and Sardone et al. (2021) found that the cumulative H /i.pc mass stays mostly flat below a critical column density of a few times 10 19 cm -2 , and defined the diffuse neutral fraction as the fraction of H /i.pc below column densities of 10 19 cm -2 , which is of the same order as the sensitivity level of the MIGHTEE-H /i.pc 32k data with r = 0.5 over a  \n20 kms -1 velocity width. Thus we adopt a 3 𝜎 clipping approach to define the galaxy boundary, and compare the total flux within that boundary to the flux within a 1 ' -diameter aperture (up to 2 ' for a few exceptions) in Figure 9. The differences between the blue dots and orange stars are minimal for both the robust-0.0 and robust-0.5 data in the left and right panels, indicating that the 3 𝜎 clipping approach is capable of recovering most of the H /i.pc gas that is detectable by the MeerKAT.  \nIn Figure 10, we show an example of extracting the FAST-H /i.pc and MIGHTEE-H/i.pc fluxes in a spiral galaxy with a stellar mass of 1 . 7 × 10 9 M ⊙ . The FAST-H /i.pc mass is 8 . 41 ± 0 . 54 × 10 8 M ⊙ which has a H /i.pc mass excess of ∼ 27 ± 5 percent compared to the MIGHTEE-H /i.pc mass of ∼ 6 . 17 ± 0 . 23 × 10 8 M ⊙ . The MIGHTEE-H/i.pc moment-0 map with the r = 0.0 is shown in the left panel, while the MIGHTEE robust0.5 moment-0 and the FAST-H /i.pc maps are shown in the middle panels. The robust-0 image has the highest resolution of ∼ 12 arcseconds, and can therefore partially resolve this galaxy, featuring a faint H /i.pc tail in the south-east region of the galaxy. The FAST-H /i.pc excess as a function of frequency is clearly shown in the right panel where the MIGHTEE-H/i.pc flux densities are extracted from a circular aperture indicated by the white dashed circles in the moment-0 maps with a diameter of 1 arcminute.  \nThe impact of the zero spacing issue is strong for the robust-0.0 data across wide flux scales in the left panel of Figure 9, but only appears to arise at the high flux end in the right panel, where the FAST-H/i.pc fluxes are predominately larger than the MIGHTEE-H /i.pc fluxes even for the robust-0.5 data with the maximum recoverable angular scale of ∼ 1 arcminute. This is unsurprising as galaxies with intermediate and high H /i.pc masses dominate the high flux end while dwarf H /i.pc galaxies tend to dominate the low flux end and lack the ability of hosting a large fraction of H /i.pc gas in their circumgalatic medium, therefore there is no strong FAST-H /i.pc flux excess for the dwarfs. As some galaxies in our sample are larger than 1 arcminute, we note that there could be a small fraction of flux missing from the MIGHTEE-H/i.pc robust-0.5 data when the H /i.pc masses within the galaxies are determined. However, considering that MeerKAT's antennas are densely packed and the number of galaxies larger than 1 arcminute are not dominating our sample, the impact of the zero spacing issue on the flux measurements within the galaxy disks for the robust-0.5 data is limited.", '3.3 FAST-MeerKAT synergy': "With the ∼ 3 ' beam size of FAST, we can at least partially resolve some galaxy groups and pairs, or even a few large nearby galaxies in our data (e.g. Wang et al. 2023; Lin et al. 2023). By convolving the MIGHTEE-H/i.pc data inside these group and pair members or large nearby galaxies to the same angular resolution of FAST, we are able to subtract the individual galaxy contributions from the FAST intensity map, leaving the residual that is the diffuse H /i.pc gas distributed in IGM/CGM resolved out by MeerKAT. For dwarf galaxies at low redshifts and galaxy groups at high redshifts ( 𝑧 > 0 . 2), most of the targets cannot be resolved by FAST. We will therefore concentrate efforts on measuring the total amount of H /i.pc gas in IGM/CGM within various environments (e.g. isolated galaxies or galaxy groups) by subtracting the MIGHTEE-H /i.pc fluxes from the FAST-H /i.pc images.  \nWe are also able to cross-correlate the FAST H /i.pc intensity maps with optical galaxy catalogues (e.g. Chang et al. 2010; Tramonte & Ma 2020; Wolz et al. 2022) in order to mitigate foreground and systematic effects and thus help to constrain cosmological parameters such as the total H /i.pc mass density ( Ω HI ), the dark matter power spectrum and their evolution. Moreover, the cross-correlation with  \nFigure 10. Example of extracting the FAST-H /i.pc and MIGHTEE-H /i.pc fluxes in a nearby spiral galaxy. From left to right panels, the MIGHTEE robust-0.0, robust-0.5, FAST-H /i.pc moment-0 maps, and the flux density as a function of the frequency. The white dashed circles in the moment-0 maps indicate the 1 ' boundary and the white dashed lines are tracing the same angular size on different scales between the two middle panels. The blue dots in the corners are the corresponding beam sizes. The lowest column density contour on the FAST-H /i.pc moment-0 map is about 6 × 10 17 cm -2 which is ∼ 3 times of the local noise level, and the contours increase by a factor of 2. The background optical image in g-band is from DECaLS.  \n<!-- image -->  \nFigure 9. H/i.pc flux comparison between FAST and MeerKAT for isolated galaxies. The MIGHTEE-H /i.pc robust-0.0 and 0.5 data against the FAST-H /i.pc data are on the left and right panels respectively. For both panels, the blue dots and orange stars are the measurements with a 3 𝜎 clipping approach and an extraction within a 1-arcminute aperture respectively. The diagonal grey lines are the one-to-one relations.  \n<!-- image -->  \n<!-- image -->  \n1e 3  \n1390.0  \n1390.5  \n1391.0  \n1391.5  \nFrequency [MHz]  \nr=0.5  \n11  \nRA  \nFAST  \n12  \ns  \nRA  \nthe MIGHTEE-H/i.pc intensity map from interferometric observations can give a high-SNR detection of the large scale H /i.pc structure due to the exceptional sensitivity of both radio instruments, and provide a complementary view of the H /i.pc mass function from the H /i.pc detections alone.", '4 CONCLUSIONS': "The FAST-H/i.pc COSMOS survey presented in this paper demonstrates the capability of the FAST telescope in detecting the faint H /i.pc signals. We describe the data processing of the FAST raw data across the frequency band of 1310-1420 MHz, showing that the FAST data have reached a median 3 𝜎 column density of 𝑁 HI ∼ 2 × 10 17 cm -2 over a 5 kms -1 channel width. This sensitivity allows us to study the distribution of H /i.pc gas within galaxies, CGM and IGM when combined with the high-resolution of the MeerKAT telescope.  \nWe searched the H /i.pc spectral cube and visually identified a total of 80 sources, of which ∼ 70 percent are cross-matched with the  \nMIGHTEE-H/i.pc catalogue. We examined the H /i.pc mass against redshift and the H /i.pc flux against the velocity width 𝑊 50 , and we found that the non cross-matched sample is dominated by the low-flux or narrowline width H /i.pc galaxies. We compared the H /i.pc masses between the FAST-H/i.pc and MIGHTEE-H /i.pc 32k data and find that the FAST-H /i.pc masses are systematically higher than the MIGHTEE-H /i.pc masses by a median fraction of ∼ 7% for isolated galaxies alone.  \nWe further studied the 𝑀 HI -𝑀 ★ relation in the last billion years with the cross-matched sources, and find that the FAST-H /i.pc sources exhibit higher H /i.pc masses on average by a mean value of ∼ 0.04 dex (ranging up to 0.32 dex) for the isolated galaxies and galaxy pairs when compared to the H /i.pc masses from the MIGHTEE-H /i.pc COSMOS sample. We find that the major FAST-H /i.pc mass excess is contributed by galaxies with intermediate stellar masses based on our relatively small sample size.  \nBy separating the galaxy pairs from the isolated galaxies, we find the total H /i.pc gas fraction in the IGM and CGM together in the galaxy pairs is statistically higher than the gas fraction in the CGM (and/or  \ns  \n10  \nh  \n02  \nm  \n14  \ns  \n13  \ns  \n12  \ns  \n10  \ns  \n09  \ns  \n2°06'  \n03'  \n00'  \n1°57'  \n10  \nh  \n02  \nm  \n36  \ns  \n24  \ns  \n00  \ns  \n01  \nm  \n48  \ns  \nFlux density [Jy]  \n4  \n3  \n2  \n1  \n0  \nr=0.0  \nr=0.5  \nFAST  \nIGM)of isolated galaxies by ∼ 13%. The former being ∼ 15 ± 4%and the latter being ∼ 2 ± 2% suggest that the CGM and IGM associated with interacting galactic systems are H /i.pc gas richer than those surrounding the isolated galaxies (possibly due to the gas accretion from the nearby reservoir or tidal interaction between the pair galaxies), albeit with large uncertainties in the measurements as a result of the small sample size of the galaxies. Taken together, the diffuse H /i.pc in the CGM and the IGM accounts for ∼ 5 ± 2% of the total H /i.pc gas in relatively simple galactic environments (i.e. without considering the complex galaxy groups).  \nBy combining the capabilities of FAST and MeerKAT telescopes, the full potential of both instruments can be explored. We highlight the prospects of detecting the faint diffuse H /i.pc gas in IGM and CGM, and studying the evolution of the Ω HI . This FAST-H /i.pc COSMOS survey acts as a pilot project to form a foundation for further collaboration between MeerKAT and FAST on various astrophysical and cosmological applications with much wider area surveys.", 'DATA AVAILABILITY': 'The FAST-H/i.pc COSMOS spectral data are available upon reasonable request. All the raw data of the FAST observations will be freely accessible to any user after a twelve-month privileged period as per the FAST data release policy 2 .', 'ACKNOWLEDGEMENTS': 'We thank the anonymous referee for their suggestions that have improved this paper. We thank Di Li, Bradley Frank, Chun Sun, Tobias Westmeier and Hongwei Xi for useful discussions. This work has used the data from the Five-hundred-meter Aperture Spherical radio Telescope (FAST). FAST is a Chinese national mega-science facility, operated by the National Astronomical Observatories of Chinese Academy of Sciences (NAOC).  \nThe MeerKAT telescope is operated by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Innovation. We acknowledge use of the IDIA data intensive research cloud for data processing. The IDIA is a South African university partnership involving the University of Cape Town, the University of Pretoria and the University of the Western Cape. The authors acknowledge the Centre for High Performance Computing (CHPC), South Africa, for providing computational resources to this research project.  \nWe acknowledge the use of the Ilifu cloud computing facility www.ilifu.ac.za , a partnership between the University of Cape Town, the University of the Western Cape, the University of Stellenbosch, Sol Plaatje University, the Cape Peninsula University of Technology and the South African Radio Astronomy Observatory. The Ilifu facility is supported by contributions from IDIA and the Computational Biology division at UCT and the Data Intensive Research Initiative of South Africa (DIRISA).  \nHP, MJJ and IH acknowledge support from a UKRI Frontiers Research Grant [EP/X026639/1], which was selected by the European Research Council. HP, MJJ, MGS and IH acknowledge support from the South African Radio Astronomy Observatory which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Innovation. MJJ acknowledge generous support from the Hintze Family Charitable Foundation through the Oxford  \nHintze Centre for Astrophysical Surveys and the UK Science and Technology Facilities Council [ST/S000488/1]. HP also acknowledges the warm hospitality at the Purple Mountain Observatory in 2021. YZM is supported by the National Research Foundation of South Africa under Grants No. 150580, No. ERC23040389081, No. CHN22111069370 and No. RA211125652169. MGS also acknowledges support from the National Research Foundation (Grant No. 84156). AAP, MJJ and IH acknowledge support of the STFC consolidated grant [ST/S000488/1] and [ST/W000903/1]. YC thanks the Center for Astronomical Mega-Science, Chinese Academy of Sciences, for the FAST distinguished young researcher fellowship (19-FAST-02). YC also acknowledges the support from the National Natural Science Foundation of China (NSFC) under grant No. 12050410259 and the Ministry of Science and Technology (MOST) of China grant no. QNJ2021061003L. YPJ is supported by NSFC (12133006) and by 111 project No. B20019.', 'APPENDIX A: INJECTION AND RECOVERY': 'Since the noise in our FAST-H /i.pc data cube does not follow a perfect Gaussian distribution (as shown in Figure 2) due to complex systematic effects present in our data, we use fake sources with known H/i.pc masses to verify that optimal input parameters are employed by SoFiA-2 on our FAST data.  \nWe first replace the visually-identified sources in the H /i.pc spectral cube with noise from their neighboring channels, and inject fake H /i.pc sources with similar H /i.pc masses and velocity widths of the real detections to random positions centered at 𝑧 HI with a deviation of 0.001. We then repeat the injection 10 times for each detection to avoid low-number statistics. For each iteration, we randomly generate the flux density profile, following a form of the busy function (Westmeier et al. 2013; Pan et al. 2020).  \nWeshow the recovered H /i.pc mass 𝑀 Re HI against the injected H /i.pc mass 𝑀 In HI in the top panel of Figure A1 with three flux thresholds, i.e. scfind.threshold =[4, 5, 6], used by the Smooth+Clip finder of SoFiA2. The median ratio of 𝑀 Re HI to 𝑀 In HI as a function of the injected H /i.pc mass is shown in the bottom panel. As our sample is dominated by the bright H /i.pc sources, the observed median mass ratios of [1.026, 1.01, 0.994] are very close to 1, in a decreasing order for the increasing threshold from the left to right panels as expected. We also model a Gaussian background noise in comparison with the real noise, and find that the recovered mass ratio follows a similar trend to the case of real noise but with smaller uncertainties due to the well-behaved fake noise.  \nTo further examinate the accuracy of the recovered fluxes and their associated uncertainties, we show a histogram of the flux measurement error ( 𝑒𝑟𝑟 = 𝑆 Re -𝑆 In ) divided by the flux measurement uncertainty ( 𝜎 𝑠 ) from SoFiA-2. Under ideal conditions, 𝑒𝑟𝑟 / 𝜎 𝑠 should have a Gaussian distribution with a centroid of 0 and a standard deviation of 1. A Gaussian fit to the histogram with scfind.threshold =5 yields a centroid of 0.3 ± 0.1 and a standard deviation of 1.01 ± 0.06, which is close to the expected values and indicates that the flux uncertainties reported by SoFiA-2 are accurate as demonstrated by Westmeier et al. (2021). The small offset agrees with the median ratio of 1.01 as we see in the middle panel of Figure A1. When the threshold is increased to 6, the offset decreases slightly, as the green histogram shows in Figure A2. There appears to be a positive flux bias of 2-3 percent from SoFiA-2, which is assumed to be related to a default parameter of scfind.replacement as reported by Westmeier et al. (2021). Considering that the returned uncertainties are far larger than this small positive flux bias and we would like to retain a certain sensitivity to the diffuse H /i.pc gas surrounding the galaxies, we chose the scfind.threshold =5 as the best parameter to measure the total fluxes. We note that a choice of slightly higher threshold will not change our main results.', 'APPENDIX B: NON-DETECTIONS': 'Wecross-matched the FAST-H /i.pc and MIGHTEE-H /i.pc detections in the COSMOS field, and found that there are 9 sources that are present in MIGHTEE-H/i.pc catalogue but not visually identified in our FAST data. To examine the cause, we extract the fluxes from the FAST-H /i.pc cube at the location of MIGHTEE-H /i.pc detections, and compare their flux density profiles with those from MeerKAT data in Figure B1 and their H /i.pc masses in Figure B2, indicating that the FAST-H /i.pc fluxes are largely in good agreement with the MeerKAT fluxes. The fact that these sources did not pass our visual filtering is due to their low signal to noise ratios and/or narrower velocity widths, and we were  \nScfind.threshold=4  \n<!-- image -->  \n<!-- image -->  \nFigure A1. Top: H /i.pc mass comparison between injected and recovered sources. The H /i.pc masses are estimated by employing the SoFiA-2 with the flux threshold of scfind.threshold =[4, 5, 6] from left to right panels, respectively. The fake sources are injected to H /i.pc cubes with Gaussian noise in blue and Real noise in orange where the genuine sources are removed. The diagonal grey dashed lines are the one-to-one relations. Bottom: median ratio of the recovered H /i.pc mass 𝑀 Re HI to the injected H /i.pc mass 𝑀 In HI with 1 𝜎 uncertainties as a function of the injected H /i.pc mass.  \n<!-- image -->  \nFigure A2. Histogram of flux measurement error ( 𝑒𝑟𝑟 ) divided by the flux measurement uncertainty ( 𝜎 𝑠 ) derived by SoFiA-2 for the mock galaxies. The dashed, orange line shows the result of a Gaussian fit to the histogram with a peak position of 0.3 ± 0.1 and a standard deviation of 1.01 ± 0.06, when we chose the best threshold scfind.threshold =5 for the FAST-H /i.pc flux measurements.  \n<!-- image -->  \ndoing visual source finding on a coarser H /i.pc cube with a velocity resolution of 26 kms -1 . Although these faint sources show up when we zoom in the data cube, we decide to list them in Table B1 as nondetections, where the FAST-H /i.pc masses are estimated by integrating the flux densities within the vertical grey dashed lines in Figure B1.  \nTable B1. FAST-H/i.pc COSMOS catalogue of non-detections. The FAST-H /i.pc redshifts and masses are displayed in the second and third columns with the signal to noise ratio (SNR) in the fourth column.  \n<!-- image -->  \n1e 3  \nJ100016.3+022528  \nFAST  \nMeerKAT  \nFlux density [Jy]  \nFlux density [Jy]  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n1310  \n1311  \n1312  \n1e 3  \n1316  \n1317  \n1318  \nFrequency [MHz]  \n1e 3  \n1312  \n1314  \n1316  \n1e 3  \n1317.6  \n1318.0  \n1318.4  \n1318.8  \nFrequency [MHz]  \n1320  \n1322  \nFrequency [MHz]  \nFigure B1. Comparison of the H /i.pc flux density profiles between FAST and MeerKAT for 9 galaxies that are not visually identified in the FAST data. The grey dashed lines indicate the spectral limits used for calculating the integrated fluxes.Figure B2. H/i.pc mass comparison between FAST and MeerKAT for 9 galaxies that are not visually identified in the FAST data. The diagonal grey dashed line is the one-to-one relation.  \n<!-- image -->  \nFrequency [MHz]  \nJ100011.7+020441  \nFAST  \nMeerKAT  \nFrequency [MHz]  \nJ100031.3+020752  \nFAST  \nMeerKAT  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n1e 3  \nJ100023.8+022205  \nFAST  \nMeerKAT  \n1375  \n1376  \nFrequency [MHz]  \nJ100031.4+022553  \nFAST  \nMeerKAT  \n1e 3  \n1e 3  \nJ100206.3+022053  \nFAST  \nMeerKAT  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n1315.2  \n1315.6  \n1316.0  \n1316.4  \n1e 3  \nFrequency [MHz]  \nJ100115.5+022858  \nFAST  \nMeerKAT  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \n0.5  \n1.0  \n1337  \n1338  \nFrequency [MHz]  \nJ095957.4+015914  \nFAST  \nMeerKAT'}