## Prospects for cosmological research with the FAST array: 21-cm intensity mapping survey observation strategies

Jun-Da Pan, a Peng-Ju Wu, a Guo-Hong Du, a Yichao Li a and Xin Zhang a,b,c, ∗

- a Liaoning Key Laboratory of Cosmology and Astrophysics, College of Sciences, Northeastern University, Shenyang 110819, China
- b MOE Key Laboratory of Data Analytics and Optimization for Smart Industry, Northeast-
- ern University, Shenyang 110819, China
- c National Frontiers Science Center for Industrial Intelligence and Systems Optimization, Northeastern University, Shenyang 110819, China

E-mail: panjunda@stumail.neu.edu.cn, wupengju@stumail.neu.edu.cn, duguohong@stumail.neu.edu.cn, liyichao@mail.neu.edu.cn, zhangxin@mail.neu.edu.cn

Abstract. Precise cosmological measurements are essential for understanding the evolution of the universe and the nature of dark energy. The Five-hundred-meter Aperture Spherical Telescope (FAST), the most sensitive single-dish radio telescope, has the potential to provide the precise cosmological measurements through neutral hydrogen 21 cm intensity mapping sky survey. This paper primarily explores the potential of technological upgrades for FAST in cosmology. The most crucial upgrade begins with equipping FAST with a wide-band receiver (0 &lt; z &lt; 2 5). . This upgrade can enable FAST to achieve higher precision in cosmological parameter estimation than the Square Kilometre Array Phase 1 Mid-Frequency Array. On this basis, expanding to a FAST array (FASTA) consisting of six identical FASTs would offer significant improvements in precision compared to FAST. Additionally, compared with the current results from the data combination of cosmic microwave background, baryon acoustic oscillations (optical galaxy surveys), and type Ia supernovae, FASTA can provide comparable constraints. Specifically, for the dark-energy equation-of-state parameters, FASTA can achieve σ w ( 0 ) = 0 09 and . σ w ( a ) = 0 33. .

∗ Corresponding author.

## Contents

| 1 Introduction   | 1 Introduction         | 1 Introduction                        |   1 |
|------------------|------------------------|---------------------------------------|-----|
| 2                | 21 cm intensity        | mapping                               |   3 |
|                  | 2.1                    | Signal                                |   3 |
|                  | 2.2                    | Noise                                 |   3 |
|                  | 2.3                    | Foreground                            |   4 |
|                  | 2.4                    | Instrument parameters                 |   4 |
|                  | 2.5                    | The Fisher matrix                     |   5 |
|                  | 2.6                    | Cosmological parameter constraints    |   6 |
| 3                | Results and discussion | Results and discussion                |   6 |
|                  | 3.1                    | Combination of FAST and SKA-Mid       |   6 |
|                  | 3.2                    | Combination of FAST and 40-m antennas |   7 |
|                  | 3.3                    | FASTA                                 |   8 |
|                  | 3.4                    | Phased array feed                     |   9 |
| 4                | Conclusion             | Conclusion                            |  11 |

## 1 Introduction

Dark energy is a component of the universe that causes the late-time cosmic acceleration. Although it is known to comprise 68% of the total energy density of the universe, the nature of dark energy remains unknown. The Λ cold dark matter (ΛCDM) model, the standard model of cosmology, treats dark energy as a constant energy density filling space homogeneously. However, recent observations of baryon acoustic oscillations (BAO) from the Dark Energy Spectroscopic Instrument (DESI) have shown deviations from the ΛCDM model [18]. The result suggests that dark energy may be a dynamical energy density evolving over time or indicate the presence of other new physics. Nonetheless, further observations and advanced instruments are needed to confirm these possibilities and to better understand the evolutionary dynamics of dark energy.

BAO from the large-scale structure (LSS) is a late-time cosmological probe that can be used as a standard ruler to study the expansion history of the universe [9-16]. Currently, BAOstudies primarily rely on optical galaxy surveys [1, 17-20]. While effective, these surveys are limited in terms of efficiency and coverage. To overcome these limitations and achieve more precise cosmological measurements, advancements in radio astronomy have developed a promising method known as 21 cm intensity mapping (21 cm IM). After reionization, neutral hydrogen (H I ) mainly exists within galaxies. This technique measures the intensity of the characteristic 21 cm emission line from H I to study the LSS [21-25]. This method does not require resolving individual galaxies, thus it has high survey efficiency. The characteristics of the 21 cm signal allow us to probe deeper into the universe and directly obtain redshift information. Recent forecasts and observational studies have confirmed the potential of 21 cm IM in cosmology [26-39]. In addition to the 21 cm line, intensity mapping of other spectral lines, such as CO and CII, has also attracted attention. These lines can provide information about the LSS and the distribution of molecular and ionized gas in galaxies [40-42].

To achieve more precise cosmological measurements, advanced radio telescopes are utilized in 21 cm IM. These telescopes include those currently in use, such as the Five-hundredmeter Aperture Spherical Telescope (FAST) [43], the MeerKAT [44-46], the Tianlai [47-52], the Canadian Hydrogen Intensity Mapping Experiment (CHIME) [53]. Currently, 21 cm signals have been observed through cross-correlation with galaxy surveys [22, 24, 46, 54-57], and the MeerKAT telescope has also detected the autocorrelation spectrum on Mpc scales [58]. However, the large scale autocorrelation power spectrum still has not been observed [23]. Future instruments planned for use include the Square Kilometre Array (SKA) [5963], the Baryon Acoustic Oscillations from Integrated Neutral Gas Observations (BINGO) [64, 65], the Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) [66], and the Canadian Hydrogen Observatory and Radio-transient Detector (CHORD) [67]. These advanced facilities are expected to further enhance the capabilities of 21 cm IM, enabling more detailed and accurate cosmological studies.

Among these, FAST is particularly significant due to its large aperture and sensitivity. The current FAST has achieved substantial success in many areas, particularly in fast radio burst and gravitational wave researches [68-70]. However, to harness its potential for 21 cm IM and to extend its capability for cosmological studies, several upgrades are planned for FAST [71-74].

This paper aims to explore how potential future upgrades to FAST could promote cosmological research by using 21 cm IM. For the expansion of FAST itself, equipping it with a wide-band receiver aims to extend its 21 cm signal detection capability, enabling broader redshift coverage [71]. This will allow us to observe a larger volume of the universe, aiding in constructing a more precise three-dimensional map of the LSS. Moreover, a phased array feed (PAF) is planned to enhance its survey efficiency [75, 76]. Additionally, the expansion of the construction around FAST is also considered. The plan involves building 40-meter aperture radio telescopes (40-m antennas) in two phases around FAST: the first phase involves constructing six 40-m antennas, and the second phase will increase this number to twentyfour [74]. The most significant upgrade is the FAST array (FASTA), which plans to add five more identical 500-meter diameter telescopes to join the current FAST [72]. This network will significantly improve sensitivity and survey efficiency, enhancing the ability to explore LSS and providing higher precision in cosmological parameter estimation.

In this paper, we present a forecast for estimating cosmological parameters using simulated 21 cm IM data and discuss how the various extensions of FAST will improve the precision of these constraints.

In Section 2, we introduce the 21 cm IM simulation and the methods used to constrain cosmological parameters. In Section 3, we discuss the performance of various configurations for 21 cm IM and the resulting cosmological constraints. In Section 4, we conclude with perspectives for future research.

In our analysis, we adopt the Planck best-fit ΛCDM model as fiducial cosmology [77]. The parameters are H 0 = 67 3 km s . -1 Mpc -1 , Ω m = 0 317, Ω . b = 0 0495, Ω . K = 0, σ 8 = 0 812, . and n s = 0 965. .

## 2 21 cm intensity mapping

## 2.1 Signal

21 cm IM measures the intensity of 21 cm signals and converts it into brightness temperature. The mean temperature of H I can be expressed as [26]

$$\overline { T } _ { \text{b} } = \frac { 3 } { 3 2 \pi } \frac { \tilde { \not c } ^ { 3 } A _ { 1 0 } } { k _ { \text{B} } m _ { \text{p} } \nu _ { 2 1 } ^ { 2 } } \frac { ( 1 + z ) ^ { 2 } } { H ( z ) } \Omega _ { \text{H} 1 } ( z ) \rho _ { \text{c,0} }$$

where h is the Planck constant, c is the speed of light, A 10 is the Einstein coefficient for spontaneous emission [21], k B is the Boltzmann constant, m p is the proton mass, ν 21 is the H I line frequency, H z ( ) is the Hubble parameter, Ω H I ( z ) is the fractional density of H I , and ρ c 0 , is the current critical density of the universe.

The signal covariance is expressed in terms of the power spectrum [9, 26]

$$C ^ { \text{S} } ( q, y ) = \frac { T _ { b } ^ { 2 } ( z _ { i } ) } { r ^ { 2 } r _ { \nu } } ( b _ { \text{H1} } + f \mu ^ { 2 } ) ^ { 2 } \exp ( - k ^ { 2 } \mu ^ { 2 } \sigma _ { \text{NL} } ^ { 2 } ) P ( k, z ),$$

where q = k ⊥ r and y = k r ∥ ν . r denotes the comoving distance to z , and r ν corresponds to the comoving distance for the frequency interval. b H I represents the bias factor of H I [78]. f is the linear growth rate of the structure. µ = k /k ∥ is the cosine of the angle between the wavevector k and the line of sight. σ NL = 7 Mpc represents the non-linear dispersion scale [79]. P k, z ( ) is the matter power spectrum at wavenumber k and redshift z , which can be calculated by CAMB [80].

## 2.2 Noise

The noise covariance can be written as [26]

$$C ^ { \mathbb { N } } ( q, y ) = \frac { T _ { s y s } ^ { 2 } } { t _ { t o t } \ \Delta \nu \ N _ { \mathrm b } N _ { \mathrm d } } U _ { \mathrm b } \ B _ { \perp } ^ { - 2 } B _ { \| } ^ { - 1 },$$

where T sys includes the instrumental temperature T inst and the sky background radiation T sky ≈ 60 K × ( ν/ 300 MHz) -2 5 . . t tot represents the total integration time, N b is the number of beams, N d is the number of dishes in the array, and ∆ ν is the total bandwidth. U bin = S area ∆˜ represents the survey volume at a given redshift, where ν S area is the survey area and ∆˜ = ∆ ν ν/ν 21 is the dimensionless bandwidth for a given redshift bin. B ∥ is the radial beam response

$$B _ { \| } ( y ) = \exp \left ( - \frac { \left ( y \delta \nu / \nu _ { 2 1 } \right ) ^ { 2 } } { 1 6 \ln 2 } \right ),$$

where δν is the bandwidth of an individual frequency channel. B ⊥ is the transverse beam response in single dish mode

$$B _ { \perp } ( q ) = \exp \left ( - \frac { \left ( q \frac { \lambda } { D _ { \text{dish} } } \right ) ^ { 2 } } { 1 6 \ln 2 } \right ),$$

where D dish is the diameter of the dish.

## 2.3 Foreground

The foreground is an important component of the signal received by the telescope. It primarily originates from our own galaxy and extragalactic point sources. At this stage, the effectiveness of methods for subtracting foregrounds from observational data has been very limited. The methods include Polynomial Fitting, Principal Component Analysis, Independent Component Analysis, and others [81-92]. We assume a method that can effectively subtract the foreground, and the residual foreground model can be written as [26, 93]

$$C ^ { \text{F} } ( q, y ) = \varepsilon _ { \text{FG} } ^ { 2 } \sum _ { X } A _ { X } \left ( \frac { l _ { p } } { 2 \pi q } \right ) ^ { n _ { X } } \left ( \frac { \nu _ { p } } { \nu _ { i } } \right ) ^ { m _ { X } },$$

where ε FG is a factor that describes the efficiency of foreground subtraction. ε FG = 1 represents no foreground subtraction, while ε FG = 0 represents perfect foreground subtraction. We assume ε FG = 10 -6 , as this level is effective for extracting cosmological signals. X is the varies foreground component shown in Table 1.

Table 1 : Foreground model parameters at l p = 1000 and ν p = 130 MHz [93].

| Foreground                  |   A X [mK 2 |   n X |   m X |
|-----------------------------|-------------|-------|-------|
| Extragalactic point sources |      57     |   1.1 |  2.07 |
| Extragalactic free-free     |       0.014 |   1   |  2.1  |
| Galactic synchrotron        |     700     |   2.4 |  2.8  |
| Galactic free-free          |       0.088 |   3   |  2.15 |

## 2.4 Instrument parameters

In this paper, we consider several different configurations of FAST and the Square Kilometre Array Phase 1 Mid-Frequency Array (SKA1-Mid) for simulating the 21 cm IM. The specific instrument parameters are listed in Table 2.

The current capability of FAST to detect 21 cm signals covers the redshift range 0 &lt; z &lt; 0 35. . A new wide-band receiver under test extends this range to 0 &lt; z &lt; 1 84 [71]. . Currently, DESI measurements have mapped the cosmological distance-redshift relation with an accuracy of 1%-3% across redshifts below z &lt; 2 5 . [2]. This paper assumes that the upgraded new wide-band receiver on FAST can detect 21 cm signals up to 0 &lt; z &lt; 2 5. .

Additionally, constructing 6 × 40-m antennas or 24 × 40-m antennas around FAST is considered. These dishes will use the same receiver and feed system as FAST, so their frequency bands and N b will be identical. Given the outstanding performance of FAST, we also consider constructing five additional identical FAST telescopes to form FASTA, significantly improving survey efficiency. Furthermore, we consider equipping FAST with a PAF to increase the value of N b to 100.

Regarding instrument temperature, in Ref. [94] the authors measured the temperature of FAST near boresight to be 20 K in the L-band. We assume that the frequency band can be broadened while maintaining the instrument temperature at 20 K. For the PAF, we assume the temperature to be approximately 30 K due to the increase in the number of beams. The total observation time is assumed to be 10000 hours. The calculation for the survey area is as follows

$$S _ { \text{area} } = \eta \int _ { \phi } \int _ { \theta } \sin \theta \text{d} \theta \text{d} \phi$$

where η represents the fraction of the initial survey area that remains usable for scientific observation after excluding the Milky Way. For FAST/FASTA, we consider η = 0 7 and the . range of declination from -14 ◦ to +66 , yielding an area of approximately 8100 deg . ◦ 2

For SKA1-Mid, we consider Band 1 (0 35 . &lt; z &lt; 3) planned for 21 cm IM, with the survey area of 20000 deg 2 and the total observation time of 10000 hours. The instrument temperature is modeled based on actual measurements and simulation data from the SKA Observatory, and the form is as follows [62, 95]

$$T _ { \text{inst} } = 1 5 \, K + 3 0 \, K \left ( \frac { f } { \text{GHz} } - 0. 7 5 \right ) ^ { 2 }.$$

Table 2 : Instrument parameters for the 21 cm IM forecast.

| Instrument    |   z min |   z max | N d     |   N b |   D dish [m] |   S area [deg 2 |   t tot [h] | T inst [K]   |
|---------------|---------|---------|---------|-------|--------------|-----------------|-------------|--------------|
| FAST          |    0    |     2.5 | 1       |    19 |          300 |            8100 |       10000 | 20           |
| SKA1-Mid      |    0.35 |     3   | 197     |     1 |           15 |           20000 |       10000 | 18-28        |
| 40-m antennas |    0    |     2.5 | 6 or 24 |    19 |           40 |            8100 |       10000 | 20           |
| FASTA         |    0    |     2.5 | 6       |    19 |          300 |            8100 |       10000 | 20           |
| FAST(PAF)     |    0    |     2.5 | 1       |   100 |          300 |            8100 |       10000 | 30           |

## 2.5 The Fisher matrix

To achieve optimal observational results, different telescopes require different survey strategies, so the methods for calculating C tot also vary. For a single dish or identical dishes observing the same area of the sky (such as FAST or FASTA), the total covariance matrix is expressed as

$$C ^ { t o t } = C ^ { S } + C ^ { N } + C ^ { F }.$$

For the combination of different telescopes, we consider two scenarios. The first scenario involves two telescopes in close geographic proximity observing the same area of the sky, such as combining 40-m antennas and FAST. This approach maximizes the reduction of thermal noise impact on observations. Since C S and C F are only related to the observed sky area, they remain fixed. Different telescopes have varying C N due to their differing thermal noise integration effects. Therefore, the total covariance matrix can be expressed as

$$C ^ { t o t } = C ^ { S } + \left ( ( C _ { 1 } ^ { N } ) ^ { - 1 } + ( C _ { 2 } ^ { N } ) ^ { - 1 } \right ) ^ { - 1 } + C ^ { F },$$

where C N 1 and C N 2 are the noise covariance matrices of the different dishes, with ( C N ) -1 denoting their inverse matrices. The second scenario involves two telescopes with significant differences in latitude, resulting in a very small common observable sky area. Therefore, they must observe different areas of the sky, which also maximizes the survey volume. Since both the observed sky areas and the telescopes are different, the total covariance matrix is given by

$$C ^ { t o t } & = \left ( ( C _ { 1 } ^ { t o t } ) ^ { - 1 } + ( C _ { 2 } ^ { t o t } ) ^ { - 1 } \right ) ^ { - 1 },$$

where C tot 1 and C tot 2 are the total covariance matrices of the different telescopes.

After obtaining the total covariance matrix C tot , we can transform it into the Fisher matrix of observables { p i } [26]

$$F _ { i j } = \frac { 1 } { 2 } U _ { \text{bin} } \int \frac { \mathrm d ^ { 2 } q \, \mathrm d y } { ( 2 \pi ) ^ { 3 } } \left [ \frac { \partial \ln C ^ { \mathrm t o t } ( q, y ) } { \partial p _ { i } } \frac { \partial \ln C ^ { \mathrm t o t } ( q, y ) } { \partial p _ { j } } \right ].$$

In this paper, we mainly calculate the Fisher matrix for several key observables to constain the evolution of LSS, as well as the growth and geometry of the universe. These include the transverse distance indicator D A ( z ), radial distance indicator H z ( ) and the parameter combination fσ 8 ( z ) from redshift space distortion (RSD), where f ( z ) is growth rate and σ 8 is the amplitude of matter density fluctuations on a scale of 8 h -1 Mpc (with h the dimensionless Hubble constant).

## 2.6 Cosmological parameter constraints

We use the Fisher matrix to propagate errors from observational errors to the BAO parameters D A ( z ), H z ( ), and fσ 8 ( z ). Then, we use the errors of BAO parameters to form the likelihood function for cosmological models. Finally, we employ Markov Chain Monte Carlo (MCMC) analysis to constrain the cosmological parameters [96]. Our considered parameter set is { H , 0 Ω m , Ω b , σ 8 , n s , w 0 , w a } . We assume a flat universe with Ω K = 0 and fix the total neutrino mass ∑ m ν = 0 06 eV. . For the free parameters, we assume the priors within the following ranges: H 0 ∈ [50 , 80] km s -1 Mpc -1 , Ω m ∈ [0 , 0 7], . Ω b ∈ [0 03 . , 0 07], . σ 8 ∈ [0 , 2], n s ∈ [0 8 . , 1 2], . w 0 ∈ -[ 3 5], and , w a ∈ -[ 5 5]. , Although BAO alone cannot constrain n s and Ω , we consider them as free parameters to maintain consistency of the parameter space with b the cases when considering CMB and supernova (SN) data.

To compare with the mainstream observations, we calculated the constraints from the latest observations. The data used include: CMB Planck 2018 highℓ TT TE EE , , , and lowℓ TT EE , , and PR4 lensing data [97-101]; BAO - DESI data [2]; and SN - Pantheon+ data [102]. The data combination of CMB, BAO, and SN in this paper is also abbreviated as CBS.

The equation of state (EoS) of dark energy is given by w = p/ρ , where p is the pressure and ρ is the energy density. We employ different models to describe its dynamics. (i) ΛCDM: w = -1. (ii) w CDM: w z ( ) = w , where w is a constant EoS for dark energy. (iii) w w 0 a CDM: w varies with redshift according to the parameterization w z ( ) = w 0 + w z/ a (1 + z ), where w 0 is the present-day value of the EoS parameter and w a describes its evolutionary behavior.

## 3 Results and discussion

In this section, we first discuss the capabilities of various instruments and their combinations in 21 cm IM using the relative errors of BAO parameters { ε D ( A ( z )) , ε ( H z ( )) , ε ( fσ 8 ( z )) } , where ε ξ ( ) = σ ξ /ξ ( ) . We then present the cosmological parameter constraints employing FAST, SKA1-Mid, FASTA, and CBS, followed by comparison and discussion of the results.

## 3.1 Combination of FAST and SKA-Mid

In this subsection, we discuss the capabilities of SKA1-Mid and FAST with a wide-band receiver, as well as their combination in 21 cm IM. Due to geographical factors, the overlapping observable sky area for FAST and SKA1-Mid is very limited, which necessitates a strategy where each telescope observes different sky regions.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-bf13100f-ced6-44ed-9eed-53ea87b22810.jpg)

z

z

z

Figure 1 : Relative errors of D A ( z ), H z ( ), and fσ 8 ( z ) for 21 cm IM with FAST, SKA1-Mid, and FAST+SKA1-Mid.

It should be noted that for SKA1-Mid, we consider the single dish mode. This is because, in interferometer mode, the maximum observable scale is determined by the shortest baseline. Although SKA1-Mid has an advantage in angular resolution in interferometer mode, it lacks very short baselines, which limits its ability to observe large angular scales. Its sparse array design results in low baseline density and does not achieve good uv coverage, making it unsuitable for 21 cm IM. Improving its performance in interferometer mode would require a denser antenna array, like the CHIME design. In single dish mode, the maximum observable scale is determined by the aperture of the telescope, allowing SKA1-Mid to observe larger scales. Single dish mode also offers other advantages: a wider field of view, more beams, enabling broader sky coverage within limited observation time. Therefore, single dish mode is an optimal choice for 21 cm IM on cosmological scales [26, 37].

From Figure 1, it is evident that FAST and SKA1-Mid have distinct advantages across different redshift ranges. SKA1-Mid achieves tight constraints within the range of 0 35 . &lt; z &lt; 0 8 due to its high survey efficiency. . In this range, the relative errors in BAO parameters are typically below 3% (see also Ref. [33]). On the other hand, FAST attains robust constraints within the range of 0 8 . &lt; z &lt; 2 5 because of its high transverse resolution and sensitivity. . In this range, the relative errors in BAO parameters are typically around or below 4%.

Their combined observations enable precise measurement of BAO parameters across the redshift range of 0 &lt; z &lt; 2 5, achieving constraints on . D A ( z ), H z ( ), and fσ 8 ( z ) to below 2% for most of this range, and to as low as 1% at z ≈ 1.

Upgrading the existing receiver to detect 21 cm signals up to 0 &lt; z &lt; 2 5 will enable . FAST to utilize its high-redshift advantage in 21 cm IM. Combined with SKA1-Mid, this complementary strategy will improve measurements of BAO parameters across a wide redshift range.

## 3.2 Combination of FAST and 40-m antennas

In this subsection, we discuss the extent to which incorporating 40-m antennas with FAST enhances the constraints on BAO parameters in 21 cm IM. An array of 40-m antennas, located in close geographical proximity to FAST, enhances the survey efficiency by enabling observations of the same sky regions.

Figure 2 illustrates the fractional constraints achievable on D A ( z ), H z ( ), and fσ 8 ( z ) using 21 cm IM with an array of 40-m antennas, FAST, and their combined configuration.

Figure 2 : Relative errors of D A ( z ), H z ( ), and fσ 8 ( z ) for 21 cm IM with FAST, 40-m antennas, FAST+6 × 40-m antennas, and FAST+24 × 40-m antennas.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-80b686dc-d6ee-431a-a2db-d65b36bdb0a1.jpg)

Table 3 : The errors for the cosmological parameters in the ΛCDM, w CDM, and w w 0 a CDM models using CBS and 21 cm IM. Note that here H 0 is in units of km s -1 Mpc -1 .

| Data     | ΛCDM              | ΛCDM               | w CDM             | w CDM              | w CDM            | w 0 w a CDM       | w 0 w a CDM        | w 0 w a CDM        | w 0 w a CDM        |
|----------|-------------------|--------------------|-------------------|--------------------|------------------|-------------------|--------------------|--------------------|--------------------|
| Data     | σ (Ω m ) / 10 - 3 | σ ( H 0 ) / 10 - 1 | σ (Ω m ) / 10 - 3 | σ ( H 0 ) / 10 - 1 | σ ( w ) / 10 - 2 | σ (Ω m ) / 10 - 3 | σ ( H 0 ) / 10 - 1 | σ ( w 0 ) / 10 - 2 | σ ( w a ) / 10 - 1 |
| CBS      | 4.6               | 3.4                | 6.5               | 6.4                | 2.4              | 6.4               | 6.5                | 6.8                | 3.0                |
| SKA1-Mid | 8.7               | 7.0                | 9.0               | 9.1                | 4.5              | 26.0              | 15.0               | 16.0               | 7.3                |
| FAST     | 7.6               | 5.1                | 7.8               | 7.2                | 3.6              | 21.0              | 14.0               | 14.0               | 5.7                |
| FASTA    | 3.6               | 2.6                | 3.7               | 3.9                | 2.0              | 12.0              | 8.9                | 9.0                | 3.3                |

The 40-m antennas perform better at low redshifts because they have higher survey efficiency. FAST is more effective at higher redshifts due to its greater sensitivity and resolution. The combination of FAST and 40-m antennas significantly reduces observational errors, particularly at high redshifts. The relative errors in BAO parameters are typically below 4%, with errors reaching as low as 2% at z ≈ 1.

For the configurations of 6 × 40-m antennas and 24 × 40-m antennas combined with FAST, Figure 2 shows that the larger array provides only marginal improvements. This result suggests that the combination of FAST with 6 × 40-m antennas offers a more optimal configuration for the 21 cm IM survey.

## 3.3 FASTA

In this subsection, we present the cosmological parameters constraints by using simulations of 21 cm IM data with different instruments and current observational data from CBS. Figure 3 shows the parameter constraints within the 68.3% and 95.4% confidence intervals for the ΛCDM, w CDM, and w w 0 a CDM models. The specific values are listed in Table 3.

Compared to the results in Ref. [33] for FAST (0 &lt; z &lt; 0 35), . a wide-band receiver proposed in this paper utilizes the resolution advantage of the FAST at high redshifts, significantly enhancing the constrains on cosmological parameters. The main specific parameter constraints are as follows: for the ΛCDM model, FAST achieves σ (Ω m ) = 0 0076 and . σ H ( 0 ) = 0 51 km s . -1 Mpc -1 . For the w CDM model, the constraint is σ w ( ) = 0 036. . For the w w 0 a CDM model, the constraints are σ w ( 0 ) = 0 14, and . σ w ( a ) = 0 57. . Compared to SKA1Mid, FAST also has higher sensitivity and lower receiver noise, which enables it to achieve

Figure 3 : Constraints (68.3% and 95.4% confidence level) on the dark energy models using simulated data of 21 cm IM with SKA1-Mid, FAST, FASTA, and observational data of CBS. Left, middle, and right panels show the constraints on the ΛCDM, w CDM, and w w 0 a CDM models, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-ccea56cd-e9b5-4695-84ee-d1ad98231107.jpg)

better parameter constraints. For the dark-energy EoS parameters, there is an improvement of approximately 20% for w , 12.5% for w 0 , and 21.9% for w a .

Further advances are observed with FASTA. From Figure 4, FASTA performs well in observing BAO parameters across the entire band, achieving an observational error as low as 1%. This precise observation of BAO parameters translates into tight constraints on cosmological parameters.

In the ΛCDM model, FASTA provides best constraints, σ (Ω m ) = 3 6 . × 10 -3 and σ H ( 0 ) = 0 26 km s . -1 Mpc -1 , which surpasses other instruments and CBS. In the w CDM model, FASTA also provides best constraints, σ (Ω m ) = 0 0037, . σ H ( 0 ) = 0 39 km s . -1 Mpc -1 , and σ w ( ) = 0 02. . For the parameter w , FASTA offers an improvement of approximately 55.56% compared to SKA1-Mid, 44.44% compared to FAST, and 16.67% compared to CBS. In the w w 0 a CDMmodel, FASTA achieves highly precise parameter constraints, with σ (Ω m ) = 0 012, . σ H ( 0 ) = 0 89 km s . -1 Mpc -1 , σ w ( 0 ) = 0 09, and . σ w ( a ) = 0 33. For . w 0 and w a , FASTA provides an improvement of approximately 43.75% and 54.79% compared to SKA1-Mid, and 35.71% and 42.11% compared to FAST, while achieving constraints comparable to CBS using 21 cm IM alone.

In summary, FASTA significantly improves cosmological parameter constraints across various models, surpassing other instruments and achieving results comparable to CBS using 21 cm IM alone.

## 3.4 Phased array feed

In this subsection, we discuss the enhancements and challenges of upgrading the feed of FAST to PAF for 21 cm IM. PAF can achieve various beamforming modalities by adjusting the phase and amplitude of each element in the array, thereby generating multiple beams simultaneously. This capability covers a larger field of view (FOV) and enables rapidly beam steering and scanning [75, 76]. At the same time, different beamforming strategies significantly impact noise characteristics, and data calibration and processing. We considered two beamforming strategies in terms of survey efficiency.

In the first scenario, we consider a fixed FOV and maximize the N b within this field. The size of a beam is given by θ = c/ νD ( dish ). If there is no overlap between the beams, the

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-16e28e90-80e2-4b87-bed5-5ffb557b9b12.jpg)

Figure 4 : Relative errors of D A ( z ), H z ( ), and fσ 8 ( z ) for 21 cm IM with SKA1-Mid, FAST, and FASTA.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-22ef7f2b-826e-41f0-9067-23da39b741d3.jpg)

z

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-9243ad64-6052-44e7-aa4a-de8ef4392dde.jpg)

z

Figure 5 : Relative errors of D A ( z ), H z ( ), and fσ 8 ( z ) for 21 cm IM with FAST and FAST(PAF), where FAST(PAF)-1 represents the fixed FOV scenario and FAST(PAF)-2 represents the fixed N b scenario.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00268v2-figure-3ab69a8a-b2ac-4c03-8d14-befe9bcf5fcc.jpg)

maximum N b can be written as

$$N _ { \mathrm b } = \text{FOV} \cdot \left ( \frac { \nu D } { c } \right ) ^ { 2 }.$$

In this case, C N follows the same pattern as described in Ref. [26]

$$C ^ { \mathbb { N } } ( q, y ) \to C ^ { \mathbb { N } } ( q, y ) \times \begin{cases} 1, & \nu > \nu _ { \text{crit} } \\ \left ( \nu _ { \text{crit} } / \nu \right ) ^ { 2 }, & \nu \leqslant \nu _ { \text{crit} } \end{cases} \,.$$

Here we have ν crit = 1420 MHz. As shown in Figure 5, the errors increase significantly at low frequencies. This is because at low frequencies, the beam size increases, resulting in a decrease in N b .

Another scenario is a fixed N b . We assume the case of 100 beams. As shown in Figure 5, the increase in N b significantly enhances survey efficiency, achieving best constraint result of 2% around z = 1. However, in practice, the FOV does not increase endlessly as the beam size increases at lower frequencies. As the FOV certainly deviates from the boresight, the survey efficiency decreases. This requires actual measurements to PAF, but here we make

an ideal estimate. The real situation is expected to lie between the fixed FOV and the fixed N b scenarios.

## 4 Conclusion

In this paper, we explore several potential upgrades to FAST using 21 cm IM to enhance cosmological research. Those upgrades include FAST with a wide-band receiver, the combination of 40-m antennas and FAST, FASTA, and FAST with a PAF. We simulated the full 21 cm IM power spectrum to constrain the BAO and RSD parameters D A ( z ), H z ( ), and fσ 8 ( z ), and used these results to evaluate the capabilities of these different configurations. Finally, we used the simulated BAO and RSD data to constrain the cosmological parameters of the ΛCDM, w CDM, and w w 0 a CDM models.

Firstly, the wide-band receiver utilizes the capabilities of FAST in the redshift range 0 8 . &lt; z &lt; 2 5, achieving better cosmological parameter constraints than SKA1-Mid. . It also enables FAST to form a complementary observation with SKA1-Mid, which has strengths in the redshift range of 0 &lt; z &lt; 0 8, . providing high-precision observations across a wider redshift range.

Secondly, we discuss the combination of the 40-m antennas and FAST. Incorporating 40-m antennas with FAST significantly enhances survey efficiency, thereby reducing BAO parameter errors to typically below 2%. For the combination with FAST, 6 × 40-m antennas are a more balanced configuration compared to 24 × 40-m antennas for 21 cm IM.

Next, we discuss the cosmological parameter constraints for FASTA. Through BAO parameter measurements approaching 1% accuracy, it shows excellent performance in constraining cosmological parameters. Specifically, in the w w 0 a CDM model, FASTA achieves highly precise parameter constraints, with improvements of 43.75% for w 0 and 54.79% for w a over SKA1-Mid, and 35.71% for w 0 and 42.11% for w a over FAST, achieving results comparable to CBS using 21 cm IM alone.

Lastly, we analyze FAST with a PAF. We consider two beamforming strategies, a fixed FOV and a fixed N b . The real situation is likely to fall between these two cases. This provides a basis for understanding the advantages and limitations of different beamforming strategies for 21 cm IM. A balanced beamforming strategy will be beneficial for effectively utilizing PAF.

In summary, these upgrades to FAST have enormous potential to improve the precision of cosmological measurements. Despite facing some technical challenges, they can greatly enhance our understanding of the evolution of the universe and the nature of dark energy.

## Acknowledgments

We thank Guang-Chen Sun, Ji-Guo Zhang, Zi-Qiang Zhao, and Lu Feng for helpful discussions and suggestions. We are grateful for the support from the National SKA Program of China (Grants Nos. 2022SKA0110200 and 2022SKA0110203), the National Natural Science Foundation of China (Grants Nos. 12473001, 11975072, 11875102, 11835009, and 12305069), the National 111 Project (Grant No. B16009), and the Program of the Education Department of Liaoning Province (Grant No. JYTMS20231695).

## References

- [1] DESI collaboration, The DESI experiment part I: science, targeting, and survey design , 1611.00036 .
- [2] DESI collaboration, DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars , 2404.03000 .
- [3] DESI collaboration, DESI 2024 VI: Cosmological constraints from the measurements of baryon acoustic oscillations , 2404.03002 .
- [4] N. Craig, D. Green, J. Meyers and S. Rajendran, No ν s is Good News , Journal of High Energy Physics 09 (2024) 097.
- [5] D. Green and J. Meyers, The Cosmological Preference for Negative Neutrino Mass , 2407.07878 .
- [6] P.-J. Wu and X. Zhang, Measuring cosmic curvature with non-CMB observations , 2411.06356 .
- [7] G.-H. Du, P.-J. Wu, T.-N. Li and X. Zhang, Impacts of dark energy on weighing neutrinos after DESI BAO , 2407.15640 .
- [8] T.-N. Li, P.-J. Wu, G.-H. Du, S.-J. Jin, H.-L. Li, J.-F. Zhang et al., Constraints on Interacting Dark Energy Models from the DESI Baryon Acoustic Oscillation and DES Supernovae Data , The Astrophysical Journal 976 (2024) 1.
- [9] C. Blake and K. Glazebrook, Probing Dark Energy Using Baryonic Oscillations in the Galaxy Power Spectrum as a Cosmological Ruler , Astrophysical Journal 594 (2003) 665.
- [10] S. Cole, W.J. Percival, J. Peacock, P. Norberg et al., The 2dF Galaxy Redshift Survey: Power-spectrum analysis of the final dataset and cosmological implications , Monthly Notices of the Royal Astronomical Society 362 (2005) 505.
- [11] W.J. Percival, B. Reid, D. Eisenstein, N. Bahcall et al., Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample , Monthly Notices of the Royal Astronomical Society 401 (2010) 2148.
- [12] F. Beutler, C. Blake, M. Colless, D. Jones et al., The 6dF Galaxy Survey: Baryon Acoustic Oscillations and the Local Hubble Constant , Monthly Notices of the Royal Astronomical Society 416 (2011) 3017.
- [13] C. Blake, S. Brough, M. Colless, C. Contreras et al., The WiggleZ Dark Energy Survey: Testing the Cosmological Model with Baryon Acoustic Oscillations at z=0.6 , Monthly Notices of the Royal Astronomical Society 415 (2011) 2892.
- [14] A.J. Ross, L. Samushia, C. Howlett, W.J. Percival et al., The Clustering of the SDSS DR7 Main Galaxy Sample - I. A 4 Per Cent Distance Measure at z = 0.15 , Monthly Notices of the Royal Astronomical Society 449 (2015) 835.
- [15] S. Alam, F. Beutler, C. Blake, J.E. Bautista et al., The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample , Monthly Notices of the Royal Astronomical Society 470 (2017) 2617.
- [16] M. Ata, F. Beutler, C. Blake, J.E. Bautista et al., The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: First measurement of baryon acoustic oscillations between redshift 0.8 and 2.2 , Monthly Notices of the Royal Astronomical Society 473 (2018) 4773.
- [17] D.G. York, J. Adelman, J. Anderson, S. Anderson et al., The Sloan Digital Sky Survey: Technical Summary , The Astronomical Journal 120 (2000) 1579-1587.

- [18] W.J. Percival, S. Cole, G. Efstathiou, C.S. Frenk et al., The 2dF Galaxy Redshift Survey: The power spectrum and the matter content of the Universe , Monthly Notices of the Royal Astronomical Society 327 (2001) 1297-1306.
- [19] N. Kaiser et al., Pan-STARRS: A Large Synoptic Survey Telescope Array , Proc. SPIE Int. Soc. Opt. Eng. 4836 (2002) 154.
- [20] DES collaboration, The Dark Energy Survey: more than dark energy-an overview , Monthly Notices of the Royal Astronomical Society 460 (2016) 1270.
- [21] S.R. Furlanetto, S.P. Oh and F.H. Briggs, Cosmology at low frequencies: The 21 cm transition and the high-redshift Universe , Physics Reports 433 (2006) 181.
- [22] T.-C. Chang, U.-L. Pen, J.B. Peterson and P. McDonald, An Intensity Map of Hydrogen 21cm Emission at Redshift z ≈ 0 8, . Nature 466 (2010) 463.
- [23] E. Switzer, K. Masui, K. Bandura, L.-M. Calin et al., Determination of z ∼ 0 8 . neutral hydrogen fluctuations using the 21 cm intensity mapping autocorrelation , Monthly Notices of the Royal Astronomical Society: Letters 434 (2013) L46.
- [24] K. Masui, E. Switzer, N. Banavar, K. Bandura et al., Measurement of 21 cm brightness fluctuations at z ≈ 0 8 . in cross-correlation , The Astrophysical Journal Letters 763 (2013) L20.
- [25] J.R. Pritchard and A. Loeb, 21 cm cosmology in the 21st century , Reports on Progress in Physics 75 (2012) 086901.
- [26] P. Bull, P.G. Ferreira, P. Patel and M.G. Santos, Late-time cosmology with 21 cm intensity mapping experiments , The Astrophysical Journal 803 (2015) 21.
- [27] J.-F. Zhang, L.-Y. Gao, D.-Z. He and X. Zhang, Improving cosmological parameter estimation with the future 21 cm observation from SKA , Physics Letters B 799 (2019) 135064.
- [28] J.-F. Zhang, B. Wang and X. Zhang, Forecast for weighing neutrinos in cosmology with SKA , Science China Physics, Mechanics &amp; Astronomy 63 (2020) 280411.
- [29] S.-J. Jin, D.-Z. He, Y. Xu, J.-F. Zhang and X. Zhang, Forecast for cosmological parameter estimation with gravitational-wave standard siren observation from the Cosmic Explorer , Journal of Cosmology and Astroparticle Physics 03 (2020) 051.
- [30] W. Hu, X. Wang, F. Wu, Y. Wang, P. Zhang and X. Chen, Forecast for FAST: from galaxies survey to intensity mapping , Monthly Notices of the Royal Astronomical Society 493 (2020) 5854-5870.
- [31] M. Zhang, B. Wang, P.-J. Wu, J.-Z. Qi, Y. Xu, J.-F. Zhang et al., Prospects for Constraining Interacting Dark Energy Models with 21 cm Intensity Mapping Experiments , The Astrophysical Journal 918 (2021) 56.
- [32] S.-J. Jin, L.-F. Wang, P.-J. Wu, J.-F. Zhang and X. Zhang, How can gravitational-wave standard sirens and 21-cm intensity mapping jointly provide a precise late-universe cosmological probe? , Physical Review D 104 (2021) 103507.
- [33] P.-J. Wu and X. Zhang, Prospects for measuring dark energy with 21 cm intensity mapping experiments , Journal of Cosmology and Astroparticle Physics 2022 (2022) 060.
- [34] P.-J. Wu, Y. Shao, S.-J. Jin and X. Zhang, A path to precision cosmology: synergy between four promising late-universe cosmological probes , Journal of Cosmology and Astroparticle Physics 06 (2023) 052.
- [35] P.-J. Wu, Y. Li, J.-F. Zhang and X. Zhang, Prospects for measuring dark energy with 21 cm intensity mapping experiments: A joint survey strategy , Science China Physics, Mechanics &amp; Astronomy 66 (2023) 270413.
- [36] X. Chen, A joint survey strategy helps 21 cm intensity mapping become a powerful cosmological probe , Science China Physics, Mechanics &amp; Astronomy 66 (2023) 270431.

- [37] M. Zhang, Y. Li, J.-F. Zhang and X. Zhang, Hi intensity mapping with MeerKAT: forecast for delay power spectrum measurement using interferometer mode , Monthly Notices of the Royal Astronomical Society 524 (2023) 2420.
- [38] Y. Li, Y. Wang, F. Deng et al., FAST Drift Scan Survey for Hi Intensity Mapping: I. Preliminary Data Analysis , The Astrophysical Journal 954 (2023) 139.
- [39] S. Cunnington, Y. Li, M.G. Santos, J. Wang et al., Hi intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies , Monthly Notices of the Royal Astronomical Society 518 (2022) 6262-6272.
- [40] J.L. Bernal, P.C. Breysse and E.D. Kovetz, Cosmic Expansion History from Line-Intensity Mapping , Physical Review Letters 123 (2019) 251301.
- [41] J.L. Bernal, P.C. Breysse, H. Gil-Mar´ ın and E.D. Kovetz, User's guide to extracting cosmological information from line-intensity maps , Physical Review D 100 (2019) 123522.
- [42] H. Fronenberg, A.S. Maniyar, A. Liu and A.R. Pullen, New Probe of the High-z Baryon Acoustic Oscillation Scale: BAO Tomography with CMB × LIM-Nulling Convergence , Physical Review Letters 132 (2024) 241001.
- [43] R. Nan, D. Li, C. Jin et al., The Five-Hundred-Meter Aperture Spherical Radio Telescope (FAST) Project , International Journal of Modern Physics D 20 (2011) 989.
- [44] M. Santos, P. Bull, S. Camera, S. Chen et al., A Large Sky Survey with MeerKAT , Proceedings of Science MeerKAT2016 (2018) 032.
- [45] J. Wang, M.G. Santos, P. Bull, K. Grainge et al., H I intensity mapping with MeerKAT: calibration pipeline for multidish autocorrelation observations , Monthly Notices of the Royal Astronomical Society 505 (2021) 3698.
- [46] Y. Li, M.G. Santos, K. Grainge, S. Harper and J. Wang, Hi intensity mapping with MeerKAT: 1/f noise analysis , Monthly Notices of the Royal Astronomical Society 501 (2020) 4344.
- [47] X. Chen, Radio detection of dark energy-the Tianlai project , Scientia Sinica Physica, Mechanica &amp; Astronomica 41 (2011) 1358.
- [48] X. Chen, The Tianlai Project: a 21CM Cosmology Experiment , International Journal of Modern Physics: Conference Series 12 (2012) 256.
- [49] J. Li, S. Zuo, F. Wu, Y. Wang et al., The Tianlai Cylinder Pathfinder array: System functions and basic performance analysis , Science China Physics, Mechanics &amp; Astronomy 63 (2020) 129862.
- [50] F. Wu, J. Li, S. Zuo, X. Chen et al., The Tianlai dish pathfinder array: design, operation, and performance of a prototype transit radio interferometer , Monthly Notices of the Royal Astronomical Society 506 (2021) 3455.
- [51] O. Perdereau, R. Ansari, A. Stebbins, P.T. Timbie, X. Chen et al., The Tianlai dish array low-z surveys forecasts , Monthly Notices of the Royal Astronomical Society 517 (2022) 4637.
- [52] S. Sun, J. Li, F. Wu et al., The Electromagnetic Characteristics of the Tianlai Cylindrical Pathfinder Array , Research in Astronomy and Astrophysics 22 (2022) 065020.
- [53] L.B. Newburgh et al., Calibrating CHIME, A New Radio Interferometer to Probe Dark Energy , Proceedings of the SPIE International Society for Optical Engineering 9145 (2014) 4V.
- [54] C.J. Anderson et al., Low-amplitude clustering in low-redshift 21-cm intensity maps cross-correlated with 2dF galaxy densities , Monthly Notices of the Royal Astronomical Society 476 (2018) 3382.
- [55] D. Tramonte, Y.-Z. Ma et al., The neutral hydrogen distribution in large-scale haloes from 21-cm intensity maps , Monthly Notices of the Royal Astronomical Society 498 (2020) 5916.

- [56] L. Wolz, A. Pourtsidou et al., HI constraints from the cross-correlation of eBOSS galaxies and Green Bank Telescope intensity maps , Monthly Notices of the Royal Astronomical Society 510 (2021) 3495.
- [57] M. Amiri, K. Bandura, T. Chen, M. Deng et al., Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment , The Astrophysical Journal 947 (2023) 16.
- [58] S. Paul, M.G. Santos, Z. Chen and L. Wolz, A first detection of neutral hydrogen intensity mapping on Mpc scales at z ≈ 0 32 . and z ≈ 0 44, . 2301.11943 .
- [59] P.E. Dewdney, P.J. Hall et al., The Square Kilometre Array , Proceedings of the IEEE 97 (2009) 1482.
- [60] M.G. Santos, P. Bull, D. Alonso et al., Cosmology with a SKA HI intensity mapping survey , 1501.03989 .
- [61] R. Braun, T. Bourke, J.A. Green, E. Keane and J. Wagg, Advancing Astrophysics with the Square Kilometre Array , Proceedings of Science AASKA14 (2015) 174.
- [62] D.J. Bacon, R.A. Battye, P. Bull, S. Camera et al., Cosmology with Phase 1 of the Square Kilometre Array Red Book 2018: Technical specifications and performance forecasts , Publications of the Astronomical Society of Australia 37 (2020) e007.
- [63] T. An, X. Wu, B. Lao, S. Guo et al., Status and progress of China SKA Regional Centre prototype , Science China Physics, Mechanics &amp; Astronomy 65 (2022) 129501.
- [64] R.A. Battye, M.L. Brown, I.W.A. Browne et al., BINGO: a single dish approach to 21cm intensity mapping , 1209.1041 .
- [65] C. Wuensche, The BINGO telescope: a new instrument exploring the new 21 cm cosmology window , Journal of Physics: Conference Series 1269 (2019) 012002.
- [66] L.B. Newburgh et al., HIRAX: A Probe of Dark Energy and Radio Transients , Proceedings of the SPIE International Society for Optical Engineering 9906 (2016) 99065X.
- [67] K. Vanderlinde, A. Liu, B. Gaensler, D. Bond et al., The Canadian Hydrogen Observatory and Radio-transient Detector (CHORD) , 1911.01777 .
- [68] D. Li, P. Wang, W.W. Zhu et al., A bimodal burst energy distribution of a repeating fast radio burst source , Nature 598 (2021) 267-271.
- [69] C.-H. Niu, K. Aggarwal, D. Li, X. Zhang et al., A repeating fast radio burst associated with a persistent radio source , Nature 606 (2022) 873-877.
- [70] H. Xu, S. Chen, Y. Guo et al., Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I , Research in Astronomy and Astrophysics 23 (2023) 075024.
- [71] C.-P. Zhang, P. Jiang, M. Zhu, J. Pan, C. Cheng et al., Performance of FAST with an Ultra-Wide Bandwidth Receiver at 500-3300 MHz , Research in Astronomy and Astrophysics 23 (2023) 075016.
- [72] M. Xue, W. Zhu, X. Wu, R. Xu and H. Wang, Pulsar Discovery Prospect of FASTA , Research in Astronomy and Astrophysics 23 (2023) 095005.
- [73] J. Yin, P. Jiang and R. Yao, Pose optimization of the FAST feed support system based on the new feed cabin mechanism , Science China Physics, Mechanics &amp; Astronomy 66 (2023) 239513.
- [74] P. Jiang, R. Chen, H. Gan, J. Sun, B. Zhu and othes, The fast core array , Astronomical Techniques and Instruments 1 (2024) 84.
- [75] J. Landon, M. Elmer, J. Waldron, D. Jones et al., Phased array feed calibration, beamforming, and imaging , The Astronomical Journal 139 (2010) 1154.

- [76] A. Hotan, J. Bunton, A. Chippendale et al., Australian square kilometre array pathfinder: I. system description , Publications of the Astronomical Society of Australia 38 (2021) e009.
- [77] Planck collaboration, Planck 2018 results. VI. Cosmological parameters , Astronomy &amp; Astrophysics 641 (2020) A6.
- [78] Y. Xu, X. Wang and X. Chen, Forecasts on the Dark Energy and Primordial Non-Gaussianity Observations with the Tianlai Cylinder Array , The Astrophysical Journal 798 (2015) 40.
- [79] C. Li, Y.P. Jing, G. Kauffmann, G. Borner, X. Kang and L. Wang, Luminosity dependence of the spatial and velocity distributions of galaxies: semi-analytic models versus the Sloan Digital Sky Survey , Monthly Notices of the Royal Astronomical Society 376 (2007) 984-996.
- [80] A. Lewis, A. Challinor and A. Lasenby, Efficient Computation of Cosmic Microwave Background Anisotropies in Closed Friedmann-Robertson-Walker Models , The Astrophysical Journal 538 (2000) 473-476.
- [81] X. Wang, M. Tegmark, M.G. Santos and L. Knox, 21 cm Tomography with Foregrounds , The Astrophysical Journal 650 (2006) 529-537.
- [82] M.F. Morales, B. Hazelton, I. Sullivan and A. Beardsley, Four Fundamental Foreground Power Spectrum Shapes for 21 cm Cosmology Observations , The Astrophysical Journal 752 (2012) 137.
- [83] A.R. Parsons, J.C. Pober, J.E. Aguirre, C.L. Carilli, D.C. Jacobs and D.F. Moore, A Per-Baseline, Delay-Spectrum Technique for Accessing the 21cm Cosmic Reionization Signature , The Astrophysical Journal 756 (2012) 165.
- [84] A. Liu, A.R. Parsons and C.M. Trott, Epoch of reionization window. I. Mathematical formalism , Physical Review D 90 (2014) 023018.
- [85] J.R. Shaw, K. Sigurdson, M. Sitwell, A. Stebbins and U.-L. Pen, Coaxing cosmic 21 cm fluctuations from the polarized sky using m-mode analysis , Physical Review D 91 (2015) 083514.
- [86] H.-M. Zhu, U.-L. Pen, Y. Yu and X. Chen, Recovering lost 21cm radial modes via cosmic tidal reconstruction , Physical Review D 98 (2018) 043511.
- [87] S. Zuo, X. Chen, R. Ansari and Y. Lu, 21 cm Signal Recovery via Robust Principal Component Analysis , The Astronomical Journal 157 (2018) 4.
- [88] I.P. Carucci, M.O. Irfan and J. Bobin, Recovery of 21-cm intensity maps with sparse component separation , Monthly Notices of the Royal Astronomical Society 499 (2020) 304-319.
- [89] S. Cunnington, M.O. Irfan, I.P. Carucci, A. Pourtsidou and J. Bobin, 21-cm foregrounds and polarization leakage: cleaning and mitigation strategies , Monthly Notices of the Royal Astronomical Society 504 (2021) 208-227.
- [90] S. Ni, Y. Li, L.-Y. Gao and X. Zhang, Eliminating Primary Beam Effect in Foreground Subtraction of Neutral Hydrogen Intensity Mapping Survey with Deep Learning , The Astrophysical Journal 934 (2022) 83.
- [91] L.-Y. Gao, Y. Li, S. Ni and X. Zhang, Eliminating polarization leakage effect for neutral hydrogen intensity mapping with deep learning , Monthly Notices of the Royal Astronomical Society 525 (2023) 5278.
- [92] M. Spinelli, I.P. Carucci, S. Cunnington, S.E. Harper, M.O. Irfan, J. Fonseca et al., SKAO H I intensity mapping: blind foreground subtraction challenge , Monthly Notices of the Royal Astronomical Society 509 (2022) 2048.
- [93] M.G. Santos, A. Cooray and L. Knox, Multifrequency Analysis of 21 Centimeter Fluctuations from the Era of Reionization , The Astrophysical Journal 625 (2005) 575-587.

- [94] P. Jiang, N.-Y. Tang et al., The fundamental performance of FAST with 19-beam receiver at L band , Research in Astronomy and Astrophysics 20 (2020) 064.
- [95] R. Braun, A. Bonaldi, T. Bourke, E. Keane and J. Wagg, Anticipated Performance of the Square Kilometre Array - Phase 1 (SKA1) , 1912.12699 .
- [96] A. Witzemann, P. Bull, C. Clarkson, M.G. Santos, M. Spinelli and A. Weltman, Model-independent curvature determination with 21cm intensity mapping experiments , Monthly Notices of the Royal Astronomical Society: Letters 477 (2018) L122-L127.
- [97] N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi et al., Planck 2018 results: V. CMB power spectra and likelihoods , Astronomy &amp; Astrophysics 641 (2020) A5.
- [98] G. Efstathiou and S. Gratton, A Detailed Description of the CAMSPEC Likelihood Pipeline and a Reanalysis of the Planck High Frequency Maps , The Open Journal of Astrophysics 4 (2021) 8.
- [99] J. Carron, M. Mirmelstein and A. Lewis, CMB lensing from Planck PR4maps , Journal of Cosmology and Astroparticle Physics 2022 (2022) 039.
- [100] J. Carron, A. Lewis and G. Fabbian, Planck integrated Sachs-Wolfe-lensing likelihood and the CMB temperature , Physical Review D 106 (2022) 103507.
- [101] E. Rosenberg, S. Gratton and G. Efstathiou, CMB power spectra and cosmological parameters from Planck PR4 with CamSpec , Monthly Notices of the Royal Astronomical Society 517 (2022) 4620-4636.
- [102] D. Scolnic, D. Brout, A. Carr, A.G. Riess et al., The Pantheon+ Analysis: The Full Data Set and Light-curve Release , The Astrophysical Journal 938 (2022) 113.
Jun-Da Pan Peng-Ju Wu Guo-Hong Du Yichao Li Xin Zhang
Prospects for cosmological research with the FAST array: 21-cm intensity mapping survey observation strategies
Precise cosmological measurements are essential for understanding the
evolution of the universe and the nature of dark energy. The Five-hundred-meter
Aperture Spherical Telescope (FAST), the most sensitive single-dish radio
telescope, has the potential to provide the precise cosmological measurements
through neutral hydrogen 21-cm intensity mapping sky survey. This paper
primarily explores the potential of technological upgrades for FAST in
cosmology. The most crucial upgrade begins with equipping FAST with a wide-band
receiver ($0 < z < 2.5$). This upgrade can enable FAST to achieve higher
precision in cosmological parameter estimation than the Square Kilometre Array
Phase-1 Mid frequency. On this basis, expanding to a FAST array (FASTA)
consisting of six identical FASTs would offer significant improvements in
precision compared to FAST. Additionally, compared with the current results
from the data combination of cosmic microwave background, baryon acoustic
oscillations (optical galaxy surveys), and type Ia supernovae, FASTA can
provide comparable constraints. Specifically, for the dark-energy
equation-of-state parameters, FASTA can achieve $\sigma(w_0) = 0.09$ and
$\sigma(w_a) = 0.33$.