Forecast measurement of the 21 cm global spectrum from Lunar orbit with the Vari-Zeroth-Order Polynomial (VZOP) method

Tianyang Liu,Jiajun Zhang,Yuan Shi,Junhua Gu,Quan Guo,Yidong Xu,Furen Deng,Fengquan Wu,Yanping Cong,Xuelei Chen
2024-06-24
Abstract:The cosmic 21 cm signal serves as a crucial probe for studying the evolutionary history of the Universe. However, detecting the 21 cm signal poses significant challenges due to its extremely faint nature. To mitigate the interference from the Earth's radio frequency interference (RFI), the ground and the ionospheric effects, the Discovering the Sky at the Longest Wavelength (DSL) project will deploy a constellation of satellites in Lunar orbit, with its high-frequency daughter satellite tasked with detecting the global 21 cm signal from cosmic dawn and reionization era (CD/EoR). We intend to employ the Vari-Zeroth-Order Polynomial (VZOP) for foreground fitting and subtracting. We have studied the effect of thermal noise, thermal radiation from the Moon, the Lunar reflection, anisotropic frequency-dependent beam, inaccurate antenna beam pattern, and RFI contamination. We discovered that the RFI contamination can significantly affect the fitting process and thus prevent us from detecting the signal. Therefore, experimenting on the far side of the moon is crucial. We also discovered that using VZOP together with DSL, after 1080 orbits around the Moon, which takes about 103 days, we can successfully detect the CD/EoR 21 cm signal.
Instrumentation and Methods for Astrophysics,Cosmology and Nongalactic Astrophysics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to successfully detect the global 21 - centimeter signal from Cosmic Dawn (CD) and the Epoch of Reionization (EoR) using the Vari - Zeroth - Order Polynomial (VZOP) method in lunar orbit.** Specifically, the main challenges faced by the researchers include: 1. **The 21 - centimeter signal is extremely weak**: The intensity of the 21 - centimeter signal is expected to be less than 200 mK, while the foreground temperature generated by the synchrotron radiation in the Milky Way is about 4 - 5 orders of magnitude higher than the signal. 2. **Interference on Earth**: The ionosphere on Earth will strongly refract, reflect and absorb low - frequency electromagnetic waves, resulting in an antenna temperature deviation that is 2 - 3 orders of magnitude larger than the global 21 - centimeter signal. In addition, there is a large amount of natural and artificial radio - frequency interference (RFI) on Earth, and these interferences will also seriously affect the observation. 3. **Instrument noise and other errors**: Antenna chromaticity effects, ground reflections, thermal noise, etc. will all introduce additional errors, making signal detection more difficult. To overcome these challenges, the researchers proposed the following solutions: 1. **Select the far side of the moon as the observation site**: The far side of the moon has no atmosphere and ionosphere, and can effectively shield the electromagnetic radiation from the earth and the sun, providing a relatively clean radio environment. 2. **Use the VZOP method for foreground fitting**: The VZOP method can better remove the additional spectral structure introduced by the antenna chromaticity effect by incorporating the antenna beam pattern information into the polynomial fitting model, so as to recover the 21 - centimeter signal more accurately. Through simulation experiments, the researchers verified the effectiveness of the VZOP method and found that after 1,080 lunar orbits (about 103 days), the 21 - centimeter signal from the CD/EoR period can be successfully detected. ### Key formulas 1. **Antenna temperature integration formula**: \[ \bar{T}_A(\nu) = \frac{1}{n_0 p} \sum_{n = 1}^{n_0} \int_{t_n - p/2}^{t_n + p/2} \left[ \int B(\nu, n, t) S(n) T_b(\nu, n, t) d\Omega \right] dt \] where \( T_b(\nu, n, t) \) represents the true sky temperature at frequency \(\nu\), direction \(n\) and time \(t\), \( S(n) \) is the occlusion function, \( B(\nu, n, t) \) is the antenna gain, and \(\Omega\) is the solid angle. 2. **Average antenna temperature formula**: \[ \bar{T}_A(\nu) \approx \sum_{n = 1}^{n_0} \frac{1}{2\pi n_0} \int S(\theta) \bar{T}_b(\nu, \theta, t_n) \bar{B}(\nu, \theta, t_n) \cos\theta d\theta \] where \(\bar{B}(\nu, \theta, t_n)\) is the beam model averaged along the azimuth, and \(\bar{T}_b(\nu, \theta, t_n)\) is the sky temperature model averaged along the azimuth. 3. **Gaussian model representing the 21 - centimeter signal**: \[ T_{\text{eor}}(\nu) = A \exp \left( -\frac{(\nu - \nu_c)^2}{2\omega^2} \right) \] where the amplitude \( A=- 0.150 \, \text{K}\), the center frequency \(\nu_c = 78.3 \, \text{MHz}\), and the width \(\omega = 5 \, \text{MHz}\). 4. **Thermal noise formula**: \[ \sigma_n(\nu)=\frac{