Search

We present the first results of the HI intensity mapping power spectrum analysis with the MeerKAT International GigaHertz Tiered Extragalactic Exploration (MIGHTEE) survey. We use data covering $\sim$ 4 square degrees in the COSMOS field using a frequency range 962.5 MHz to 1008.42 MHz, equivalent to HI emission in $0.4

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

Removing contaminants is a delicate yet crucial step in neutral hydrogen (HI) intensity mapping, often considered the technique's greatest challenge. Here, we address this challenge by analysing HI intensity maps of about $100$ deg$^2$ at redshift $z\approx0.4$ collected by the MeerKAT radio telescope, a SKA Observatory (SKAO) precursor, with a combined 10.5-hour observation. Using unsupervised statistical methods, we remove the contaminating foreground emission and systematically test step-by-step common pre-processing choices to facilitate the cleaning process. We also introduce and test a novel multiscale approach, where data is redundantly decomposed into subsets referring to different spatial scales (large and small), and the cleaning procedure is performed independently. We confirm the detection of the HI cosmological signal in cross-correlation with an ancillary galactic data set without the need to correct for signal loss. In the best set-up reached, we constrain the HI distribution through the combination of its cosmic abundance ($\Omega_{\rm HI}$) and linear clustering bias ($b_{\rm HI}$) up to a cross-correlation coefficient ($r$) and measure $\Omega_{\rm HI}b_{\rm HI}r = [0.93 \pm 0.17]\,\times\,10^{-3}$ with $\approx6\sigma$ confidence. The measurement is independent of scale cuts at both edges of the probed scale range ($0.04 \lesssim k \lesssim 0.3 \,h$Mpc$^{-1}$), corroborating its robustness. Our new pipeline has successfully found an optimal compromise in separating contaminants without incurring a catastrophic signal loss, instilling more confidence in the outstanding science we can deliver with MeerKAT on the path towards HI intensity mapping surveys with the full SKAO., Comment: 24 pages (+ appendix), 24 figures. Summary of results in Table 1 and Figure 17 on page 16. Comments welcome

We present results from MeerKAT single-dish HI intensity maps, the final observations to be performed in L-band in the MeerKAT Large Area Synoptic Survey (MeerKLASS) campaign. The observations represent the deepest single-dish HI intensity maps to date, produced from 41 repeated scans over $236\,{\rm deg}^2$, providing 62 hours of observational data for each of the 64 dishes before flagging. By introducing an iterative self-calibration process, the estimated thermal noise of the reconstructed maps is limited to ${\sim}\,1.21\,$mK ($1.2\,\times$ the theoretical noise level). This thermal noise will be sub-dominant relative to the HI fluctuations on large scales ($k\,{\lesssim}\,0.15\,h\,\text{Mpc}^{-1}$), which demands upgrades to power spectrum analysis techniques, particularly for covariance estimation. In this work, we present the improved MeerKLASS analysis pipeline, validating it on both a suite of mock simulations and a small sample of overlapping spectroscopic galaxies from the Galaxy And Mass Assembly (GAMA) survey. Despite only overlapping with ${\sim}\,25\%$ of the MeerKLASS deep field, and a conservative approach to covariance estimation, we still obtain a ${>}\,4\,\sigma$ detection of the cross-power spectrum between the intensity maps and the 2269 galaxies at the narrow redshift range $0.39\,{<}\,z\,{<}\,0.46$. We briefly discuss the HI auto-power spectrum from this data, the detection of which will be the focus of follow-up work. For the first time with MeerKAT single-dish intensity maps, we also present evidence of HI emission from stacking the maps onto the positions of the GAMA galaxies., Comment: 30 pages, 27 figures. See Fig.14 and 15 for auto-HI and HI-galaxy cross power spectra. Fig 19 and 20 for stacking detections. Accepted to MNRAS

We present a simulation-based framework to forecast the HI power spectrum on non-linear scales ($k\gtrsim 1\ {\rm Mpc^{-1}}$), as measured by interferometer arrays like MeerKAT in the low-redshift ($z\leq 1.0$) universe. Building on a galaxy-based HI mock catalog, we meticulously consider various factors, including the emission line profiles of HI discs and some observational settings, and explore their impacts on the HI power spectrum. While it is relatively insensitive to the profile shape of HI emission line at these scales, we identify a strong correlation with the profile width, that is, the Full Width at Half Maxima (FWHM, also known as $W_{\rm 50}$ in observations) in this work. By modeling the width function of $W_{50}$ as a function of $v_{\rm max}$, we assign each HI source a emission line profile and find that the resulting HI power spectrum is comparatively close to results from particles in the IllustrisTNG hydrodynamical simulation. After implementing $k$-space cuts matching the MeerKAT data, our prediction replicates the trend of the measurements obtained by MeerKAT at $z\approx 0.44$, though with a significantly lower amplitude. Utilizing a Monte Carlo Markov Chain sampling method, we constrain the parameter $A_{W_{\rm 50}}$ in the $W_{\rm 50}$ models and $\Omega_{\rm HI}$ with the MeerKAT measurements and find that a strong degeneracy exists between these two parameters.

Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital sources in the Radio Navigational Satellite System within the 1100-1350 MHz frequency. This simulation can be split into two parts: (I) satellite positioning, emission power, and beam response on the telescope and (II) fitting of the satellite signal to data in order to improve the original model. We use previously observed single dish MeerKAT L-band data which needs to be specially calibrated to include data contaminated by satellite-based RFI. We find that due to non-linearity effects, it becomes non-trivial to fit the satellite power. However, when masking regions where this non-linearity is problematic, we can recreate the satellite contamination with high accuracy around its peak frequencies. The simulation can predict satellite movements and signal for past and future observations, which can help in RFI avoidance and testing novel cleaning methods. The predicted signal from simulations sits below the noise in the target cosmology window for the L-band (970 - 1015 MHz) making it difficult to confirm any out-of-band emission from satellites. However, a power spectrum analysis shows that such signal can still contaminate the 21cm power spectrum at these frequencies. In our simulations, this contamination overwhelms the auto-power spectrum but still allows for a clean detection of the signal in cross-correlations with mild foreground cleaning. Whether such contamination does exist one will require further characterization of the satellite signals far away from their peak frequencies., Comment: 20 pages, 27 figures, 4 tables, Comments welcome

Estimators for $n$-point clustering statistics in Fourier-space demand that modern surveys of large-scale structure be transformed to Cartesian coordinates to perform Fast Fourier Transforms (FFTs). In this work, we explore this transformation in the context of pixelised line intensity maps (LIM), highlighting potential biasing effects on power spectrum measurements. Current analyses often avoid a complete resampling of the data by approximating survey geometry as rectangular in Cartesian space, an increasingly inaccurate assumption for modern wide-sky surveys. Our simulations of a $20\,{\times}\,20\,\text{deg}^2$ 21cm LIM survey at $0.34\,{<}\,z\,{<}\,0.54$ show this assumption biases power spectrum measurements by ${>}\,20\%$ across all scales. We therefore present a more robust framework for regridding the voxel intensities onto a 3D FFT field by coordinate transforming large numbers of Monte-Carlo sampling particles. Whilst this unbiases power spectrum measurements on large scales, smaller-scale discrepancies remain, caused by structure smoothing and aliasing from separations unresolved by the grid. To correct these effects, we introduce modelling techniques, higher-order particle assignments, and interlaced FFT grids to suppress the aliased power. Using a Piecewise Cubic Spline (PCS) particle assignment and an interlaced FFT field, we achieve sub-percent accuracy up to 80% of the Nyquist frequency for our 21cm LIM simulations. We find a more subtle hierarchical improvement in results for higher-order assignment schemes, relative to the gains made for galaxy surveys, which we attribute to the extra complexity in LIM from additional discretising steps. Python code accompanying this paper is available at github.com/stevecunnington/gridimp., Comment: 15 pages, 11 figures. See Fig.9 for headline result showing sub-percent accurate power spectrum estimation. 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

We present a survey strategy to detect the neutral hydrogen (HI) power spectrum at $5

Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensity mapping experiments. However, assessing whether this approach is viable for future intensity mapping surveys where precision cosmology is the aim, remains unexplored. In this work, using simulations, we validate for the first time the use of a foreground transfer function to reconstruct power spectra of foreground-cleaned low-redshift intensity maps and look to expose any limitations. We reveal that even when aggressive foreground cleaning is required, which causes ${>}\,50\%$ negative bias on the largest scales, the power spectrum can be reconstructed using a transfer function to within sub-percent accuracy. We specifically outline the recipe for constructing an unbiased transfer function, highlighting the pitfalls if one deviates from this recipe, and also correctly identify how a transfer function should be applied in an auto-correlation power spectrum. We validate a method that utilises the transfer function variance for error estimation in foreground-cleaned power spectra. Finally, we demonstrate how incorrect fiducial parameter assumptions (up to ${\pm}100\%$ bias) in the generation of mocks, used in the construction of the transfer function, do not significantly bias signal reconstruction or parameter inference (inducing ${<}\,5\%$ bias in recovered values)., Comment: 25 pages, 20 figures. See Figure 4 for the main demonstration of the transfer function's performance for reconstructing signal loss from foreground cleaning. Accepted for publication in MNRAS

We report the first direct detection of the cosmological power spectrum using the intensity signal from 21-cm emission of neutral hydrogen (HI), derived from interferometric observations with the L-band receivers of the new MeerKAT radio telescope. Intensity mapping is a promising technique to map the three-dimensional matter distribution of the Universe at radio frequencies and probe the underlying Cosmology. So far, detections have only been achieved through cross-correlations with galaxy surveys. Here we present independent measurements of the HI power spectrum at redshifts $0.32$ and $0.44$ with high statistical significance using a foreground avoidance method (at $8.0\sigma$ and $11.5\sigma$ respectively). We constrain the rms of the fluctuations of the HI distribution to be $\sigma_{\rm HI} = (0.44\pm 0.04)\,{\rm mK}$ and $\sigma_{\rm HI} = (0.63\pm 0.03)\,{\rm mK}$ respectively at scales of 1.0 Mpc. The information contained in the power spectrum measurements allows us to probe the parameters of the HI mass function and HI halo model. These results are a significant step towards precision cosmology with HI intensity mapping using the new generation of radio telescopes., Comment: 19 pages, 10 figures, submitted to ApJL

We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7\sigma$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($\Omega_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constraint $\Omega_{\rm HI} b_{\rm HI} r\,{=}\,[0.86\,{\pm}\,0.10\,({\rm stat})\,{\pm}\,0.12\,({\rm sys})]\,{\times}\,10^{-3}$, at an effective scale of $k_{\rm eff}\,{\sim}\,0.13\,h\,{\rm Mpc}^{-1}$. The intensity maps were obtained from a pilot survey with the MeerKAT telescope, a 64-dish pathfinder array to the SKA Observatory (SKAO). The data were collected from 10.5 hours of observations using MeerKAT's L-band receivers over six nights covering the 11hr field of WiggleZ, in the frequency range $1015-973\,{\rm MHz}$ (0.400$\,{<}\,z\,{<}\,$0.459 in redshift). This detection is the first practical demonstration of the multi-dish auto-correlation intensity mapping technique for cosmology. This marks an important milestone in the roadmap for the cosmology science case with the full SKAO., Comment: 11 pages, 7 figures. See Fig.5 for the main result, showing the cross-power spectrum between MeerKAT and WiggleZ, revealing the 7.7\sigma detection. v2 is the version accepted for publication by MNRAS. The main changes from v1 were in introducing the radio intensity mapping data

We conduct the first case study towards developing optimal foreground mitigation strategies for neutral hydrogen (HI) intensity mapping using radio interferometers at low redshifts. A pipeline for simulation, foreground mitigation and power spectrum estimation is built, which can be used for ongoing and future surveys using MeerKAT and Square Kilometre Array Observatory (SKAO). It simulates realistic sky signals to generate visibility data given instrument and observation specifications, which is subsequently used to perform foreground mitigation and power spectrum estimation. A quadratic estimator formalism is developed to estimate the temperature power spectrum in visibility space. Using MeerKAT telescope specifications for observations in the redshift range z~0.25-0.30 corresponding to the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey, we present a case study where we compare different approaches of foreground mitigation. We find that component separation in visibility space provides a more accurate estimation of HI clustering comparing to foreground avoidance, with the uncertainties being 30 per cent smaller. Power spectrum estimation from image is found to be less robust with larger bias and more information loss when compared to estimation in visibility. We conclude that for z~0.25-0.30, the MIGHTEE survey will be capable of measuring the HI power spectrum from k~0.5 Mpc$^{-1}$ to k~10 Mpc$^{-1}$ with high accuracy. We are the first to show that, at low redshift, component separation in visibility space suppresses foreground contamination at large line-of-sight scales, allowing measurement of HI power spectrum closer to the foreground wedge, crucial for data analysis towards future detections., Comment: 20 pages, 19 figures. Updated to match the published version in MNRAS