## The most distant H i galaxies discovered by the 500 m dish FAST

Hongwei Xi , 1 Bo Peng , 1 Lister Staveley-Smith , 2, 3 Bi-Qing For , 2, 3 Bin Liu , 1 Ru-Rong Chen, 1 Lei Yu, 1,4 Dejian Ding, 1,4 Wei-Jian Guo, 1 Hu Zou, 1 Suijian Xue, 1 Jing Wang, 5,6 Thomas G. Brink, 7 WeiKang Zheng, 7 Alexei V. Filippenko, 7 Yi Yang, 7,8 Jianyan Wei, 1,6 Y. Sophia Dai, 9 Zi-Jian Li, 1,4,9 Zizhao He, 10,11 Chengzi Jiang, 11,12 Alexei Moiseev, 13 and Sergey Kotov 13

1 National Astronomical Observatories (NAOC), Chinese Academy of Sciences

20 Datun Rd.,Chaoyang District, Bejing, 100101, Bejing, China

2 International Centre for Radio Astronomy Research (ICRAR), University of Western Australia

35 Stirling Hwy, Perth, 6009, WA, Australia

3 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia

4 School of Astronomy and Space Science, University of Chinese Academy of Sciences

Bejing, 100049, Beijing, China

5 Guangxi Key Laboratory for Relativistic Astrophysics, School of Physical Science and Technology, Guangxi University

Nanning, 530004, Guangxi, China

6 Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences

Beijing, 100101, Beijing, China

7 Department of Astronomy, University of California

Berkeley, CA 94720-3411, USA

8 Physics Department and Tsinghua Center for Astrophysics (THCA), Tsinghua University

Beijing, 100084, China

9 Chinese Academy of Sciences South America Center for Astronomy (CASSACA), National Astronomical Observatories of China, Chinese Academy of Sciences

20A Datun Road, Beijing, 100101, Bejing, China

10 Purple Mountain Observatory, Chinese Academy of Sciences

No. 10 Yuanhua Road, Qixia District, Nanjing, 210023, Jiangsu, China

11 School of Astronomy and Space Science, University of Science and Technology of China

No. 96 Jinzhai Road, Baohe District, Hefei, 230026, Anhui, China

12 CAS Key Laboratory of Planetary Sciences, Purple Mountain Observatory, Chinese Academy of Sciences

No. 10 Yuanhua Road, Qixia District, Nanjing, 210023, Jiangsu, China

13 Special Astrophysical Observatory, Russian Academy of Sciences

Nizhnii Arkhyz, Zelenchukskiy region, Karachai-Cherkessian Republic, 357147, Russia

## ABSTRACT

Neutral hydrogen (H ) is i the primary component of the cool interstellar medium (ISM) and is the reservoir of fuel for star formation. Owing to the sensitivity of existing radio telescopes, our understanding of the evolution of the ISM in galaxies remains limited, as it is based on only a few hundred galaxies detected in H i beyond the local Universe. With the high sensitivity of the Fivehundred-meter Aperture Spherical radio Telescope (FAST), we carried out a blind H i search, the FAST Ultra-Deep Survey (FUDS), which extends to redshifts up to 0.42 and a sensitivity of 50 µ Jy · beam -1 . Here, we report the first discovery of six galaxies in H i at z &gt; 0 38. . For these galaxies, the FAST angular resolution of ∼ 4 ′ corresponds to a mean linear size of ∼ 1 3 . h -1 70 Mpc. These galaxies are among the most distant H i emission detections known, with one having the most massive H i content (10 10 93 . ± 0 04 . h -2 70 M ). Using recent data from the DESI survey, and new observations with the Hale, ⊙ BTA, and Keck telescopes, optical counterparts are detected for all galaxies within the 3σ positional uncertainty (0 5 . h -1 70 Mpc) and 200 km s · -1 in recession velocity. Assuming that the dominant source of

Corresponding author: Bo Peng pb@nao.cas.cn

Corresponding author: Lister Staveley-Smith

H i is the identified optical counterpart, we find an evidence of evolution in the H i content of galaxies over the last 4.2 Gyr. Our new high-redshift H i galaxy sample provides the opportunity to better investigate the evolution of cool gas in galaxies. A larger sample size in the future will allow us to refine our knowledge of the formation and evolution of galaxies.

Keywords: H i line emission(690) - High-redshift galaxies(734) - Galaxy evolution(594)

## 1. INTRODUCTION

Hydrogen is the most common element in the Universe. In its neutral form (H i ), it is abundant in spiral galaxies such as our Milky Way. It also indirectly regulates star formation by feeding the molecular clouds that are responsible for initiating star formation (Sancisi et al. 2008). H i commonly extends far beyond the stellar disk of galaxies, making it an excellent probe for tracing interactions between galaxies, or between galaxies and the circumgalactic medium (CGM).

suited to the exploration of gas evolution over 30% of the age of the Universe. Here, we report six H i galaxy candidates directly detected in the 21 cm emission line in the pilot survey of FUDS (FUDS0) that have higher redshifts (0 38 . &lt; z &lt; 0 40) . than the highest redshift CHILES detections, and are unlensed so that their intrinsic properties are well understood. In this paper, we use the flat universe model with parameters of H 0 = 70 h 70 kms -1 Mpc -1 , Ω M = 0 3, and Ω . Λ = 0 7. .

One of the most significant questions in astronomy is how the Universe evolved into its present form. To answer this, it is necessary to study galaxies over a range of cosmic times and therefore redshifts. The Hubble Space Telescope (HST) reported an oldest galaxy in optical band at z = 11 1 (Oesch et al. . 2016), which is only ∼ 400 Myr after the Big Bang. Such detections allow measurements of the evolution of the molecular gas fraction of galaxies (Heintz et al. 2023), the evolution of supermassive black holes (Behroozi &amp; Silk 2018), and tests of the galaxy formation theory inside the standard ΛCDM cosmological framework (Behroozi &amp; Silk 2018). However, by comparison with the detection of stellar components of galaxies in optical or infrared (IR) bands, detection of cool H i is much harder owing to the weakness of the 21 cm emission line. Almost all H i detections (Barnes et al. 2001; Meyer et al. 2004; Giovanelli et al. 2005; Freudling et al. 2011; Hoppmann et al. 2015; Xi et al. 2021; Catinella &amp; Cortese 2015; Verheijen et al. 2007; Koribalski et al. 2020; Meyer et al. 2009; Fern´ndez et al. 2016; Jarvis et a al. 2016; Blyth et al. 2016; Baker et al. 2018; Li et al. 2018; Zhang et al. 2019) are confined to z &lt; 0 35, with only two detections be-. yond this, the most distant unlensed detection being at z = 0 376 from the VLA CHILES survey (Fern´ andez . et al. 2016), and the most distant lensed H i galaxy detection being at z = 1 3 (Chakraborty &amp; Roy 2023). .

The FAST Ultra-Deep Survey (FUDS, Xi et al. 2022) is one of the deepest H i surveys, having a sensitivity of ∼ 50 µ Jy beam · -1 , a frequency resolution of 22.9 kHz (velocity resolution of 6.76 km · s -1 at z ∼ 0 4), . a redshift limit of z = 0 42 (a lookback . time of 4.4 Gyr), and an angular resolution of ∼ 4 ′ (physical resolution of 1.3 h -1 70 Mpc) at the redshift limit. Hence, it is well-

## 2. DETECTIONS BY FAST

FAST is not only the largest single-dish telescope in the world, but hosts a 19-beam receiver (Dunning et al. 2017), which dramatically boosts its capability to make multiple simultaneous observations and increases its ability to survey and map large areas of sky. For the first FUDS field (FUDS0), a region centred at RA(J2000) = 08 17 h m 12 , s Dec(J2000) = +22 10 48 0 ◦ ′ . ′′ was observed for a total of 95 hr. This yielded a sensitivity of 50 µ Jy · beam -1 over an area of 0.72 deg 2 . The observations were performed in on-the-fly (OTF) mode from 2019 August 25 to 2020 May 22. More details of the observations, calibration, data reduction, and noise analysis are described in Appendix A and B.

We searched for sources with a signal-to-noise ratio (S/N) larger than 7 in the final cube by averaging channels on different scales (25 kHz, 50 kHz, 100 kHz, and 200 kHz) using a custom algorithm (see section 4.1 in Xi et al. 2022). Our initial study focuses on H i detections with z &gt; 0 38 ( . f &lt; 1 03 GHz). . Candidate signals were visually inspected to remove obvious artefacts e.g., some discrete ranges of frequency with bad spectral baselines contain an excess of positive and negative signals due to radio frequency interference. The final list consists of six H i galaxy candidates with z &gt; 0 38. . We present H i column-density contours overlaid on top of DESI Legacy Imaging Survey 1 (Zou et al. 2017; Dey et al. 2019) 3-band optical images and their spectra in Figure 1. A two-dimensional Gaussian fit to the columndensity maps is used to derive positions with a 1 σ uncertainty of 0 52 . ′ (based on the offset of 62 FUDS galaxies

1 https://www.legacysurvey.org/

10 20

Figure 1. Six high-redshift H i detections in the FUDS0 field. Upper panel: the contours of column density overlaid on optical images from the DESI Legacy Imaging Surveys. The cyan crosses indicate the positions of the optical counterparts, for which zoomed-in images are shown in the bottom-right insets. The red circles are the beam size in the final cube. Lower panel: the spatially integrated spectra from the final cube. The red lines are the best-fit Busy Function, and the black lines with double arrow heads indicate the line width W Cor 20 after correcting for frequency resolution. Optical spectra from the DESI legacy survey (ID: 123, 124, 125, and 126), Keck (ID: 127), and Hale (ID: 128) are overlaid in cyan (H α for galaxies 124, 125, 127, and 128; [O II] λ 3727 (refers to the sum of [O II] λ 3726 and [O II] λ 3729) for galaxies 123 and 126 owing to the absence of H ). α The second candidate optical counterpart of galaxy 123 is also shown, both image and spectrum ([O II] λ 3727 captured by BTA), in orange.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00419v1-figure-9fb7605b-3cb1-4dad-80a0-5ecf4d3ae7cc.jpg)

from their optical counterparts with spectroscopic redshift measurements). The Busy Function (Westmeier et al. 2014) was employed to derive model spectra, which are overlaid in Figure 1.

In order to further confirm the detections, we split our raw data into two (P1 and P2). The P1 data are from observations performed between 2019 August 25 and 2020 March 01; P2 data are from observations between 2020 March 02 and 2020 May 22. Both sets of data have similar integration times, yielding two independent cubes. Each cube has higher noise by a factor of ∼ √ 2. Using similar line widths and beam sizes, HI

column-density images and spatially integrated spectra were produced from both datasets . All six detections were identified in both P1 and P2 cubes. Cubes from two individual polarization channels also show the reality of the six detections.

The H i properties of the six galaxy candidates are listed in Table 1 (more details in Appendix Table 2). They lie in the redshift range of 0 389 . &lt; z &lt; 0 402 . and have logarithmic H i masses, log( M H i /h -2 70 M ), be-⊙ tween 10.37 and 10.93. The most H -massive galaxy i candidate (ID 124) has the largest H i mass to date. Considering the 1 σ uncertainty, it is slightly more

Table 1. The H i and optical properties of the six high-redshift galaxy candidates. Column 1 shows the candidate ID in the FUDS0 catalogue. The J2000 coordinates and redshifts from the H i emission line are given in columns 2 and 3, the rest-frame line width ( W Cor 20 ), corrected for instrumental resolution, is given in column 4 (the frequency resolution is 22.9 kHz, equivalent to 6.8 km s -1 at z ≈ 0 4). . The logarithmic H i mass is given in column 5. J2000 coordinates of the identified optical counterparts are listed in column 6, while the spectroscopic redshift is given in column 7. Logarithmic stellar mass ( M ∗ ) and SFR are given in columns 8 and 9, respectively. Column 10 lists the spatial offset in h -1 70 kpc for fitted H i positions from their identified optical counterparts. Column 11 gives their velocity separation in km · s -1 . The 1σ uncertainties are given in parentheses.

|          | H i properties                              | H i properties   | H i properties   | H i properties   | H i properties                       | Optical properties                          | Optical properties   | Optical properties                 | Optical properties                        |                          |                      |
|----------|---------------------------------------------|------------------|------------------|------------------|--------------------------------------|---------------------------------------------|----------------------|------------------------------------|-------------------------------------------|--------------------------|----------------------|
| ID - (1) | R.A., Dec (J2000) HH:MM:SS.S ± DD:MM:SS (2) | z - (3)          | W 20 km s - (4)  | Cor 1            | log( M H i ) log( h - 2 70 M ⊙ ) (5) | R.A., Dec (J2000) HH:MM:SS.S ± DD:MM:SS (6) | z - (7)              | log( M ∗ ) log( h - 2 70 M ⊙ ) (8) | log( SFR ) log( h - 2 70 M ⊙ yr - 1 ) (9) | offset h - 1 70 kpc (10) | ∆ cz km · s - 1 (11) |
| 123      | 08:17:47.1 +22:17:14                        | 0.38867          | (2)              | 587 (6)          | 10.54 (0.04)                         | 08:17:45.4 +22:17:20                        | 0.38860 (21)         | 11.28 (0.09)                       | 0.18 (0.45)                               | 127 (120)                | 21 (39)              |
| 124      | 08:18:06.6 +22:25:35                        | 0.39072          | (2)              | 372 (8)          | 10.93 (0.04)                         | 08:18:06.3 +22:25:19                        | 0.39070 (10)         | 10.97 (0.09)                       | 1.10 (0.14)                               | 87 (114)                 | 7 (19)               |
| 125      | 08:16:46.2 +22:11:49                        | 0.39563          | (2)              | 427 (18)         | 10.51 (0.04)                         | 08:16:48.1 +22:11:52                        | 0.39600 (10)         | 11.12 (0.08)                       | 0.88 (0.18)                               | 148 (125)                | 112 (31)             |
| 126      | 08:17:00.7 +22:04:58                        | 0.39612          | (1)              | 366 (2)          | 10.37 (0.04)                         | 08:17:01.3 +22:06:13                        | 0.39610 (10)         | 11.25 (0.09)                       | 0.03 (0.37)                               | 401 (156)                | 7 (19)               |
| 127      | 08:17:04.2 +22:22:11                        | 0.39819          | (5)              | 478 (16)         | 10.79 (0.05)                         | 08:17:06.4 +22:22:00                        | 0.39774 (3)          | 10.60 (0.13)                       | -0.20 (0.36)                              | 170 (129)                | 136 (19)             |
| 128      | 08:17:53.2 +22:09:44                        | 0.40199          | (4)              | 258 (17)         | 10.56 (0.05)                         | 08:17:55.0 +22:09:51                        | 0.40158 (10)         | 10.92 (0.09)                       | 1.10 (0.22)                               | 145 (126)                | 134 (33)             |

massive than Malin 1 (log( M H i /h -2 70 M ) = 10 82 ⊙ . ± 0 06, . Lelli et al. 2010; Huang et al. 2014), HIZOA J0836-43 (log( M H i /h -2 70 M ) ⊙ = 10 88, . KraanKorteweg et al. 2005; Cluver et al. 2010), or the two highest-mass HIGHz (Catinella &amp; Cortese 2015) galaxies (J160938.00+312958.5, J165940.12+344307.8, both with an H i mass of log( M H i /h -2 70 M ) = 10 89). ⊙ . A caveat is that the H i mass may have contributions from other galaxies at the same redshift within the FAST beam. This requires high angular resolution observations and prohibitively long integration times (e.g., a few hundred hours with the Karl G. Jansky Very Large Array (VLA) to detect them at S/N = 5).

## 3. DISCUSSION

The comoving volume in the FUDS0 survey at 0 38 . &lt; z &lt; 0 42 is 5 81 . . × 10 4 h -3 70 Mpc . We computed the com3 pleteness, cosmic variance, and detectable volume for the six galaxies by following the methods used for the Arecibo Ultra-Deep Survey (AUDS, Hoppmann et al. 2015; Xi et al. 2021). The Σ 1 V max method (Schmidt 1968) gives us a higher H i mass function (H MF) (log(Φ) = i -2 85 . ± 0 43 . at log( M H i h 2 70 M -1 ⊙ ) = 10 5) . than that from local surveys, HIPASS and ALFALFA. This indicates a more extended H MF at the high-mass end, i more massive H i galaxies than in the local Universe, and significant evolution over the last ∼ 4 2 Gyr. . The full AUDS catalogue (Xi et al. 2021) previously suggested a weak evolutionary trend for an increase of characteristic mass of H MF with increasing redshift. i Our result is consistent with this trend, and may indicate rapid gas consumption following the epoch of peak star-formation density at z ≈ 2. Again, our result is an upper limit on H MF taking into account the confusing mass. i

In our sample, all six optical counterparts are found within a spatial separation of 0.5 h -1 70 Mpc (the 3σ uncertainty of fitted H i positions) and a velocity separation &lt; 200 km · s -1 . Four candidates (ID: 123, 124, 125 and 126) have spectroscopic confirmations from the unreleased data of the Dark Energy Spectroscopic Instrument (DESI, DESI Collaboration et al. 2016a,b). The optical counterparts of the two remaining galaxy candidates (ID: 127 and 128) were identified in spectroscopic observations using the Keck 2 and Hale 3 telescopes, respectively. For galaxy 123, our Big Telescope Alt-azimuth (BTA) data reveals a second candidate optical counterpart at a similar redshift. Being further from the beam centre, this is not the preferred counterpart. Nevertheless, it may contribute confusing flux to the H i spectrum. The H α (or [OII], if both Lyman and Balmer series are not detected) emission lines are overlaid on top of H i spectra in Figure 1 (refer to Section D in the Appendix for the whole optical spectra). To derive the stellar mass ( M ∗ ) and star-formation rate ( SFR ), we employed MAGPHYS 4 (da Cunha et al. 2008, 2012) to fit spectral energy distributions (SEDs). The optical data cover five bands ( u , g , r , i , and z ) from SDSS, while the IR data cover four bands (W1, W2, W3, and W4) from the Wide-field Infrared Survey Explorer (WISE , Wright et al. 2010). 5 In Appendix D, we describe the fitting method in more detail. The optical properties are listed in Table 1 (more details given

2 https://keckobservatory.org/

3 https://sites.astro.caltech.edu/palomar/about/telescopes/hale.html

4 http://www.iap.fr/magphys/

5 https://wise2.ipac.caltech.edu/docs/release/allwise/

Figure 2. The relations between H i ( H i mass and gas fraction) and optical/IR (stellar mass, star-formation rate, stellar mass surface density, and concentration index) properties. Note that we use the H i mass as upper limits (red plus signs) for our six galaxy candidates. For comparison, the distribution of xGASS galaxies is also shown with the grey dots. The light-green triangles represent H i upper limits in xGASS. The galaxies from HIGHz and CHILES are shown by green squares and the orange dot, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00419v1-figure-12522f88-fb7f-4841-8215-3d8bca248a54.jpg)

in Appendix Table 2). The optical images and SEDs are also shown in Appendix Figure 6.

We now compare the H i properties with the derived stellar properties, namely M ∗ , SFR , stellar mass surface density ( µ ∗ = M / πR ∗ 2 2 50 ), and concentration index ( R 90 /R 50 ), comparing with galaxies from the extended GALEX Arecibo SDSS survey 6 (xGASS, Catinella et al. 2018), CHILES, and the HIGHz (Catinella &amp; Cortese 2015) sample of 39 H i galaxies with 0 17 . &lt; z &lt; 0 25. . This gives us an insight into possible evolution in the properties of the most H i -massive galaxies. Considering possible confusion, our results are upper limits for M H i and gas fraction ( f H i = M H i /M ∗ ). In the upper panels of Figure 2, M H i is shown as a function of M ∗ and SFR . The FUDS galaxy candidates clearly have large values of M H i and SFR , but consistent with the properties of HIGHz galaxies. Unsurprisingly, the FUDS candidates also have high M ∗ , though not quite as high as the highest in xGASS. Given the tightness of the M H i -SFR correlation, depletion times (the ratio of the two) are also consistent with xGASS and HIGHz galaxies. In the lower panels of Figure 2, the gas fraction is shown as a function of µ ∗ and R 90 /R 50 . The FUDS candidates show reasonable consistency with the xGASS scaling re- lations. The gas fractions for the FUDS candidates are high, though well below the most gas-rich galaxies in xGASS. Values for µ ∗ and R 90 /R 50 are low, but consistent with HIGHz galaxies and the trend for xGASS galaxies. Low values for µ ∗ , indicative of a large halo spin parameter, are also consistent with the ALFALFA HIghMass (Huang et al. 2014) sample. Therefore, the FUDS galaxies appear to be extreme high-H i mass analogues of galaxies in the local Universe, with properties similar to the lower-redshift, SFR -selected HIGHz sample. The results here indicate that such galaxies are more numerous at z ≈ 0 4 . than in the local Universe, therefore suggesting significant evolution over the last ∼ 4 2 Gyr. . Although there is only a single CHILES galaxy at a similar redshift ( z = 0 376), it has similar . properties to the FUDS candidates, albeit with a higher SFR , and is therefore consistent with this picture.

Furthermore, using the stacking technique at higher redshift ( z ≈ 1), optically selected blue galaxies appear to have H i masses that are a factor of 3.5 higher for their stellar mass than their z = 0 counterparts (Chowdhury et al. 2021, 2022). Depletion timescales have also been reported to be lower at this redshift (Chowdhury et al. 2021, 2022) and at intermediate redshifts ( z ≈ 0 35) . (Sinigaglia et al. 2022; Bera et al. 2023). Considering the possible confusion, the real depletion time in our sample could be lower, and more consistent with

CHILES. These results therefore also seem consistent with a galaxy evolution picture where gas consumption exceeds gas replenishment for high-mass galaxies following the peak of SFR density at 'cosmic noon', following which there is an almost exponential decrease of their H i mass with time and a gradual truncation of the highmass cutoff in the H i mass function from high redshift to low redshift. In the immediate future, the ongoing 'blind' FUDS survey will provide even more direct detections of individual galaxies in a larger survey volume, and therefore a more accurate measurement of the evolution of the H i MF, allowing further exploration of the cool-gas evolution of galaxies out to z ≈ 0 4. . Future observations with the Square Kilometer Array (SKA) will offer the possibility of both high sensitivity and high angular resolution, so will avoid the requirement to confirm optical counterparts through deep optical spectroscopy.

## 4. SUMMARY

The FAST Ultra-Deep Survey is a deep blind H i galaxy survey designed to explore the evolutionary trend in the cool gas content of galaxies. We have now completed the observation and data reduction for the first of the six target fields, FUDS0. In this paper, we focus on the highest-redshift detections detected in the field. The main findings are summarized below:

- · Six galaxies are discovered through blind detection in 21-cm emission line observations with redshifts z &gt; 0 38. . These are at the redshift limit of the 19-beam FAST receiver and first such direct detections at this redshift to date, other than through gravitational lensing.
- · All have optical counterparts identified using data from DESI or from new spectroscopic observations, which we conducted with the Hale, BTA and Keck telescopes.
- · Comparison with galaxies in xGASS reveals that the six galaxies are extreme high-H i mass analogues of galaxies in the local Universe.
- · The significantly larger HIMF that is implied for high-mass galaxies at z ≈ 0 4 indicates their sig-. nificant consumption of cool gas over the last ∼ 4 2 Gyr. .

## 5. ACKNOWLEDGMENTS

This work made use of the data from FAST (Fivehundred-meter Aperture Spherical radio Telescope). FAST is a Chinese national mega-science facility, operated by the National Astronomical Observatories, Chinese Academy of Sciences. The work is supported by the National Key R&amp;D Program of China under grant number 2018YFA0404703, and the FAST Collaboration. Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. This research used data obtained with the Dark Energy Spectroscopic Instrument (DESI). DESI construction and operations are managed by the Lawrence Berkeley National Laboratory. This material is based upon work supported by the U.S. Department of Energy (DOE), Office of Science, Office of High-Energy Physics, under Contract No. DE-AC02-05CH11231, and by the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility under the same contract. Additional support for DESI was provided by the U.S. National Science Foundation (NSF), Division of Astronomical Sciences under contract AST-0950945 to the NSF's National Optical-Infrared Astronomy Research Laboratory; the Science and Technology Facilities Council of the United Kingdom; the Gordon and Betty Moore Foundation; the Heising-Simons Foundation; the French Alternative Energies and Atomic Energy Commission (CEA); the National Council of Science and Technology of Mexico (CONACYT); the Ministry of Science and Innovation of Spain (MICINN); and the DESI Member Institutions 7 . Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the NSF, the DOE, or any of the listed funding agencies. Hu Zou and Suijian Xue acknowledge support from the National Natural Science Foundation of China with grant 12120101003 and the National Key R&amp;D Program of China with grant 2022YFA1602902. J. Wang is supported by the National Natural Science Foundation of China under grants 12173009 and the Natural Science Foundation of Guangxi (2020GXNSFDA238018). Y. Yang appreciates the generous financial support provided to the supernova group at U.C. Berkeley (PI: A.V. Filippenko) by Gary and Cynthia Bengier, Clark and Sharon Winslow, Sanford Robertson, and numerous other donors. A.V. Filippenko's group has received additional financial assistance from the Christopher R. Redlich Fund, Alan Eustace (W.K. Zheng is a Eustace Specialist in Astronomy), Briggs and Kathleen Wood (T.G. Brink is a Wood Specialist in Astronomy), and many other donors. The data presented herein were obtained in part at the W. M. Keck Observatory, which is operated as a scientific partnership among the Califor-

7 https://www.desi.lbl.gov/collaborating-institutions

nia Institute of Technology, the University of California, and the National Aeronautics and Space Administration (NASA); the observatory was made possible by the generous financial support of the W. M. Keck Foundation. This research also uses data obtained through the Telescope Access Program (TAP). Z. He would like to acknowledge support received from the CAS Project for Young Scientists in Basic Research (No. YSBR-063). We obtained part of the observed data on the unique scientific facility 'Big Telescope Alt-azimuthal' of SAO RAS, as well as reduced and analysed the galaxy spectra with the financial support of grant 075-15-2022-262

(13.MNPMU.21.0003) of the Ministry of Science and Higher Education of the Russian Federation. The authors would like to thank Guo Chen for help with the optical observations by the Palomar P200 telescope, and Dr. Marat Musin for timely coordinating the Director's Discretionary Time observation using the 6 m BTA telescope at SAO.

Facilities:

FAST, Hale, BTA, Keck

Software: astropy (Astropy Collaboration et al. 2013, 2018), LPipe (Perley 2019), Pypeit (Prochaska et al. 2020a,b), slinefit (Schreiber et al. 2018)

## APPENDIX

## A. OBSERVATIONS

The FUDS survey consists of six target fields aiming at exploring the gas evolution in galaxies. The FUDS0 field covers an area of 0.72 deg 2 , and was chosen to overlap with the GAL2577 field of the Arecibo Ultra-Deep Survey (AUDS, Freudling et al. 2011; Hoppmann et al. 2015; Xi et al. 2021). AUDS has a lower sensitivity ( ∼ 75 µ Jybeam -1 ) and a smaller redshift coverage ( z &lt; 0 16) than FUDS, but nevertheless provides a useful comparison with the FAST . results in order to verify the calibration and data-reduction methods. An important selection criterion is the avoidance of strong continuum sources. In the NVSS catalogue (Condon et al. 1998), there are four continuum sources with flux density greater than 50 mJy in or close to the FUDS0 field. The strongest source is located at RA(J2000) = 08 h 17 m 35 , s Dec(J2000) = +22 37 11 9 with a 20 cm flux density of 1.28 Jy. ◦ ′ . ′′

The observations were taken between 2019 Aug 25 and 2020 May 22, mostly during the FAST commissioning phase. The FUDS0 field was scanned in both the RA and Dec directions with the 19-beam receiver in on-the-fly (OTF) mode (Xi et al. 2022). The 500MHz bandwidth (1.0-1.5 GHz) was split into 65,536 channels, which gives a frequency resolution of 7.63 kHz (with equivalent velocity resolution 2.25 km s -1 at z = 0 4). . The spectra were recorded with a 1 s integration time. The total observation time was 129 hr, consisting of 15 hr on the calibrator, 95 hr on FUDS0, and 19 hr overhead.

## B. DATA REDUCTION

The data-reduction method has been previously described (Xi et al. 2022). The final cube for the FUDS0 field covers 1 deg × 1 deg, as shown in the left panel of Figure 3, with a pixel size of 1 ′ × 1 . ′ The gridded beam size is 4 11 . ′ at z = 0 4 and its associated frequency dependence is given in our previous study (Xi et al. 2022). The Hanning smoothed . cube has a frequency resolution of 22.9 kHz (6.76 km s -1 at z = 0 4). . Since the field is not uniformly sampled, the noise rises from the centre to the edge, with the lowest noise being about 50 µ Jybeam -1 in the RFI-free frequency range in the centre of the field.

We used the flagging procedures developed in our previous study (Xi et al. 2022). Most flagging is associated with external radio frequency interference (RFI), such as radar ( ∼ 1 09 GHz), GNSS (1.15-1.30 GHz), and the geostationary . AsiaStar satellite ( ∼ 1 48 GHz). FAST also suffered from strong RFI generated by the compressor for the refrigerating . dewar before Aug 2021. Our flagging method was carefully designed to be robust against this RFI. The total flagged fraction is about 24%, including the H i emission line from the Milky Way.

In order to examine the noise distribution in the final data cube, we followed the procedures described previously (Hoppmann et al. 2015). For each voxel, we calculated the root-mean square (RMS) over the channels within ± 150 kms -1 , excluding any channels containing detections. The noise cube was then smoothed with a Gaussian function ( σ = 100kms -1 ) within ± 200 kms -1 . The resulting spatial noise distribution at 1.02 GHz is shown in the left panel of Figure 3. The lowest noise is distributed around the field centre. There is excess residual noise from a strong continuum source ( S ν = 1 28 Jy) at the upper edge. This impacts the detectability of nearby sources. . The right panel shows the central RMS noise as a function of frequency between 1.0 and 1.029 GHz ( z = 0 38). . Internal RFI has a

Figure 3. Left panel: the spatial RMS distribution at 1.02 GHz. The red plus sign indicates the position of the strongest continuum source in the field. The dotted red circle shows the extent of pixels strongly influenced by the continuum source. The red stars indicate the locations of the six galaxy candidates. Right panel: The RMS as a function of frequency for the central pixel of the final cube. The red bands show the frequency ranges strongly impacted by RFI. The frequencies of the six candidates are given by red dotted lines.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00419v1-figure-a14972d5-c00d-431e-a547-07c0f1fc507b.jpg)

modest impact on the RMS noise level, but its presence causes poor baselines and difficulty in identifying faint sources. Frequencies where source identification was strongly affected are shown in red in the panel.

## C. DETECTIONS

The algorithm used for the source finding is similar to that of our RFI finder (see Section 4.1 in Xi et al. 2022). The finder provides masks for detected signals and flattens the spectral baselines using polynomial functions. Detections with small spatial ( &lt; 2 pixels in either RA or Dec) or frequency ( &lt; 5 channels) extent are removed from our candidate list. Criteria for true detections are (1) peak flux density | S peak | ≥ 7 σ , where σ is the local RMS in the smoothed spectrum; (2) | S peak | ≥ 5 σ plus a distinct double-horn feature in the line profile, or with a beam-sized spatial extent; or (3) | S peak | ≥ 5 σ with matched redshifts from other surveys, including SDSS DR16 (Ahumada et al. 2020), AUDSOC (Hoppmann 2014), and AUDS (Xi et al. 2021). Negative detections that satisfy the above criteria were included for reliability analysis. The source finder did not separate galaxies joined in both spatial extent and frequency. These were later separated manually wherever possible based on their moment-zero maps.

Many false-negative detections were found around the strong continuum source or close to RFI frequencies. These detections were removed if they were within 5 ′ of the strongest continuum source (the dotted red line in the left panel of Figure 3) or within the red bands in the right panel of Figure 3, unless they satisfied at least one of these criteria: distinct double-horn features, beam-sized spatial extent, or spectroscopic redshift. Finally, we detected six sources at 0 38 . &lt; z &lt; 0 42. . No negative signals satisfying our selection criteria were detected in this redshift range; the reliability for the six galaxies is therefore high. Figure 3 shows their positions in space and frequency in the FUDS0 field. The contours of column density and spectra of the six galaxies are given in Figure 1 in the main text.

Since the physical beam size of FAST is about 1.32 h -1 70 Mpc at z = 0 4, it is reasonable to treat these high-redshift . galaxies as point sources. We used a two-dimensional (2D) Gaussian function with a fixed full width at half-maximum intensity (FWHM) corresponding to the model beam size in the final cube (Xi et al. 2022) to fit the moment-zero map (see Figure 4) to derive the position (R.A. and Dec). The spatial-integrated spectrum is derived from

$$S ( \nu ) = \frac { \sum _ { i j } I _ { i j } ( \nu ) w _ { i j } \Delta x \Delta y } { \sum _ { i j } w _ { i j } ^ { 2 } \Delta x \Delta y } \,,$$

where I ij ( ν ) is the intensity in pixel with indices of i , j , w ij is the amplitude of the 2D Gaussian fit, and ∆ x ∆ y is the solid angle of the pixel. The spatial-integrated spectra are shown in the lower panel of Figure 1 in the main text.

We used the Busy Function (Westmeier et al. 2014) to fit the spatially integrated spectrum to derive the integrated flux S int , flux density weighted frequency ν cen , peak flux density S peak ( ν ), and line width W 50 and W 20 . In order to

Figure 4. The moment-zero maps of the six galaxies from the final cube. The contours of column density are displayed by white lines. The red circles are the beam size in the final cube. And the red stars indicate the positions of optical counterparts.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00419v1-figure-08aee38f-c4f6-4830-aff1-04b82f562302.jpg)

compare with previous work, we correct for the instrumental effect on the line width using

$$W ^ { C o r } = \sqrt { W ^ { 2 } - f _ { i n s } ^ { 2 } } \,,$$

where W Cor is the corrected line width, W is the line width in units of Hz derived by fitting the Busy function, and f ins is the frequency resolution (22.9 kHz; equivalent velocity resolution 6.76 km s -1 at z = 0 4). .

Table 2. The H i properties of the six high-redshift galaxy candidates. Column 1 shows the galaxy ID in the FUDS0 catalogue. The coordinates in J2000 are given in column 2, while the redshift from H i emission line in column 3. We list the integrated flux ( S int ), peak flux density ( S peak ), and corrected line width ( W Cor 20 ) in columns 4-6, respectively. Column 7 provides the logarithm of the H i mass. The 1 σ uncertainties are given in parentheses.

|   ID DDD (1) | R.A. H i , Dec H i (J2000) HH:MM:SS.S, ± DD:MM:SS (2)   | z - (3)     | S int mJy MHz (4)   | S peak mJy (5)   | W Cor 20 MHz (6)   | log( M H i ) log( h - 2 70 M ⊙ ) (7)   |
|--------------|---------------------------------------------------------|-------------|---------------------|------------------|--------------------|----------------------------------------|
|          123 | 08:17:47.1, +22:17:14                                   | 0.38867 (2) | 0.156 (0.016)       | 0.098 (0.013)    | 2.003 (0.020)      | 10.54 (0.04)                           |
|          124 | 08:18:06.6, +22:25:35                                   | 0.39072 (2) | 0.381 (0.039)       | 0.444 (0.046)    | 1.269 (0.028)      | 10.93 (0.04)                           |
|          125 | 08:16:46.2, +22:11:49                                   | 0.39563 (2) | 0.143 (0.015)       | 0.132 (0.015)    | 1.450 (0.060)      | 10.51 (0.04)                           |
|          126 | 08:17:00.7, +22:04:58                                   | 0.39612 (1) | 0.103 (0.010)       | 0.118 (0.013)    | 1.241 (0.006)      | 10.37 (0.04)                           |
|          127 | 08:17:04.2, +22:22:11                                   | 0.39819 (5) | 0.266 (0.028)       | 0.234 (0.026)    | 1.618 (0.055)      | 10.79 (0.05)                           |
|          128 | 08:17:53.2, +22:09:44                                   | 0.40199 (4) | 0.154 (0.018)       | 0.266 (0.029)    | 0.871 (0.056)      | 10.56 (0.05)                           |

Considering that some galaxy candidates are detected close to RFI where the spectral baseline may not be flat, the Busy function can underestimate parameter uncertainties. Hence we employed the following jack-knife method to estimate their uncertainties by taking into account thermal noise and baseline noise around the galaxy spectral profiles:

- 1. The noise spectrum is extracted by subtracting the best-fit Busy function from the spatially integrated spectrum.
- 2. The noise spectrum is fractionally shifted by 1/30 of its length. An artificial spectrum is generated by summing the noise spectrum and the best-fit Busy function.
- 3. The artificial spectrum is refitted with a Busy function.

- 4. The above process is repeated 30 times to generate a set of 30 values for each parameter. The parameter uncertainty is estimated from the standard deviation of each set.

For flux uncertainties, we conservatively introduce a further uncertainty of 10% to account for other sources of error, including calibration, gridding, and frequency dependence. The relative errors of flux parameters (including S int and S peak ( ν )) are given by

$$\sigma _ { \text{rel} } ^ { \prime } & = \sqrt { \sigma _ { \text{rel} } ^ { 2 } + ( 1 0 \% ) ^ { 2 } } \,,$$

where σ ′ rel is the relative error after correction and σ rel is the relative error derived by the noise-shift method. The H i masses for the detections are derived with (Meyer et al. 2017):

$$M _ { H 1 } = 4 9. 7 \, D _ { \text{Lum} } ^ { 2 } \, S _ { \text{int} },$$

where M H i is the H i mass in h -2 70 M , ⊙ D Lum is the luminosity distance of the galaxy in h -1 70 Mpc, and S int is the integrated flux in Jy Hz. The derived properties of the H i galaxies are listed in Table 2.

## D. OPTICAL COUNTERPARTS

The Dark Energy Spectroscopic Instrument (DESI) is a multiplexed spectroscopic survey to explore dark energy by employing the Mayall 4 m telescope at Kitt Peak National Observatory and 5000 robotic fiber positioners (Levi et al. 2013; DESI Collaboration et al. 2022, 2016a,b). The survey observes galaxies and stars selected from the DESI Legacy Imaging Surveys (Cooper et al. 2023; Hahn et al. 2023; Zhou et al. 2023; Raichoor et al. 2023; Chaussidon et al. 2023; Zou et al. 2017; Dey et al. 2019; Moustakas et al. 2023; Myers et al. 2023). The target selection was validated during the DESI Survey Validation phase (Lan et al. 2023; Alexander et al. 2023; DESI Collaboration et al. 2023a,b). The DESI experiment necessitates multiple supporting software pipelines and products (Raichoor et al. 2023; Schlafly et al. 2023; Guy et al. 2023; Brodzeller et al. 2023). The data are processed and analysed by an extensive spectroscopic reduction pipeline and a template-fitting pipeline to derive classifications and redshifts for each targeted source (Guy et al. 2023; Brodzeller et al. 2023).

Figure 5. The optical spectra of the optical counterparts. The model continuum and lines are displayed in cyan. The detected emission/absorption lines are indicated by dotted lines. Note that the wavelength is in vacuum.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00419v1-figure-a44c81b3-02a2-47cd-964f-0a0dcc0b3309.jpg)

We found the optical counterparts of four H i galaxy candidates (ID: 123, 124, 125, and 126) in DESI survey data. The optical spectrum for galaxy candidate 127 was obtained with the Keck-I 10 m telescope and the Low Resolution Imaging Spectrometer (LRIS, Oke et al. 1995) on 2023 Oct 19. The 5600 ˚ dichroic, the 600/4000 blue A grism (dispersion 0.63 ˚ pixel A -1 ), and the 400/8500 red grating (dispersion 1.20 ˚ pixel A -1 ) provided a wavelength coverage of 3150-10,270 ˚ . The slit width was set to be 1 A ′′ during the observations and the exposure time was 900 s. The raw data were reduced by using the LPipe pipeline (Perley 2019), which performs a completely automated, endto-end reduction of LRIS spectra. A linear combination of a linear local continuum and two Gaussian profiles (H α and [N II] λ 6583 is adopted to model the emission-line profile in the H α region. Galaxy candidate 128 was observed on 2023 Oct 16 with the long-slit Double Spectrograph (DBSP) on Palomar P200 telescope. A total of 35 min on target exposure time yielded an optical spectrum, in which the H α , H , and [O III] β λ 5007 emission lines were detected at S/N = 41.0, 5.6, and 3.9, respectively. The spectroscopic data were reduced following the Pypeit pipeline (Prochaska et al. 2020a,b), and the redshift is derived using slinefit (Schreiber et al. 2018).

In our spectroscopic observations with the BTA, we discovered a second optical source (R.A.: 08:17:48.26, Dec: +22:17:50.7) with a similar redshift to galaxy 123 within the FAST beam. The spectrum was captured using SCORPRIO-2 (Spectral Camera with Optical Reducer for Photometric and Interferometric Observations, Afanasiev &amp; Moiseev 2011) at the 6 m BTA telescope of the Special Astrophysical Observatory (SAO) on 2023 Dec 20. The [O II] λλ 3726 3729 emission and H&amp;K Ca II absorption corresponding to a redshift of 0.3894 (5) was clearly identified. , Although it is at a larger separation ( ∼ 1 ), it may contribute additional flux to the H ′ i spectrum. All of the galaxies in our sample are confirmed with optical spectroscopic redshifts. Methods used in the SDSS survey 8 were employed to remove continuum and fit the emission/absorption lines. Figure 5 displays the whole optical spectra of the optical counterparts with detected lines marked. We did not detect any Lyman or Balmer lines in the spectra of the galaxies with IDs 123 and 126, which indirectly implies less ionized gas. Our measurements from radio observations also shows lower H i masses for these two galaxies.

Figure 6. The SEDs for the six galaxies. The red dots represent data from SDSS, while the orange 'x' represent data from WISE. The black lines show the best-fit SEDs by MAGPHYS, while the blue lines indicate the stellar spectra without attenuation by dust.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00419v1-figure-20bf8df2-a891-450a-b6c9-6defe1fca3ad.jpg)

We employed MAGPHYS 9 (da Cunha et al. 2008, 2012) to fit optical/IR photometry from SDSS ( u , g , r , i , and z ) and WISE ( W 1, W 2, W 3, and W 4) to derive the stellar mass ( M ∗ ), star-formation rate ( SFR ), and specific SFR ( sSFR ). MAGPHYS is designed for fitting the spectral energy distribution (SED) in the ultraviolet, optical, and near-IR bands. The code takes into account the stellar component, as well as the dust component which attenuates stellar-population spectra with emission at IR bands. In our data, we removed the bands with relative uncertainties larger than unity. The Milky Way extinction at each band was removed by using a dust-reddening map (Schlafly

8 https://classic.sdss.org/dr7/algorithms/speclinefits.php

9 http://www.iap.fr/magphys/

&amp; Finkbeiner 2011) and a dust-extinction function (Fitzpatrick 1999). Considering the possibility of contamination of the W 4 band of WISE by active galactic nuclei, we performed a refit without this band when the fit failed. This happened when best-fit parameters with an uncertainty of zero were derived, or when the boundary of parameter space was reached. The final fit results are given in Figure 6. The best-fit parameters of the optical counterparts are listed in Table 3.

Table 3. Properties of the optical counterparts of the six H i galaxy candidates. Column 1 gives the ID of the detections in FUDS0 survey. The coordinates of the SDSS galaxies are given in column 2. The spectroscopic redshift is given in column 3. The M ∗ , SFR and sSFR from fitting the SED by MAGPHYS are listed in columns 4-6, respectively. Column 7 indicates the status of the fit. The 1σ uncertainties are given in parentheses.

|   ID DDD (1) | R.A. Opt , Dec Opt (J2000) HH:MM:SS.S, ± DD:MM:SS (2)   | z - (3)      | log( M ∗ ) log( h - 2 70 M ⊙ ) (4)   | log( SFR ) log( h - 2 70 M ⊙ yr - 1 ) (5)   | log( sSFR ) log(yr - 1 ) (6)   | fit flag - (7)   |
|--------------|---------------------------------------------------------|--------------|--------------------------------------|---------------------------------------------|--------------------------------|------------------|
|          123 | 08:17:45.4, +22:17:20                                   | 0.38860 (21) | 11.28 (0.09)                         | 0.18 (0.45)                                 | -11.12 (0.48)                  | success          |
|          124 | 08:18:06.3, +22:25:19                                   | 0.39070 (10) | 10.97 (0.09)                         | 1.10 (0.14)                                 | -9.93 (0.23)                   | success          |
|          125 | 08:16:48.1 +22:11:52                                    | 0.39600 (10) | 11.12 (0.08)                         | 0.88 (0.18)                                 | -10.22 (0.21)                  | success          |
|          126 | 08:17:01.3 +22:06:13                                    | 0.39610 (10) | 11.25 (0.09)                         | 0.03 (0.37)                                 | -11.17 (0.36)                  | success          |
|          127 | 08:17:06.4 +22:22:00                                    | 0.39774 (3)  | 10.60 (0.13)                         | -0.20 (0.36)                                | -10.77 (0.38)                  | success          |
|          128 | 08:17:55.0, +22:09:51                                   | 0.40158 (10) | 10.92 (0.09)                         | 1.10 (0.22)                                 | -9.82 (0.26)                   | success          |

## REFERENCES

Afanasiev, V. L., &amp; Moiseev, A. V. 2011, Baltic Astronomy, 20, 363, doi: 10.1515/astro-2017-0305 Ahumada, R., Prieto, C. A., Almeida, A., et al. 2020, ApJS, 249, 3, doi: 10.3847/1538-4365/ab929e Alexander, D. M., Davis, T. M., Chaussidon, E., et al. 2023, AJ, 165, 124, doi: 10.3847/1538-3881/acacfc Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&amp;A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝cz, B. M., o et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Baker, A. J., Blyth, S., Holwerda, B. W., &amp; LADUMA Team. 2018, in American Astronomical Society Meeting Abstracts, Vol. 231, American Astronomical Society Meeting Abstracts #231, 231.07 Barnes, D. G., Staveley-Smith, L., de Blok, W. J. G., et al. 2001, MNRAS, 322, 486, doi: 10.1046/j.1365-8711.2001.04102.x Behroozi, P., &amp; Silk, J. 2018, MNRAS, 477, 5382, doi: 10.1093/mnras/sty945 Bera, A., Kanekar, N., Chengalur, J. N., &amp; Bagla, J. S. 2023, ApJL, 950, L18, doi: 10.3847/2041-8213/acd0b3 Blyth, S., Baker, A. J., Holwerda, B., et al. 2016, in MeerKAT Science: On the Pathway to the SKA, 4

Brodzeller, A., Dawson, K., Bailey, S., et al. 2023, AJ, 166, 66, doi: 10.3847/1538-3881/ace35d

Catinella, B., &amp; Cortese, L. 2015, MNRAS, 446, 3526, doi: 10.1093/mnras/stu2241

Catinella, B., Saintonge, A., Janowiecki, S., et al. 2018, MNRAS, 476, 875, doi: 10.1093/mnras/sty089 Chakraborty, A., &amp; Roy, N. 2023, MNRAS, 519, 4074, doi: 10.1093/mnras/stac3696 Chaussidon, E., Y`che, C., Palanque-Delabrouille, N., et al. e 2023, ApJ, 944, 107, doi: 10.3847/1538-4357/acb3c2 Chowdhury, A., Kanekar, N., &amp; Chengalur, J. N. 2022, ApJL, 941, L6, doi: 10.3847/2041-8213/ac9d8a Chowdhury, A., Kanekar, N., Das, B., Dwarakanath, K. S., &amp; Sethi, S. 2021, ApJL, 913, L24, doi: 10.3847/2041-8213/abfcc7 Cluver, M. E., Jarrett, T. H., Kraan-Korteweg, R. C., et al. 2010, ApJ, 725, 1550, doi: 10.1088/0004-637X/725/2/1550 Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693, doi: 10.1086/300337 Cooper, A. P., Koposov, S. E., Allende Prieto, C., et al. 2023, ApJ, 947, 37, doi: 10.3847/1538-4357/acb3c0 of Galaxies - SED 2011, ed. R. J. Tuffs &amp; C. C. Popescu, da Cunha, E., Charlot, S., Dunne, L., Smith, D., &amp; Rowlands, K. 2012, in The Spectral Energy Distribution Vol. 284, 292-296, doi: 10.1017/S1743921312009283

da Cunha, E., Charlot, S., &amp; Elbaz, D. 2008, MNRAS, 388, 1595, doi: 10.1111/j.1365-2966.2008.13535.x

DESI Collaboration, Aghamousa, A., Aguilar, J., et al. 2016a, arXiv e-prints, arXiv:1611.00036,

- doi: 10.48550/arXiv.1611.00036
- -. 2016b, arXiv e-prints, arXiv:1611.00037, doi: 10.48550/arXiv.1611.00037
- DESI Collaboration, Abareshi, B., Aguilar, J., et al. 2022, AJ, 164, 207, doi: 10.3847/1538-3881/ac882b

DESI Collaboration, Adame, A. G., Aguilar, J., et al.

2023a, arXiv e-prints, arXiv:2306.06307, doi: 10.48550/arXiv.2306.06307

- -. 2023b, arXiv e-prints, arXiv:2306.06308,

doi: 10.48550/arXiv.2306.06308

- Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168, doi: 10.3847/1538-3881/ab089d
- Dunning, A., Bowen, M., Castillo, S., et al. 2017, in 2017 XXXIInd General Assembly and Scientific Symposium of the International Union of Radio Science (URSI GASS), 1-4, doi: 10.23919/URSIGASS.2017.8105012
- Fern´ndez, X., Gim, H. B., van Gorkom, J. H., et al. 2016, a ApJL, 824, L1, doi: 10.3847/2041-8205/824/1/L1
- Fitzpatrick, E. L. 1999, PASP, 111, 63, doi: 10.1086/316293 Freudling, W., Staveley-Smith, L., Catinella, B., et al. 2011, ApJ, 727, 40, doi: 10.1088/0004-637X/727/1/40
- Giovanelli, R., Haynes, M. P., Kent, B. R., et al. 2005, AJ, 130, 2598, doi: 10.1086/497431
- Guy, J., Bailey, S., Kremin, A., et al. 2023, AJ, 165, 144, doi: 10.3847/1538-3881/acb212
- Hahn, C., Wilson, M. J., Ruiz-Macias, O., et al. 2023, AJ, 165, 253, doi: 10.3847/1538-3881/accff8
- Heintz, K. E., Gim´nez-Arteaga, C., Fujimoto, S., et al. e 2023, ApJL, 944, L30, doi: 10.3847/2041-8213/acb2cf Hoppmann, L. 2014, PhD thesis, University of Western
- Australia
- Hoppmann, L., Staveley-Smith, L., Freudling, W., et al. 2015, MNRAS, 452, 3726, doi: 10.1093/mnras/stv1084 Huang, S., Haynes, M. P., Giovanelli, R., et al. 2014, ApJ, 793, 40, doi: 10.1088/0004-637X/793/1/40
- Jarvis, M., Taylor, R., Agudo, I., et al. 2016, in MeerKAT Science: On the Pathway to the SKA, 6. https://arxiv.org/abs/1709.01901

Koribalski, B. S., Staveley-Smith, L., Westmeier, T., et al.

2020, Ap&amp;SS, 365, 118, doi: 10.1007/s10509-020-03831-4

Kraan-Korteweg, R. C., Staveley-Smith, L., Donley, J.,

Koribalski, B., &amp; Henning, P. A. 2005, in Maps of the

Cosmos, ed. M. Colless, L. Staveley-Smith, &amp; R. A.

Stathakis, Vol. 216, 203, doi: 10.1017/S0074180900196640

- Lan, T.-W., Tojeiro, R., Armengaud, E., et al. 2023, ApJ, 943, 68, doi: 10.3847/1538-4357/aca5fa
- Lelli, F., Fraternali, F., &amp; Sancisi, R. 2010, A&amp;A, 516, A11, doi: 10.1051/0004-6361/200913808

| Levi, M., Bebek, C., Beers, T., et al. 2013, arXiv e-prints, arXiv:1308.0847. https://arxiv.org/abs/1308.0847                            |
|------------------------------------------------------------------------------------------------------------------------------------------|
| Li, D., Wang, P., Qian, L., et al. 2018, IEEE Microwave Magazine, 19, 112, doi: 10.1109/MMM.2018.2802178                                 |
| Meyer, M., Heald, G., & Serra, P. 2009, Proceedings of Panoramic Radio Astronomy. https://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=89 |
| Meyer, M., Robotham, A., Obreschkow, D., et al. 2017, PASA, 34, 52, doi: 10.1017/pasa.2017.31                                            |
| Meyer, M. J., Zwaan, M. A., Webster, R. L., et al. 2004, MNRAS, 350, 1195, doi: 10.1111/j.1365-2966.2004.07710.x                         |
| Moustakas, J., Lang, D., Dey, A., et al. 2023, ApJS, 269, 3, doi: 10.3847/1538-4365/acfaa2                                               |
| Myers, A. D., Moustakas, J., Bailey, S., et al. 2023, AJ, 165, 50, doi: 10.3847/1538-3881/aca5f9                                         |
| Oesch, P. A., Brammer, G., van Dokkum, P. G., et al. 2016, ApJ, 819, 129, doi: 10.3847/0004-637X/819/2/129                               |
| Oke, J. B., Cohen, J. G., Carr, M., et al. 1995, PASP, 107, 375, doi: 10.1086/133562                                                     |
| Perley, D. A. 2019, PASP, 131, 084503, doi: 10.1088/1538-3873/ab215d                                                                     |
| Prochaska, J., Hennawi, J., Westfall, K., et al. 2020a, The Journal of Open Source Software, 5, 2308, doi: 10.21105/joss.02308           |
| Prochaska, J. X., Hennawi, J., Cooke, R., et al. 2020b, pypeit/PypeIt: Release 1.0.0, v1.0.0, Zenodo, doi: 10.5281/zenodo.3743493        |
| Raichoor, A., Moustakas, J., Newman, J. A., et al. 2023, AJ, 165, 126, doi: 10.3847/1538-3881/acb213                                     |
| R., Fraternali, F., Oosterloo, T., & van der Hulst, T. 2008, A&A Rv, 15, 189,                                                            |
| E. F., & Finkbeiner, D. P. 2011, ApJ, 737, doi: 10.1088/0004-637X/737/2/103                                                              |
| Schlafly, E. F., Kirkby, D., Schlegel, D. J., et al. 2023, AJ, 166, 259, doi: 10.3847/1538-3881/ad0832                                   |
| C., Glazebrook, K., Nanayakkara, T., et al. 2018, A&A, 618, A85, doi: 10.1051/0004-6361/201833070                                        |
| F., Rodighiero, G., Elson, E., et al. 2022, ApJL, 935, L13, doi: 10.3847/2041-8213/ac85ae                                                |
| Verheijen, M., van Gorkom, J. H., Szomoru, A., et al. 2007,                                                                              |
| Westmeier, T., Jurek, R., Obreschkow, D., Koribalski, B. S., & Staveley-Smith, L. 2014, MNRAS, 438,                                      |
| 1176, doi: 10.1093/mnras/stt2266 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al.                                             |
| 2010, AJ, 140, 1868, doi: 10.1088/0004-6256/140/6/1868                                                                                   |
| Sinigaglia, ApJL, 668, L9, doi: 10.1086/522621                                                                                           |
| M. 1968, ApJ, 151, 393, doi:                                                                                                             |
| Schmidt, Schreiber,                                                                                                                      |
| Raichoor et al. 2023 Sancisi,                                                                                                            |
| doi: 10.1007/s00159-008-0010-0                                                                                                           |
| Schlafly,                                                                                                                                |
| 103,                                                                                                                                     |
| 10.1086/149446                                                                                                                           |

Xi, H., Peng, B., Staveley-Smith, L., For, B.-Q., &amp; Liu, B. 2022, PASA, 39, e019, doi: 10.1017/pasa.2022.16

Xi, H., Staveley-Smith, L., For, B.-Q., et al. 2021, MNRAS, 501, 4550, doi: 10.1093/mnras/staa3931

- Zhang, K., Wu, J., Li, D., et al. 2019, Science China Physics, Mechanics, and Astronomy, 62, 959506, doi: 10.1007/s11433-019-9383-y
- Zhou, R., Dey, B., Newman, J. A., et al. 2023, AJ, 165, 58, doi: 10.3847/1538-3881/aca5fb
- Zou, H., Zhou, X., Fan, X., et al. 2017, PASP, 129, 064101, doi: 10.1088/1538-3873/aa65ba
Hongwei Xi Bo Peng Lister Staveley-Smith Bi-Qing For Bin Liu Ru-Rong Chen Lei Yu Dejian Ding Wei-Jian Guo Hu Zou Suijian Xue Jing Wang Thomas G. Brink WeiKang Zheng Alexei V. Filippenko Yi Yang Jianyan Wei Y. Sophia Dai Zi-Jian Li Zizhao He Chengzi Jiang Alexei Moiseev Sergey Kotov
The most distant HI galaxies discovered by the 500 m dish FAST
Neutral hydrogen (HI) is the primary component of the cool interstellar
medium (ISM) and is the reservoir of fuel for star formation. Owing to the
sensitivity of existing radio telescopes, our understanding of the evolution of
the ISM in galaxies remains limited, as it is based on only a few hundred
galaxies detected in HI beyond the local Universe. With the high sensitivity of
the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we carried
out a blind HI search, the FAST Ultra-Deep Survey (FUDS), which extends to
redshifts up to 0.42 and a sensitivity of 50 $\rm \mu Jy \cdot beam^{-1}$.
Here, we report the first discovery of six galaxies in HI at $z>0.38$. For
these galaxies, the FAST angular resolution of $\sim\,4'$ corresponds to a mean
linear size of $\sim1.3\,h_{70}^{-1}\,$Mpc. These galaxies are among the most
distant HI emission detections known, with one having the most massive HI
content ($10^{10.93 \pm 0.04}~h_{70}^{-2}\, \rm M_\odot$). Using recent data
from the DESI survey, and new observations with the Hale, BTA, and Keck
telescopes, optical counterparts are detected for all galaxies within the
3-$\sigma$ positional uncertainty ($0.5\,h_{70}^{-1}\,$Mpc) and $\rm 200\,km
\cdot s^{-1}$ in recession velocity. Assuming that the dominant source of HI is
the identified optical counterpart, we find an evidence of evolution in the HI
content of galaxies over the last 4.2 Gyr. Our new high-redshift HI galaxy
sample provides the opportunity to better investigate the evolution of cool gas
in galaxies. A larger sample size in the future will allow us to refine our
knowledge of the formation and evolution of galaxies.