Emission-line Stacking of 21cm Intensity Maps with MeerKLASS: Inference Pipeline and Application to the L-band Deep-field Data

Zhaoting Chen; Steven Cunnington; Alkistis Pourtsidou; Laura Wolz; Marta Spinelli; José Luis Bernal; Matilde Barberi-Squarotti; Stefano Camera; Isabella P. Carucci; José Fonseca; Keith Grainge; Melis O. Irfan; Mario G. Santos; Jingying Wang

doi:10.3847/1538-4365/add897

paper

Emission-line Stacking of 21cm Intensity Maps with MeerKLASS: Inference Pipeline and Application to the L-band Deep-field Data

arXiv:2504.03908 · doi:10.3847/1538-4365/add897

Abstract

We present a novel analysis of observational systematics through the emission-line stacking of the MeerKLASS L-band deep-field intensity maps, following the detection in arXiv:2407.21626. A stacking signal is obtained by stacking the 21cm intensity map cubelets around the galaxy positions from the GAMA survey at $0.39 ≲ z ≲ 0.46$ . An extensive simulation framework is built to study the viability of the stacking detection, the covariance estimation, and the model inference, which are then applied to the data. The statistical significance of the detection is $8.66 σ$ when averaged into an angular map, and $7.45 σ$ when averaged into a spectrum. The stacked spectrum exhibits an oscillating component of systematics, and we provide evidence that these systematics are a convolutional effect on the map data. The oscillation frequency matches the diffraction from the secondary reflector into the primary beam of the MeerKAT telescope. Bayesian inference can be used to constrain the systematics and the average HI emission of the galaxies. The fitting of the parameters gives a constraint on the systematics frequency $ν_{sys} [MHz] = 17.9 0_{- 4.27}^{+ 6.53}$ . The posterior of the systematics amplitude reaches the wide prior and gives $A_{sys} = 0.5 0_{- 0.33}^{+ 0.33}$ . A tentative measurement of the average HI mass of the sources is achieved at $lo g_{10} [⟨ M_{H I} ⟩ / M_{⊙}] = 9.8 4_{- 0.59}^{+ 0.48}$ , which is an underestimation limited by the narrow redshift bin, the strong degeneracy with the systematics, and the low-density galaxy sample. These shortfalls will be resolved for future MeerKLASS data to enable accurate measurements of the HI density through stacking of intensity maps.

arXiv comment39 pages, 25 figures. Updated to match the version accepted by ApJS

Emission-line Stacking of 21cm Intensity Maps with MeerKLASS: Inference Pipeline and Application to the L-band Deep-field Data

Abstract

References in corpus (91)

Cited by in corpus (5)