Search

Void number counts (VNC) indicates the number of low-density regions in the large-scale structure (LSS) of the Universe, and we propose to use it as an effective cosmological probe. By generating the galaxy mock catalog based on Jiutian simulations and considering the spectroscopic survey strategy and instrumental design of the China Space Station Telescope (CSST), which can reach a magnitude limit $\sim$23 AB mag and spectral resolution $R\gtrsim200$ with a sky coverage 17,500 deg$^2$, we identify voids using the watershed algorithm without any assumption of void shape, and obtain the mock void catalog and data of the VNC in six redshift bins from $z=0.3$ to1.3. We use the Markov Chain Monte Carlo (MCMC) method to constrain the cosmological and VNC parameters. The void linear underdensity threshold $\delta_{\rm v}$ in the theoretical model is set to be a free parameter at a given redshift to fit the VNC data and explore its redshift evolution. We find that, the VNC can correctly derive the cosmological information, and the constraint strength on the cosmological parameters is comparable to that from the void size function (VSF) method, which can reach a few percentage levels in the CSST full spectroscopic survey. This is because that, since the VNC is not sensitive to void shape, the modified theoretical model can match the data better by integrating over void features, and more voids could be included in the VNC analysis by applying simpler selection criteria, which will improve the statistical significance. It indicates that the VNC can be an effective cosmological probe for exploring the LSS., Comment: 8 pages, 5 figures, 2 tables. Accepted for publication in MNRAS

The void power spectrum is related to the clustering of low-density regions in the large-scale structure (LSS) of the Universe, and can be used as an effective cosmological probe to extract the information of the LSS. We generate the galaxy mock catalogs from Jiutian simulation, and identify voids using the watershed algorithm for studying the cosmological constraint strength of the China Space Station Telescope (CSST) spectroscopic survey. The galaxy and void auto power spectra and void-galaxy cross power spectra at $z=0.3$, 0.6, and 0.9 are derived from the mock catalogs. To fit the full power spectra, we propose to use the void average effective radius at a given redshift to simplify the theoretical model, and adopt the Markov Chain Monte Carlo (MCMC) technique to implement the constraints on the cosmological and void parameters. The systematical parameters, such as galaxy and void biases, and noise terms in the power spectra are also included in the fitting process. We find that our theoretical model can correctly extract the cosmological information from the galaxy and void power spectra, which demonstrates its feasibility and effectivity. The joint constraint accuracy of the cosmological parameters can be improved by $\sim20\%$ compared to that from the galaxy power spectrum only. The fitting results of the void density profile and systematical parameters are also well constrained and consistent with the expectation. This indicates that the void clustering measurement can be an effective complement to the galaxy clustering probe, especially for the next generation galaxy surveys., Comment: 11 pages, 5 figures, 2 tables

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., Comment: 15 pages, 17 figures, to be submitted to SCPMA

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

At the Royal Society meeting in 2023, we have mainly presented our lunar orbit array concept called DSL, and also briefly introduced a concept of a lunar surface array, LARAF. As the DSL concept had been presented before, in this article we introduce the LARAF. We propose to build an array in the far side of the Moon, with a master station which handles the data collection and processing, and 20 stations with maximum baseline of 10 km. Each station consists 12 membrane antenna units, and the stations are connected to the master station by power line and optical fiber. The array will make interferometric observation in the 0.1-50 MHz band during the lunar night, powered by regenerated fuel cells (RFCs). The whole array can be carried to the lunar surface with a heavy rocket mission, and deployed with a rover in 8 months. Such an array would be an important step in the long term development of lunar based ultralong wavelength radio astronomy. It has a sufficiently high sensitivity to observe many radio sources in the sky, though still short of the dark age fluctuations. We discuss the possible options in the power supply, data communication, deployment, etc., Comment: final submission version, 30 pages, 16 figures

Void size function (VSF) contains information of the cosmic large-scale structure (LSS), and can be used to derive the properties of dark energy and dark matter. We predict the VSFs measured from the spectroscopic galaxy survey operated by the China Space Station Telescope (CSST), and study the strength of cosmological constraint. We employ a high-resolution Jiutian simulation to get CSST galaxy mock samples based on an improved semi-analytical model. We identify voids from this galaxy catalog using the watershed algorithm without assuming a spherical shape, and estimate the VSFs at different redshift bins from $z=0.5$ to 1.1. We propose a void selection method based on the ellipticity, and assume the void linear underdensity threshold $\delta_{\rm v}$ in the theoretical model is redshift-dependent and set it as a free parameter in each redshift bin. The Markov Chain Monte Carlo (MCMC) method is adopted to implement the constraints on the cosmological and void parameters. We find that the CSST VSF measurement can constrain the cosmological parameters to a few percent level. The best-fit values of $\delta_{\rm v}$ are ranging from $\sim-0.4$ to $-0.1$ as the redshift increases from 0.5 to 1.1, which has a distinct difference from the theoretical calculation with $\delta_{\rm v}\simeq-2.7$ assuming the spherical evolution and using particles as tracer. Our method can provide a good reference for void identification and selection in the VSF analysis of the spectroscopic galaxy surveys., Comment: 10 pages, 7 figures, 3 tables. Accepted for publication in MNRAS

Type Ia supernova (SN Ia) as a standard candle is an ideal tool to measure cosmic distance and expansion history of the Universe. Here we investigate the SN Ia photometric measurement in the China Space Station Telescope Ultra Deep Field (CSST-UDF) survey, and study the constraint power on the cosmological parameters, such as the equation of state of dark energy. The CSST-UDF survey is expected to cover a 9 deg$^2$ sky area in two years with 250 s $\times$ 60 exposures for each band. The magnitude limit can reach $i\simeq26$ AB mag for 5$\sigma$ point source detection with a single exposure. We generate light curve mock data for SNe Ia and different types of core-collapse supernovae (CCSNe). {\tt SNCosmo} is chosen as the framework by utilizing the SALT3 model to simulate SN Ia data. After selecting high-quality data and fitting the light curves, we derive the light curve parameters and identify CCSNe as contamination, resulting in $\sim2200$ SNe with a $\sim7\%$ CCSN contamination rate. We adopt a calibration method similar to Chauvenet's criterion, and apply it to the distance modulus data to further reduce the contamination. We find that this method is effective and can suppress the contamination fraction to $\sim3.5\%$ with 2012 SNe Ia and 73 CCSNe. In the cosmological fitting stage, we did not distinguish between SNe Ia and CCSNe. We find that the constraint accuracies on $\Omega_{\rm M}$, $\Omega_{\Lambda}$ and $w$ are about two times better than the current SN surveys, and it could be further improved by a factor of $\sim$1.4 if including the baryon acoustic oscillation (BAO) data from the CSST spectroscopic wide-field galaxy survey., Comment: 12 pages, 12 figures, 5 tables. Accepted for publication in MNRAS

A large number of binary neutron star (BNS) mergers are expected to be detected by gravitational wave (GW) detectors and the electromagnetic (EM) counterparts (e.g., kilonovae) of a fraction of these mergers may be detected in multi-bands by large area survey telescopes. For a given number of BNS mergers detected by their GW signals, the expected numbers of their EM counterparts that can be detected by a survey with given selection criteria depend on the kilonova properties, including the anisotropy. In this paper, we investigate whether the anisotropy of kilonova radiation and the kilonova model can be constrained statistically by the counting method, i.e., using the numbers of BNS mergers detected via GW and multi-band EM signals. Adopting simple models for the BNS mergers, afterglows, and a simple two (blue and red)-component model for kilonovae, we generate mock samples for GW detected BNS mergers, their associated kilonovae and afterglows detected in multi-bands. By assuming some criteria for searching the EM counterparts, we simulate the observations of these EM counterparts and obtain the EM observed samples in different bands. With the numbers of BNS mergers detected by GW detectors and EM survey telescopes in different bands, we show that the anisotropy of kilonova radiation and the kilonova model can be well constrained by using the Bayesian analysis. Our results suggest that the anisotropy of kilonova radiation may be demographically and globally constrained by simply using the detection numbers of BNS mergers by GW detectors and EM survey telescopes in multi-bands., Comment: 12 pages, 7 figures, MNRAS accepted

Ultralight axion (ULA) can be one of the potential candidates for dark matter. The extremely low mass of the ULA can lead to a de Broglie wavelength the size of galaxies which results in a suppression of the growth of structure on small scales. In this work, we forecast the constraint on the ULA particle mass $m_{\text{a}}$ and relative fraction to dark matter $f_{\text{a}} = \Omega_{\text{a}}/\Omega_{\text{d}}$ for the forthcoming Stage IV space-based optical survey equipment $\it{CSST}$ (China Space Station Telescope). We focus on the $\it{CSST}$ cosmic shear and galaxy clustering photometric surveys, and forecast the measurements of shear, galaxy, and galaxy-galaxy lensing power spectra (i.e. 3$\times$2pt). The effects of neutrino, baryonic feedback, and uncertainties of intrinsic alignment, shear calibration, galaxy bias, and photometric redshift are also included in the analysis. After performing a joint constraint on all the cosmological and systematical parameters based on the simulated data from the theoretical prediction, we obtain a lower limit of the ULA particle mass $\text{log}_{10}(m_{\text{a}}/\text{eV}) \geqslant -22.5$ and an upper limit of the ULA fraction $f_{\text{a}} \leqslant 0.83$ at 95\% confidence level, and $\text{log}_{10}(m_{\text{a}}/\text{eV}) \geqslant -21.9$ with $f_{\text{a}} \leqslant 0.77$ when ignoring the baryonic feedback. We find that the CSST photometric surveys can improve the constraint on the ULA mass by about one order of magnitude, compared to the current constraints using the same kind of observational data., Comment: 13 pages, 10 figures, and 2 tables. Accepted for publication in MNRAS

This work presents the initial results of the drift-scan observation for the neutral hydrogen (HI) intensity mapping survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The data analyzed in this work were collected in night observations from 2019 through 2021. The primary findings are based on 28 hours of drift-scan observation carried out over seven nights in 2021, which covers $60\,{\rm deg}^2$ sky area. Our main findings are: (i) Our calibration strategy can successfully correct both the temporal and bandpass gain variation over the $4$-hour drift-scan observation. (ii) The continuum maps of the surveyed region are made with frequency resolution of $28$ kHz and pixel area of $2.95\,{\rm arcmin}^2$. The pixel noise levels of the continuum maps are slightly higher than the forecast assuming $T_{\rm sys}=20\,{\rm K}$, which are $36.0$ mK (for 10.0 s integration time) at the $1050$--$1150$ MHz band, and $25.9$ mK (for 16.7 s integration time) at the $1323$--$1450$ MHz band, respectively. (iii) The flux-weighted differential number count is consistent with the NRAO-VLA Sky Survey (NVSS) catalog down to the confusion limit $\sim7\,{\rm mJy}/{\rm beam}^{-1}$. (iv) The continuum flux measurements of the sources are consistent with that found in the literature. The difference in the flux measurement of $81$ isolated NVSS sources is about $6.3\%$. Our research offers a systematic analysis for the FAST HI intensity mapping drift-scan survey and serves as a helpful resource for further cosmology and associated galaxies sciences with the FAST drift-scan survey., Comment: 26 pages, 26 figures, and 4 tables

In neutral hydrogen (HI) galaxy survey, a significant challenge is to identify and extract the HI galaxy signal from observational data contaminated by radio frequency interference (RFI). For a drift-scan survey, or more generally a survey of a spatially continuous region, in the time-ordered spectral data, the HI galaxies and RFI all appear as regions which extend an area in the time-frequency waterfall plot, so the extraction of the HI galaxies and RFI from such data can be regarded as an image segmentation problem, and machine learning methods can be applied to solve such problems. In this study, we develop a method to effectively detect and extract signals of HI galaxies based on a Mask R-CNN network combined with the PointRend method. By simulating FAST-observed galaxy signals and potential RFI impacts, we created a realistic data set for the training and testing of our neural network. We compared five different architectures and selected the best-performing one. This architecture successfully performs instance segmentation of HI galaxy signals in the RFI-contaminated time-ordered data (TOD), achieving a precision of 98.64% and a recall of 93.59%., Comment: 17 pages, 9 figures, 1 tables. Accepted for publication in RAA

Cross-correlating the data of neutral hydrogen (HI) 21cm intensity mapping with galaxy surveys is an effective method to extract astrophysical and cosmological information. In this work, we investigate the cross-correlation of MeerKAT single-dish mode HI intensity mapping and China Space Station Telescope (CSST) spectroscopic galaxy surveys. We simulate a survey area of $\sim 300$ $\mathrm{deg}^2$ of MeerKAT and CSST surveys at $z=0.5$ using Multi-Dark N-body simulation. The PCA algorithm is applied to remove the foregrounds of HI intensity mapping, and signal compensation is considered to solve the signal loss problem in the HI-galaxy cross power spectrum caused by the foreground removal process. We find that from CSST galaxy auto and MeerKAT-CSST cross power spectra, the constraint accuracy of the parameter product $\Omega_{\rm HI}b_{\rm HI}r_{{\rm HI},g}$ can reach to $\sim1\%$, which is about one order of magnitude higher than the current results. After performing the full MeerKAT HI intensity mapping survey with 5000 deg$^2$ survey area, the accuracy can be enhanced to $<0.3\%$. This implies that the MeerKAT-CSST cross-correlation can be a powerful tool to probe the cosmic HI property and the evolution of galaxies and the Universe., Comment: 18 pages, 11 figures, 3 tables. Accepted for publication in RAA, revised to match the published version

Line intensity mapping (LIM) is a promising probe to study star formation, the large-scale structure of the Universe, and the epoch of reionization (EoR). Since carbon monoxide (CO) is the second most abundant molecule in the Universe except for molecular hydrogen ${\rm H}_2$, it is suitable as a tracer for LIM surveys. However, just like other LIM surveys, CO intensity mapping also suffers strong foreground contamination that needs to be eliminated for extracting valuable astrophysical and cosmological information. In this work, we take $^{12}$CO($\it J$=1-0) emission line as an example to investigate whether deep learning method can effectively recover the signal by removing the foregrounds. The CO(1-0) intensity maps are generated by N-body simulations considering CO luminosity and halo mass relation, and we discuss two cases with median and low CO signals by comparing different relations. We add foregrounds generated from real observations, including thermal dust, spinning dust, free-free, synchrotron emission and CMB anisotropy. The beam with sidelobe effect is also considered. Our deep learning model is built upon ResUNet, which combines image generation algorithm UNet with the state-of-the-art architecture of deep learning, ResNet. The principal component analysis (PCA) method is employed to preprocess data before feeding it to the ResUNet. We find that, in the case of low instrumental noise, our UNet can efficiently reconstruct the CO signal map with correct line power spectrum by removing the foregrounds and recovering PCA signal loss and beam effects. Our method also can be applied to other intensity mappings like neutral hydrogen 21cm surveys., Comment: 12 pages, 12 figures, 2 tables. Accepted for publication in MNRAS

The Tianlai Dish Pathfinder array is a radio interferometer array consisting of 16 six meter dish antennas. The original digital backend integration time is at the seconds level, designed for HI intensity mapping experiment. A new digital backend with millisecond response is added to enable it to search for fast radio burst (FRB) during its observations. The design and calibration of this backend, and the real time search pipeline for it are described in this paper. It is capable of forming 16 digital beams for each linear polarisation, covering an area of 19.6 square degrees. The search pipeline is capable of searching for, recording and classifying FRBs automatically in real time. In commissioning, we succeeded in capturing the signal pulses from the pulsars PSR B0329+54 and B2021+51., Comment: 16 pages, 14 figures, RAA accepted

The cross-correlation of optical galaxies with the neutral hydrogen (HI) radiation intensity can enhance the signal-to-noise ratio (SNR) of the HI intensity measurement. In this paper, we investigate the cross-correlation of the galaxy samples obtained by the spectroscopic survey of the China Space Station Telescope (CSST) with the HI Intensity mapping (IM) survey of the Five-hundred-meter Aperture Spherical Telescope (FAST). Using the IllusitrisTNG simulation result at redshift $0.2 \sim 0.3$, we generate mock data of the CSST survey and a FAST L-band drift scan survey. The CSST spectroscopic survey can yield a sample of galaxies with a high comoving number density of $10^{-2} (\unit{Mpc}/h)^{-3}$ at $z \sim 0.3$. We cross-correlate the foreground-removed radio intensity with the CSST galaxies, including both the whole sample, and red and blue galaxy sub-samples separately. We find that in all cases the HI and optical galaxies are well correlated. The total HI abundance can be measured with a high precision from this correlation. A relative error of $\sim 0.6\%$ for $\Omega_{\rm HI}$ could be achieved at $z\sim 0.3$ for an overlapping survey area of $10000 \unit{deg}^2$., Comment: 16 pages, 10 figures

A redshifted 21 cm line absorption signature is commonly expected from the cosmic dawn era, when the first stars and galaxies formed. The detailed traits of this signal can provide important insight on the cosmic history. However, high precision measurement of this signal is hampered by the ionosphere refraction and absorption, as well as radio frequency interference (RFI). A space observation can solve the problem of the ionosphere, and the Moon can shield the RFI from the Earth. In this paper, we present simulations of the global spectrum measurement in the 30 -- 120 MHz frequency band on the lunar orbit, from the proposed Discovering the Sky at the Longest wavelength (DSL) project. In particular, we consider how the measured signal varies as the satellite moves along the orbit, take into account the blockage of different parts of the sky by the Moon and the antenna response. We estimate the sensitivity for such a 21 cm global spectrum experiment. An RMS noise level of $\le 0.05$ K is expected at 75 MHz after 10 orbits ($\sim$ 1 day) observation, for a frequency channel width of 0.4 MHz. We also study the influence of a frequency-dependent beam, which may generate complex spectral structures in the spectrum. Estimates of the uncertainties in the foreground and 21 cm model parameters are obtained., Comment: 14 pages, 16 figures, accepted for publication in ApJ