Search

The ultra-long wavelength sky ($\nu\lesssim 30$ MHz) is still largely unexplored, as the electromagnetic wave is heavily absorbed and distorted by the ionosphere on Earth. The far-side of the Moon, either in lunar-orbit or on lunar-surface, is the ideal site for observations in this band, and the upcoming Moon-based interferometers will obtain multi-frequency high-resolution sky maps. Making use of the lunar occultation of the sky and the anisotropy of antenna primary beam response, we propose a novel method to reconstruct the ultra-long wavelength spectral shape in multiple directions in the sky using only one antenna on lunar orbit. We apply the method to one antenna on one of the nine daughter satellites of the proposed Discovering the Sky at Longest wavelength (DSL) project. Using simulated observation data between 1 - 30 MHz from one dipole antenna, we find that the spectra for different regions on the sky can be reconstructed very well and the free-free absorption feature in each region can be derived from the reconstructed spectra. This work demonstrates the feasibility to reconstruct the anisotropic ultra-long wavelength spectra with very limited instrumentation on a lunar-orbit, with mature technologies already in place. It extends the application of such kind of satellite in revealing the distribution of free electrons in the Galactic interstellar medium from the distribution of absorption features in the ultra-long wavelength sky., Comment: 10 pages, 10 figures, 2 tables, submitted to ApJ

The China Space Station Telescope (CSST) is the next-generation Stage~IV survey telescope. It can simultaneously perform multi-band imaging and slitless spectroscopic wide- and deep-field surveys in ten years and an ultra-deep field (UDF) survey in two years, which are suitable for cosmological studies. Here we review several CSST cosmological probes, such as weak gravitational lensing, two-dimensional (2D) and three-dimensional (3D) galaxy clustering, galaxy cluster abundance, cosmic void, Type Ia supernovae (SNe Ia), and baryonic acoustic oscillations (BAO), and explore their capabilities and prospects in discovering new physics and opportunities in cosmology. We find that CSST will measure the matter distribution from small to large scales and the expansion history of the Universe with extremely high accuracy, which can provide percent-level stringent constraints on the property of dark energy and dark matter and precisely test the theories of gravity., Comment: 11 pages, 7 figures, 1 table

This study investigates the estimation of the neutral hydrogen (HI) mass function (HIMF) using a Bayesian stacking approach with simulated data for the Five-hundred-meter Aperture Spherical radio Telescope (FAST) HI intensity mapping (HIIM) drift-scan surveys. Using data from the IllustrisTNG simulation, we construct HI sky cubes at redshift $z\sim0.1$ and the corresponding optical galaxy catalogs, simulating FAST observations under various survey strategies, including pilot, deep-field, and ultradeep-field surveys. The HIMF is measured for distinct galaxy populations -- classified by optical properties into red, blue, and bluer galaxies -- and injected with systematic effects such as observational noise and flux confusion caused by the FAST beam. The results show that Bayesian stacking significantly enhances HIMF measurements. For red and blue galaxies, the HIMF can be well constrained with pilot surveys, while deeper surveys are required for the bluer galaxy population. Our analysis also reveals that sample variance dominates over observational noise, emphasizing the importance of wide-field surveys to improve constraints. Furthermore, flux confusion shifts the HIMF toward higher masses, which we address using a transfer function for correction. Finally, we explore the effects of intrinsic sample incompleteness and propose a framework to quantify its impact. This work lays the groundwork for future \hiMF studies with FAST HIIM, addressing key challenges and enabling robust analyses of HI content across galaxy populations., Comment: 14 pages, 7 figures

Loop I/North Polar Spur (NPS) is the giant arc structure above the Galactic plane observed in radio sky. There is a long-standing debate about its origin. While many people believe it is created by nearby supernova remnants (SNRs), some others consider it is a giant bubble close to the Galactic Center (GC), associated with the Fermi Bubble and eROSITA X-ray bubble. In ultralong wavelength band (wavelength $\gtrsim 10$ m or frequency $\lesssim 30$ MHz), particularly below $\sim 10$ MHz, the free-free absorption of radio signal by the diffuse electrons in interstellar medium (ISM) becomes significant, resulting in sky morphology largely differs from higher frequencies. In this paper, we predict the Loop I/NPS morphology at ultralong wavelength band. We develop emissivity models for the two Loop I/NPS origin scenarios. We find that: For SNR model, the full Loop I/NPS is still a bright arc even at frequency as low as $\sim 1$ MHz. In GC model, while above $\sim 3$ MHz it is still visible; below this frequency it is visible only at Galactic latitude $b\gtrsim 30\degree$. At $b\lesssim 30 \degree$ the Loop I/NPS is invisible due to the absorption by ISM electrons between the GC and the Sun. Upcoming space ultralong wavelentgh projects such as DSL and FARSIDE can potentially distinguish these two scenarios and provide decisive information about the origin of Loop I/NPS., Comment: 11 pages,5 figures

We present the results from calibrating the data of the Commensal Radio Astronomy FAST Survey (CRAFTS) for \HI intensity mapping by the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). Using 70 hours of drift-scan observation with the L-band (1.05-1.45GHz) 19-beam receiver, we obtain the data covering $270\,\rm deg^2$ sky area. We employ both the pulsar backend and the spectrum backend to calibrate the spectral time-ordered-data (TOD) before projecting them onto HEALPix maps. We produce calibrated TOD with frequency resolution of 30 kHz and time resolution of 1 s and the map data-cube with frequency resolution of 30kHz and spatial resolution of $2.95\,\rm arcmin^2$. We carefully examine the pointing errors, noise overflow, RFI contamination and their effect on the data quality. The resulting noise level is $\sim$ 5.7 mJy for the calibrated TOD and 1.6 mJy for the map, which is consistent with the theoretical predictions within 5\% at RFI-free channels. We also validate the data by Principal Components Analysis (PCA) and find most foreground components are concentrated in the first 30 modes. We identify 447 isolated bright continuum sources in our data matching the NRAO-VLA Sky Survey (NVSS) catalog, with relative flux error of 8.3\% for TOD and 11.9\% for the map-level. We also measure the \HI emission of 90 galaxies with redshift $z<0.07$ and compare with \HI-MaNGA spectra from the Green Bank Telescope (GBT), yielding an overall relative error of the \HI integral flux of 16.7\%. Our results confirm the feasibility of conducting cosmological \HI signal detection with CRAFTS., Comment: 30 pages, 30 figures, and 3 tables

Neutral hydrogen (HI) intensity mapping (IM) presents great promise for future cosmological large-scale structure surveys. However, a major challenge for HIIM cosmological studies is to accurately subtract the foreground contamination. An accurate beam model is crucial for improving foreground subtraction accuracy. In this work, we develop a stacking-based beam reconstruction method utilizing the radio continuum point sources within the drift-scan field. Based on the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we employ two series of drift-scan survey data and merge the measurements to construct the beam patterns of the FAST L-band 19 feeds. To model the beams, we utilize the Zernike polynomial (ZP), which effectively captures % the side lobes and asymmetric features of the side lobes. asymmetric features of the main beam and the different side lobes. Due to the symmetric location of the beams, the main features of the beams are closely related to the distance from the center of the feed array, e.g., as the distance increases, side lobes become more pronounced. This modeling pipeline leverages the stable drift-scan data to extract beam patterns exclude the reflector's changing effects, and provides a more accurate measurement beam and a more precise model beam for FAST HIIM cosmology surveys., Comment: 16 pages, 12 figures

We present a catalogue of photometric redshifts for galaxies from DESI Legacy Imaging Surveys, which includes $\sim0.18$ billion sources covering 14,000 ${\rm deg}^2$. The photometric redshifts, along with their uncertainties, are estimated through galaxy images in three optical bands ($g$, $r$ and $z$) from DESI and two near-infrared bands ($W1$ and $W2$) from WISE using a Bayesian Neural Network (BNN). The training of BNN is performed by above images and their corresponding spectroscopic redshifts given in DESI Early Data Release (EDR). Our results show that categorizing galaxies into individual groups based on their inherent characteristics and estimating their photo-$z$s within their group separately can effectively improve the performance. Specifically, the galaxies are categorized into four distinct groups based on DESI's target selection criteria: Bright Galaxy Sample (BGS), Luminous Red Galaxies (LRG), Emission Line Galaxies (ELG) and a group comprising the remaining sources, referred to as NON. As measured by outliers of $|\Delta z| > 0.15 (1 + z_{\rm true})$, accuracy $\sigma_{\rm NMAD}$ and mean uncertainty $\overline{E}$ for BNN, we achieve low outlier percentage, high accuracy and low uncertainty: 0.14%, 0.018 and 0.0212 for BGS and 0.45%, 0.026 and 0.0293 for LRG respectively, surpassing results without categorization. However, the photo-$z$s for ELG cannot be reliably estimated, showing result of $>15\%$, $\sim0.1$ and $\sim0.1$ irrespective of training strategy. On the other hand, NON sources can reach 1.9%, 0.039 and 0.0445 when a magnitude cut of $z<21.3$ is applied. Our findings demonstrate that estimating photo-$z$s directly from galaxy images is significantly potential, and to achieve high-quality photo-$z$ measurement for ongoing and future large-scale imaging survey, it is sensible to implement categorization of sources based on their characteristics., Comment: 18 pages, 11 figures, Accepted for publication in MNRAS

Cosmic shear surveys serve as a powerful tool for mapping the underlying matter density field, including non-visible dark matter. A key challenge in cosmic shear surveys is the accurate reconstruction of lensing convergence ($\kappa$) maps from shear catalogs impacted by survey boundaries and masks, which seminal Kaiser-Squires (KS) method are not designed to handle. To overcome these limitations, we previously proposed the Accurate Kappa Reconstruction Algorithm (AKRA), a prior-free maximum likelihood map-making method. Initially designed for flat sky scenarios with periodic boundary conditions, AKRA has proven successful in recovering high-precision $\kappa$ maps from masked shear catalogs. In this work, we upgrade AKRA to AKRA 2.0 by integrating the tools designed for spherical geometry. This upgrade employs spin-weighted spherical harmonic transforms to reconstruct the convergence field over the full sky. To optimize computational efficiency, we implement a scale-splitting strategy that segregates the analysis into two parts: large-scale analysis on the sphere (referred to as AKRA-sphere) and small-scale analysis on the flat sky (referred to as AKRA-flat); the results from both analyses are then combined to produce final reconstructed $\kappa$ map. We tested AKRA 2.0 using simulated shear catalogs with various masks, demonstrating that the reconstructed $\kappa$ map by AKRA 2.0 maintains high accuracy. For the reconstructed $\kappa$ map in unmasked regions, the reconstructed convergence power spectrum $C_\kappa^{\rm{rec}}$ and the correlation coefficient with the true $\kappa$ map $r_\ell$ achieve accuracies of $(1-C_\ell^{\rm{rec}}/C_\ell^{\rm{true}}) \lesssim 1\%$ and $(1-r_\ell) \lesssim 1\%$, respectively. Our algorithm is capable of straightforwardly handling further issues such as inhomogeneous shape measurement noise, which we will address in subsequent analysis., Comment: 19 pages, 16 figures. Comments are welcome!

The 21 cm emission from neutral hydrogen surveys holds great potential as a valuable method for exploring the large-scale structure of the Universe. In this paper, we forecast for the cross-correlation between the Thermal Sunyaev-Zel'dovich (SZ) fluctuations as probed by the Planck satellite, and fluctuations in the HI brightness temperature as probed by the ground-based Five-hundred-meter Aperture Spherical Telescope (FAST), to trace the connection between galaxy clusters and the HI large-scale structure. Assuming that the measurement is limited by instrumental noise rather than by foreground, we estimate the potential detectability of the cross-correlation signal and their improvement in the measurement of the HI cosmic density, the hydrostatic mass bias parameter, and the universal pressure profile (UPP) parameters. We obtain a constraint on the cosmic neutral hydrogen density parameter significantly to $\sigma(\Omega_{\rm HI}) = 1.0 \times 10^{-6}$. We also find that the average halo masses contributing to the ${{\rm HI}-y}$ cross-power spectrum in the one-halo regime is $\sim 1.5\times 10^{14} M_{\odot}$. Our results also show that the HI-SZ cross-correlation has great potential to probe the distribution of neutral hydrogen (HI) within halos at low redshift., Comment: 18 pages, 4 figures, 3 tables, accepted in the Astrophysical Journal

The Gravitation Wave (GW) signals from a large number of double white dwarfs (DWDs) in the Galaxy are expected to be detected by space GW detectors, e.g., the Laser Interferometer Space Antenna (LISA), Taiji, and Tianqin in the millihertz band. In this paper, we present an alternative method by directly using the time-domain GW signal detected by space GW detectors to constrain the anisotropic structure of the Galaxy. The information of anisotropic distribution of DWDs is naturally encoded in the time-domain GW signal because of the variation of the detectors' directions and consequently the pattern functions due to their annual motion around the sun. The direct use of the time-domain GW signal enables simple calculations, such as utilizing an analytical method to assess the noise arising from the superposition of random phases of DWDs and using appropriate weights to improve the constraints. We investigate the possible constraints on the scale of the Galactic thin disk and bulge that may be obtained from LISA and Taiji by using this method with mock signals obtained from population synthesis models. We further show the different constraining capabilities of the low-frequency signal (foreground) and the high-frequency signal (resolvable-sources) via the Markov Chain Monte Carlo method, and find that the scale height and length of the Galactic thin disk and the scale radius of bulge can be constrained to a fractional accuracy of ~ 30%, 30%, 40% (or 20%, 10%, 40%) by using the low-frequency (or high-frequency) signal detected by LISA or Taiji., Comment: 20 pages, 8 figures, accepted for publication in ApJ

Filaments stand as pivotal structures within the cosmic web. However, direct detection of the cold gas content of the filaments remains challenging due to its inherent low brightness temperature. With the TNG hydrodynamical simulations, we demonstrate the effectiveness of isolating faint filament HI signal from the FAST HI intensity mapping (IM) survey through pairwise stacking of galaxies, which yields an average HI filament signal amplitude of $\sim 0.28\ {\mu{\rm K}}$ at $z\simeq 0.1$. However, our simulations reveal a non-negligible contribution from HI-rich galaxies within or near the filaments. Particularly, the faint galaxies dominantly contribute to the extra filament HI signal. Our simulation also shows that the measurement uncertainty is produced by both thermal noise and background variation caused by brightness leakage from surrounding random galaxies. Given a fixed total observation time, a wide-field HI IM survey, which includes a large number of galaxy pairs, can simultaneously reduce thermal noise to below the filament signal level and minimize background variation to a negligible level. Through the end-to-end simulation, this work demonstrates the critical role of the galaxy pairwise stacking method in future filament HI detection, outlining a road map for filament HI detection in the next-generation HI IM surveys., Comment: 16 pages, 16 figures

We explore the joint weak lensing and galaxy clustering analysis from the photometric survey operated by the China Space Station Telescope (CSST), and study the strength of the cosmological constraints. We employ a high-resolution JiuTian-1G simulation to construct a partial-sky light cone to $z=3$ covering 100 deg$^2$, and obtain the CSST galaxy mock samples based on an improved semi-analytical model. We perform a multi-lens-plane algorithm to generate corresponding synthetic weak lensing maps and catalogs. Then we generate the mock data based on these catalogs considering the instrumental and observational effects of the CSST, and use the Markov Chain Monte Carlo (MCMC) method to perform the constraints. The covariance matrix includes non-Gaussian contributions and super-sample covariance terms, and the systematics from intrinsic alignments, galaxy bias, photometric redshift uncertainties, shear calibration, and non-linear effects are considered in the analysis. We find that, for the joint analysis of the CSST weak lensing and galaxy clustering surveys, the cosmological parameters can be constrained to a few percent or even less than one percent level. This indicates the CSST photometric survey is powerful for exploring the Universe., Comment: 17 pages, 12 figures, 2 tables

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: 12 pages, 5 figures, 2 tables. Accepted for publication in ApJ

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