-
HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies
Authors:
Steven Cunnington,
Yichao Li,
Mario G. Santos,
Jingying Wang,
Isabella P. Carucci,
Melis O. Irfan,
Alkistis Pourtsidou,
Marta Spinelli,
Laura Wolz,
Paula S. Soares,
Chris Blake,
Philip Bull,
Brandon Engelbrecht,
José Fonseca,
Keith Grainge,
Yin-Zhe Ma
Abstract:
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σ$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($Ω_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constrai…
▽ More
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σ$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($Ω_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constraint $Ω_{\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.
△ Less
Submitted 25 October, 2022; v1 submitted 3 June, 2022;
originally announced June 2022.
-
Gaussian Process Regression for foreground removal in HI intensity mapping experiments
Authors:
Paula S. Soares,
Catherine A. Watkinson,
Steven Cunnington,
Alkistis Pourtsidou
Abstract:
We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source Python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use G…
▽ More
We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source Python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use GPR as a foreground removal technique in this context, and that it is better suited in some cases to recover the HI power spectrum than Principal Component Analysis (PCA), especially on small scales. GPR is especially good at recovering the radial power spectrum, outperforming PCA when considering the full bandwidth of our data. Both methods are worse at recovering the transverse power spectrum, since they rely on frequency-only covariance information. When halving our data along frequency, we find that GPR performs better in the low frequency range, where foregrounds are brighter. It performs worse than PCA when frequency channels are missing, to emulate RFI flagging. We conclude that GPR is an excellent foreground removal option for the case of single-dish, low redshift HI intensity mapping in the absence of missing frequency channels. Our Python toolkit gpr4im and the data used in this analysis are publicly available on GitHub.
△ Less
Submitted 4 October, 2021; v1 submitted 26 May, 2021;
originally announced May 2021.
-
HI intensity mapping with MeerKAT: Calibration pipeline for multi-dish autocorrelation observations
Authors:
Jingying Wang,
Mario G. Santos,
Philip Bull,
Keith Grainge,
Steven Cunnington,
Jose Fonseca,
Melis O. Irfan,
Yichao Li,
Alkistis Pourtsidou,
Paula S. Soares,
Marta Spinelli,
Gianni Bernardi,
Brandon Engelbrecht
Abstract:
While most purpose-built 21cm intensity mapping experiments are close-packed interferometer arrays, general-purpose dish arrays should also be capable of measuring the cosmological 21cm signal. This can be achieved most efficiently if the array is used as a collection of scanning autocorrelation dishes rather than as an interferometer. As a first step towards demonstrating the feasibility of this…
▽ More
While most purpose-built 21cm intensity mapping experiments are close-packed interferometer arrays, general-purpose dish arrays should also be capable of measuring the cosmological 21cm signal. This can be achieved most efficiently if the array is used as a collection of scanning autocorrelation dishes rather than as an interferometer. As a first step towards demonstrating the feasibility of this observing strategy, we show that we are able to successfully calibrate dual-polarisation autocorrelation data from 64 MeerKAT dishes in the L-band (856-1712 MHz, 4096 channels), with 10.5 hours of data retained from six nights of observing. We describe our calibration pipeline, which is based on multi-level RFI flagging, periodic noise diode injection to stabilise gain drifts and an absolute calibration based on a multi-component sky model. We show that it is sufficiently accurate to recover maps of diffuse celestial emission and point sources over a 10 deg x 30 deg patch of the sky overlapping with the WiggleZ 11hr field. The reconstructed maps have a good level of consistency between per-dish maps and external datasets, with the estimated thermal noise limited to 1.4 x the theoretical noise level (~ 2 mK). The residual maps have rms amplitudes below 0.1 K, corresponding to <1% of the model temperature. The reconstructed Galactic HI intensity map shows excellent agreement with the Effelsberg-Bonn HI Survey, and the flux of the radio galaxy 4C+03.18 is recovered to within 3.6%, which demonstrates that the autocorrelation can be successfully calibrated to give the zero-spacing flux and potentially help in the imaging of MeerKAT interferometric data. Our results provide a positive indication towards the feasibility of using MeerKAT and the future SKA to measure the HI intensity mapping signal and probe cosmology on degree scales and above.
△ Less
Submitted 13 May, 2021; v1 submitted 27 November, 2020;
originally announced November 2020.
-
Power spectrum multipole expansion for HI intensity mapping experiments: unbiased parameter estimation
Authors:
Paula S. Soares,
Steven Cunnington,
Alkistis Pourtsidou,
Chris Blake
Abstract:
We assess the performance of the multipole expansion formalism in the case of single-dish HI intensity mapping, including instrumental and foreground removal effects. This formalism is used to provide MCMC forecasts for a range of HI and cosmological parameters, including redshift space distortions and the Alcock-Paczynski effect. We first determine the range of validity of our power spectrum mode…
▽ More
We assess the performance of the multipole expansion formalism in the case of single-dish HI intensity mapping, including instrumental and foreground removal effects. This formalism is used to provide MCMC forecasts for a range of HI and cosmological parameters, including redshift space distortions and the Alcock-Paczynski effect. We first determine the range of validity of our power spectrum modelling by fitting to simulation data, concentrating on the monopole, quadrupole, and hexadecapole contributions. We then show that foreground subtraction effects can lead to severe biases in the determination of cosmological parameters, in particular the parameters relating to the transverse BAO rescaling, the growth rate and the HI bias ($α_\perp$, $\overline{T}_\text{HI} fσ_8$, and $\overline{T}_\text{HI} b_\text{HI} σ_8$, respectively). We attempt to account for these biases by constructing a 2-parameter foreground modelling prescription, and find that our prescription leads to unbiased parameter estimation at the expense of increasing the estimated uncertainties on cosmological parameters. In addition, we confirm that instrumental and foreground removal effects significantly impact the theoretical covariance matrix, and cause the covariance between different multipoles to become non-negligible. Finally, we show the effect of including higher-order multipoles in our analysis, and how these can be used to investigate the presence of instrumental and systematic effects in HI intensity mapping data.
△ Less
Submitted 5 January, 2021; v1 submitted 27 August, 2020;
originally announced August 2020.
-
Multipole expansion for HI intensity mapping experiments: simulations and modelling
Authors:
Steven Cunnington,
Alkistis Pourtsidou,
Paula S. Soares,
Chris Blake,
David Bacon
Abstract:
We present a framework and an open-source python toolkit to analyse the 2-point statistics of 3D fluctuations in the context of HI intensity maps using the multipole expansion formalism. We include simulations of the cosmological HI signal using N-body and log-normal methods, foregrounds and their removal, as well as instrumental effects. Using these simulations and analytical modelling, we invest…
▽ More
We present a framework and an open-source python toolkit to analyse the 2-point statistics of 3D fluctuations in the context of HI intensity maps using the multipole expansion formalism. We include simulations of the cosmological HI signal using N-body and log-normal methods, foregrounds and their removal, as well as instrumental effects. Using these simulations and analytical modelling, we investigate the impact of foreground cleaning and the instrumental beam on the power spectrum multipoles as well as on the Fourier space clustering wedges. We find that both the instrumental beam and the foreground removal can produce a quadrupole (and a hexadecapole) signal, and demonstrate the importance of controlling and accurately modelling these effects for precision radio cosmology. We conclude that these effects can be modelled with reasonable accuracy using our multipole expansion technique. We also perform an MCMC analysis to showcase the effect of foreground cleaning on the estimation of the HI abundance and bias parameters. The accompanying python toolkit is available at https://github.com/IntensityTools/MultipoleExpansion, and includes an interactive suite of examples to aid new users.
△ Less
Submitted 18 June, 2020; v1 submitted 13 February, 2020;
originally announced February 2020.
-
The cosmic abundance of cold gas in the local Universe
Authors:
Thomas J. Fletcher,
Amelie Saintonge,
Paula S. Soares,
Andrew Pontzen
Abstract:
We determine the cosmic abundance of molecular hydrogen (H2) in the local universe from the xCOLD GASS survey. To constrain the H2 mass function at low masses and correct for the effect of the lower stellar mass limit of 10^9 Msun in the xCOLD GASS survey, we use an empirical approach based on an observed scaling relation between star formation rate and gas mass. We also constrain the HI and HI+H2…
▽ More
We determine the cosmic abundance of molecular hydrogen (H2) in the local universe from the xCOLD GASS survey. To constrain the H2 mass function at low masses and correct for the effect of the lower stellar mass limit of 10^9 Msun in the xCOLD GASS survey, we use an empirical approach based on an observed scaling relation between star formation rate and gas mass. We also constrain the HI and HI+H2 mass functions using the xGASS survey, and compare it to the HI mass function from the ALFALFA survey. We find the cosmic abundance of molecular gas in the local Universe to be Omega_H2=(5.34+/-0.47)x10^-5 h^-1. Molecular gas accounts for 19.6 +/- 3.9% of the total abundance of cold gas, Omega_HI+H2=(4.66+/-0.70)x10^-4 h^-1. Galaxies with stellar masses in excess of 10^9 Msun account for 89% of the molecular gas in the local Universe, while in comparison such galaxies only contain 73% of the cold atomic gas as traced by the HI 21cm line. The xCOLD GASS CO, molecular gas and cold gas mass functions and Omega_H2 measurements provide constraints for models of galaxy evolution and help to anchor blind ALMA and NOEMA surveys attempting to determine the abundance of molecular gas at high redshifts.
△ Less
Submitted 12 February, 2020;
originally announced February 2020.
-
Full orbital solution for the binary system in the northern Galactic disc microlensing event Gaia16aye
Authors:
Łukasz Wyrzykowski,
P. Mróz,
K. A. Rybicki,
M. Gromadzki,
Z. Kołaczkowski,
M. Zieliński,
P. Zieliński,
N. Britavskiy,
A. Gomboc,
K. Sokolovsky,
S. T. Hodgkin,
L. Abe,
G. F. Aldi,
A. AlMannaei,
G. Altavilla,
A. Al Qasim,
G. C. Anupama,
S. Awiphan,
E. Bachelet,
V. Bakıs,
S. Baker,
S. Bartlett,
P. Bendjoya,
K. Benson,
I. F. Bikmaev
, et al. (160 additional authors not shown)
Abstract:
Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We presen…
▽ More
Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.
△ Less
Submitted 28 October, 2019; v1 submitted 22 January, 2019;
originally announced January 2019.