Showing 1–50 of 127 results for author: Santos, M G

Self-similarity of some soluble relatively free groups

Authors: Adilson A. Berlatto, Alex C. Dantas, Tulio M. G. Santos

Abstract: In this paper we prove that a free nilpotent group of finite rank is transitive self-similar. In contrast, we prove that a free metabelian group of rank $r \geq 2$ is not transitive self-similar. In this paper we prove that a free nilpotent group of finite rank is transitive self-similar. In contrast, we prove that a free metabelian group of rank $r \geq 2$ is not transitive self-similar. △ Less

Submitted 15 June, 2024; originally announced June 2024.

Investigating Mutual Coupling in the Hydrogen Epoch of Reionization Array and Mitigating its Effects on the 21-cm Power Spectrum

Authors: E. Rath, R. Pascua, A. T. Josaitis, A. Ewall-Wice, N. Fagnoni, E. de Lera Acedo, Z. E. Martinot, Z. Abdurashidova, T. Adams, J. E. Aguirre, R. Baartman, A. P. Beardsley, L. M. Berkhout, G. Bernardi, T. S. Billings, J. D. Bowman, P. Bull, J. Burba, R. Byrne, S. Carey, K. -F. Chen, S. Choudhuri, T. Cox, D. R. DeBoer, M. Dexter , et al. (56 additional authors not shown)

Abstract: Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategi… ▽ More Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 19 pages, 12 figures, submitted to MNRAS

On wreath product occurring as subgroup of automata group

Authors: Alex C. Dantas, Junio R. Oliveira, Tulio M. G. Santos

Abstract: A finitely generated group is said to be an automata group if it admits a faithful self-similar finite-state representation on some regular $m$-tree. We prove that if $G$ is a subgroup of an automata group, then for each finitely generated abelian group $A$, the wreath product $A \wr G$ is a subgroup of an automata group. We obtain, for example, that $C_2 \wr (C_{2} \wr \mathbb{Z})$,… ▽ More A finitely generated group is said to be an automata group if it admits a faithful self-similar finite-state representation on some regular $m$-tree. We prove that if $G$ is a subgroup of an automata group, then for each finitely generated abelian group $A$, the wreath product $A \wr G$ is a subgroup of an automata group. We obtain, for example, that $C_2 \wr (C_{2} \wr \mathbb{Z})$, $\mathbb{Z} \wr (C_2 \wr \mathbb{Z})$, $C_2 \wr (\mathbb{Z} \wr \mathbb{Z})$, and $\mathbb{Z} \wr (\mathbb{Z} \wr \mathbb{Z})$ are subgroups of automata groups. In the particular case $\mathbb{Z} \wr (\mathbb{Z} \wr \mathbb{Z})$, we prove that it is a subgroup of a two-letters automata group; this solves Problem 15.19 - (b) of the Kourovka Notebook proposed by A. M. Brunner and S. Sidki in 2000 [8, 17]. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Radio Frequency Interference from Radio Navigation Satellite Systems: simulations and comparison to MeerKAT single-dish data

Authors: Brandon Engelbrecht, Mario G. Santos, José Fonseca, Yichao Li, Jingying Wang, Melis O. Irfan, Stuart E. Harper, Keith Grainge, Philip Bull, Isabella P. Carucci, Steven Cunnington, Alkistis Pourtsidou, Marta Spinelli, Laura Wolz

Abstract: 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 source… ▽ More 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. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: 20 pages, 27 figures, 4 tables, Comments welcome

A demonstration of the effect of fringe-rate filtering in the Hydrogen Epoch of Reionization Array delay power spectrum pipeline

Authors: Hugh Garsden, Philip Bull, Mike Wilensky, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter , et al. (72 additional authors not shown)

Abstract: Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correl… ▽ More Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correlated with the rotating sky vs. those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes -- beam sidelobe/horizon suppression (the `mainlobe' filter), and ground-locked systematics removal (the `notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 21 pages, 18 figures, submitted to Monthly Notices of the Royal Astronomical Society

Hydrogen Epoch of Reionization Array (HERA) Phase II Deployment and Commissioning

Authors: Lindsay M. Berkhout, Daniel C. Jacobs, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (71 additional authors not shown)

Abstract: This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system an… ▽ More This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Journal ref: PASP 2024 136 045002

matvis: A matrix-based visibility simulator for fast forward modelling of many-element 21 cm arrays

Authors: Piyanat Kittiwisit, Steven G. Murray, Hugh Garsden, Philip Bull, Christopher Cain, Aaron R. Parsons, Jackson Sipple, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng , et al. (73 additional authors not shown)

Abstract: Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability… ▽ More Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 25 pages, 20 figures, submitted to RAS Techniques and Instruments, matvis is publicly available at https://github.com/HERA-Team/matvis

Bayesian estimation of cross-coupling and reflection systematics in 21cm array visibility data

Authors: Geoff G. Murphy, Philip Bull, Mario G. Santos, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Christopher Cain, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon, Nico Eksteen , et al. (54 additional authors not shown)

Abstract: Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method all… ▽ More Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages, 14 figures, submitted to MNRAS

MIGHTEE: multi-wavelength counterparts in the COSMOS field

Authors: I. H. Whittam, M. Prescott, C. L. Hale, M . J. Jarvis, I. Heywood, Fangxia An, M. Glowacki, N. Maddox, L. Marchetti, L. K. Morabito, N. J. Adams, R. A. A. Bowler, P. W. Hatfield, R. G. Varadaraj, J. Collier, B. Frank, A. R. Taylor, M. G. Santos, M. Vaccari, J. Afonso, Y. Ao, J. Delhaize, K. Knowles, S. Kolwa, S. M. Randriamampandry , et al. (4 additional authors not shown)

Abstract: In this paper we combine the Early Science radio continuum data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Survey, with optical and near-infrared data and release the cross-matched catalogues. The radio data used in this work covers $0.86$ deg$^2$ of the COSMOS field, reaches a thermal noise of $1.7$ $μ$Jy/beam and contains $6102$ radio components. We visually in… ▽ More In this paper we combine the Early Science radio continuum data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Survey, with optical and near-infrared data and release the cross-matched catalogues. The radio data used in this work covers $0.86$ deg$^2$ of the COSMOS field, reaches a thermal noise of $1.7$ $μ$Jy/beam and contains $6102$ radio components. We visually inspect and cross-match the radio sample with optical and near-infrared data from the Hyper Suprime-Cam (HSC) and UltraVISTA surveys. This allows the properties of active galactic nuclei and star-forming populations of galaxies to be probed out to $z \approx 5$. Additionally, we use the likelihood ratio method to automatically cross-match the radio and optical catalogues and compare this to the visually cross-matched catalogue. We find that 94 per cent of our radio source catalogue can be matched with this method, with a reliability of $95$ per cent. We proceed to show that visual classification will still remain an essential process for the cross-matching of complex and extended radio sources. In the near future, the MIGHTEE survey will be expanded in area to cover a total of $\sim$20~deg$^2$; thus the combination of automated and visual identification will be critical. We compare redshift distribution of SFG and AGN to the SKADS and T-RECS simulations and find more AGN than predicted at $z \sim 1$. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS

Emulating radiation transport on cosmological scale using a denoising Unet

Authors: Mosima P. Masipa, Sultan Hassan, Mario G. Santos, Gabriella Contardo, Kyunghyun Cho

Abstract: Semi-numerical simulations are the leading candidates for evolving reionization on cosmological scales. These semi-numerical models are efficient in generating large-scale maps of the 21cm signal, but they are too slow to enable inference at the field level. We present different strategies to train a U-Net to accelerate these simulations. We derive the ionization field directly from the initial de… ▽ More Semi-numerical simulations are the leading candidates for evolving reionization on cosmological scales. These semi-numerical models are efficient in generating large-scale maps of the 21cm signal, but they are too slow to enable inference at the field level. We present different strategies to train a U-Net to accelerate these simulations. We derive the ionization field directly from the initial density field without using the ionizing sources' location, and hence emulating the radiative transfer process. We find that the U-Net achieves higher accuracy in reconstructing the ionization field if the input includes either white noise or a noisy version of the ionization map beside the density field during training. Our model reconstructs the power spectrum over all scales perfectly well. This work represents a step towards generating large-scale ionization maps with a minimal cost and hence enabling rapid parameter inference at the field level. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 5 pages, 4 figures, Accepted at the ICLR 2023 Workshop on Physics for Machine Learning (Camera-ready version)

Search for the Epoch of Reionisation with HERA: Upper Limits on the Closure Phase Delay Power Spectrum

Authors: Pascal M. Keller, Bojan Nikolic, Nithyanandan Thyagarajan, Chris L. Carilli, Gianni Bernardi, Ntsikelelo Charles, Landman Bester, Oleg M. Smirnov, Nicholas S. Kern, Joshua S. Dillon, Bryna J. Hazelton, Miguel F. Morales, Daniel C. Jacobs, Aaron R. Parsons, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (58 additional authors not shown)

Abstract: Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standa… ▽ More Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 16 pages, 14 figures, accepted for publication by MNRAS

The foreground transfer function for HI intensity mapping signal reconstruction: MeerKLASS and precision cosmology applications

Authors: Steven Cunnington, Laura Wolz, Philip Bull, Isabella P. Carucci, Keith Grainge, Melis O. Irfan, Yichao Li, Alkistis Pourtsidou, Mario G. Santos, Marta Spinelli, Jingying Wang

Abstract: 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 intensi… ▽ More 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). △ Less

Submitted 23 May, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 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

Mitigating the effect of 1/f noise on the detection of the HI intensity mapping power spectrum from single-dish measurements

Authors: Melis O. Irfan, Yichao Li, Mario G. Santos, Philip Bull, Junhua Gu, Steven Cunnington, Keith Grainge, Jingying Wang

Abstract: We present and compare several methods to mitigate time-correlated (1/f) noise within the HI intensity mapping component of the MeerKAT Large Area Synoptic Survey (MeerKLASS). By simulating scan strategies, the HI signal, foreground emissions, white and correlated noise, we assess the ability of various data processing pipelines to recover the power spectrum of HI brightness temperature fluctuatio… ▽ More We present and compare several methods to mitigate time-correlated (1/f) noise within the HI intensity mapping component of the MeerKAT Large Area Synoptic Survey (MeerKLASS). By simulating scan strategies, the HI signal, foreground emissions, white and correlated noise, we assess the ability of various data processing pipelines to recover the power spectrum of HI brightness temperature fluctuations. We use MeerKAT pilot data to assess the level of 1/f noise expected for the MeerKLASS survey and use these measurements to create realistic levels of time-correlated noise for our simulations. We find the time-correlated noise component within the pilot data to be between 10 and 20 times higher than the white noise level at the scale of k = 0.04 Mpc^-1. Having determined that the MeerKAT 1/f noise is partially correlated across all the frequency channels, we employ Singular Value Decomposition (SVD) as a technique to remove both the 1/f noise and Galactic foregrounds but find that over-cleaning results in the removal of HI power at large (angular and radial) scales; a power loss of 40 per cent is seen for a 3-mode SVD clean at the scale of k = 0.04 Mpc^-1. We compare the impact of map-making using weighting by the full noise covariance (i.e. including a 1/f component), as opposed to just a simple unweighted binning, finding that including the time-correlated noise information reduces the excess power added by 1/f noise by up to 30 per cent. △ Less

Submitted 9 November, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 14 pages, 12 figures, accepted by MNRAS. This replacement features updated simulations which now include foreground emissions, telescope slew time and a receiver and elevation-dependant temperature contribution to each receiver which changes over frequency

A first detection of neutral hydrogen intensity mapping on Mpc scales at $z\approx 0.32$ and $z\approx 0.44$

Authors: Sourabh Paul, Mario G. Santos, Zhaoting Chen, Laura Wolz

Abstract: 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 underlyi… ▽ More 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σ$ and $11.5σ$ respectively). We constrain the rms of the fluctuations of the HI distribution to be $σ_{\rm HI} = (0.44\pm 0.04)\,{\rm mK}$ and $σ_{\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. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 19 pages, 10 figures, submitted to ApJL

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

Authors: The HERA Collaboration, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Rennan Barkana, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Daniela Breitman, Philip Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (70 additional authors not shown)

Abstract: We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that… ▽ More We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $Δ^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars. △ Less

Submitted 19 January, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Journal ref: 2023 ApJ 945 124

MIGHTEE-HI: The HI mass-stellar mass relation over the last billion years

Authors: Hengxing Pan, Matt J. Jarvis, Mario G. Santos, Natasha Maddox, Bradley S. Frank, Anastasia A. Ponomareva, Isabella Prandoni, Sushma Kurapati, Maarten Baes, Pavel E. Mancera Piña, Giulia Rodighiero, Martin J. Meyer, Romeel Davé, Gauri Sharma, Sambatriniaina H. A. Rajohnson, Nathan J. Adams, Rebecca A. A. Bowler, Francesco Sinigaglia, Thijs van der Hulst, Peter W. Hatfield, Srikrishna Sekhar, Jordan D. Collier

Abstract: We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and fin… ▽ More We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and find that the non-linear model is preferred over the linear one with a measured transition stellar mass of $\log_{10}(M_\star$/$M_{\odot})$ = $9.15\pm0.87$, beyond which the slope flattens. This finding supports the view that the lack of HI gas is ultimately responsible for the decreasing star formation rate observed in the massive main-sequence galaxies. For spirals alone, which are biased towards the massive galaxies in our sample, the slope beyond the transition mass is shallower than for the full sample, indicative of distinct gas processes ongoing for the spirals/high-mass galaxies from other types with lower stellar masses. We then create mock catalogues for the MIGHTEE-HI detections and non-detections with two main galaxy populations of late- and early-type galaxies to measure the underlying $M_{\rm HI}-M_{\star}$ relation. We find that the turnover in this relation persists whether considering the two galaxy populations as a whole or separately. We note that an underlying linear relation could mimic this turnover in the observed scaling relation, but a model with a turnover is strongly preferred. Measurements on the logarithmic average of HI masses against the stellar mass are provided as a benchmark for future studies. △ Less

Submitted 27 July, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: Revised manuscript accepted for publication in MNRAS

Connecting MeerKAT radio continuum properties to GAMA optical emission-line and WISE mid-infrared activity

Authors: H. F. M. Yao, M. E. Cluver, T. H. Jarrett, Gyula I. G. Jozsa, M. G. Santos, L. Marchetti, M. J. I. Brown, Y. A. Gordon, S. Brough, A. M. Hopkins, B. W. Holwerda, S. P. Driver, E. M. Sadler

Abstract: The identification of AGN in large surveys has been hampered by seemingly discordant classifications arising from differing diagnostic methods, usually tracing distinct processes specific to a particular wavelength regime. However, as shown in Yao et al. (2020), the combination of optical emission line measurements and mid-infrared photometry can be used to optimise the discrimination capability b… ▽ More The identification of AGN in large surveys has been hampered by seemingly discordant classifications arising from differing diagnostic methods, usually tracing distinct processes specific to a particular wavelength regime. However, as shown in Yao et al. (2020), the combination of optical emission line measurements and mid-infrared photometry can be used to optimise the discrimination capability between AGN and star formation activity. In this paper we test our new classification scheme by combining the existing GAMA-WISE data with high-quality MeerKAT radio continuum data covering 8 deg$^2$ of the GAMA G23 region. Using this sample of 1 841 galaxies (z < 0.25), we investigate the total infrared (derived from 12$μ$m) to radio luminosity ratio, q(TIR), and its relationship to optical-infrared AGN and star-forming (SF) classifications. We find that while q(TIR) is efficient at detecting AGN activity in massive galaxies generally appearing quiescent in the infrared, it becomes less reliable for cases where the emission from star formation in the host galaxy is dominant. However, we find that the q(TIR) can identify up to 70 % more AGNs not discernible at optical and/or infrared wavelengths. The median q(TIR) of our SF sample is 2.57 $\pm$ 0.23 consistent with previous local universe estimates. △ Less

Submitted 2 August, 2022; originally announced August 2022.

MIGHTEE-HI: Evolution of HI scaling relations of star-forming galaxies at $z<0.5$

Authors: Francesco Sinigaglia, Giulia Rodighiero, Ed Elson, Mattia Vaccari, Natasha Maddox, Bradley S. Frank, Matt J. Jarvis, Tom Oosterloo, Romeel Davé, Mara Salvato, Maarten Baes, Sabine Bellstedt, Laura Bisigello, Jordan D. Collier, Robin H. W. Cook, Luke J. M. Davies, Jacinta Delhaize, Simon P. Driver, Caroline Foster, Sushma Kurapati, Claudia del P. Lagos, Christopher Lidman, Pavel E. Mancera Piña, Martin J. Meyer, K. Moses Mogotsi , et al. (11 additional authors not shown)

Abstract: We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with… ▽ More We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with ${\rm sSFR}\equiv M_*/{\rm SFR}$), obtaining $\gtrsim 5σ$ detections in most cases, the strongest HI-stacking detections to date in this redshift range. With these detections, we are able to measure scaling relations in the probed redshift interval, finding evidence for a moderate evolution from the median redshift of our sample $z_{\rm med}\sim 0.37$ to $z\sim 0$. In particular, low-$M_*$ galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 9$) experience a strong HI depletion ($\sim 0.5$ dex in $\log_{10}(M_{\rm HI}/{\rm M}_\odot)$), while massive galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 11$) keep their HI mass nearly unchanged. When looking at the star formation activity, highly star-forming galaxies evolve significantly in $M_{\rm HI}$ ($f_{\rm HI}$, where $f_{\rm HI}\equiv M_{\rm}/M_*$) at fixed SFR (sSFR), while at the lowest probed SFR (sSFR) the scaling relations show no evolution. These findings suggest a scenario in which low-$M_*$ galaxies have experienced a strong HI depletion during the last $\sim4$ Gyr, while massive galaxies have undergone a significant HI replenishment through some accretion mechanism, possibly minor mergers. Interestingly, our results are in good agreement with the predictions of the SIMBA simulation. We conclude that this work sets novel important observational constraints on galaxy scaling relations. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures. Accepted for publication in ApJ Letters on 30 July 2022

Antenna characterization for the HIRAX experiment

Authors: Emily R. Kuhn, Benjamin R. B. Saliwanchik, Kevin Bandura, Michele Bianco, H. Cynthia Chiang, Devin Crichton, Meiling Deng, Sindhu Gaddam, Kit Gerodias, Austin Gumba, Maile Harris, Kavilan Moodley, V. Mugundhan, Laura Newburgh, Jeffrey Peterson, Elizabeth Pieters, Anna R. Polish, Alexandre Refregier, Ajith Sampath, Mario G. Santos, Onkabetse Sengate, Jonathan Sievers, Ema Smith, Will Tyndall, Anthony Walters , et al. (2 additional authors not shown)

Abstract: The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) aims to improve constraints on the dark energy equation of state through measurements of large-scale structure at high redshift ($0.8<z<2.5$), while serving as a state-of-the-art fast radio burst detector. Bright galactic foregrounds contaminate the 400--800~MHz HIRAX frequency band, so meeting the science goals will require precise… ▽ More The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) aims to improve constraints on the dark energy equation of state through measurements of large-scale structure at high redshift ($0.8<z<2.5$), while serving as a state-of-the-art fast radio burst detector. Bright galactic foregrounds contaminate the 400--800~MHz HIRAX frequency band, so meeting the science goals will require precise instrument characterization. In this paper we describe characterization of the HIRAX antenna, focusing on measurements of the antenna beam and antenna noise temperature. Beam measurements of the current HIRAX antenna design were performed in an anechoic chamber and compared to simulations. We report measurement techniques and results, which find a broad and symmetric antenna beam for $ν<$650MHz, and elevated cross-polarization levels and beam asymmetries for $ν>$700MHz. Noise temperature measurements of the HIRAX feeds were performed in a custom apparatus built at Yale. In this system, identical loads, one cryogenic and the other at room temperature, are used to take a differential (Y-factor) measurement from which the noise of the system is inferred. Several measurement sets have been conducted using the system, involving CHIME feeds as well as four of the HIRAX active feeds. These measurements give the first noise temperature measurements of the HIRAX feed, revealing a $\sim$60K noise temperature (relative to 30K target) with 40K peak- to-peak frequency-dependent features, and provide the first demonstration of feed repeatability. Both findings inform current and future feed designs. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 20 pages, 14 figures, SPIE proceedings

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.

Comments: 11 pages, 7 figures. See Fig.5 for the main result, showing the cross-power spectrum between MeerKAT and WiggleZ, revealing the 7.7σdetection. v2 is the version accepted for publication by MNRAS. The main changes from v1 were in introducing the radio intensity mapping data

Simulations of primary beam effects on the cosmic bispectrum phase observed with the Hydrogen Epoch of Reionization Array

Authors: N. Charles, G. Bernardi, H. L. Bester, O. M. Smirnov, C. Carilli, P. M. Keller, N. Kern, B. Nikolic, N. Thyagarajan, E. de Lera Acedo, N. Fagnoni, M. G. Santos

Abstract: The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly a… ▽ More The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly antenna primary beam responses. This work investigates the impact of antenna primary beams affected by mutual coupling on the closure phase and its power spectrum. Our simulations show that primary beams affected by mutual coupling lead to a leakage of foreground power into the EoR window, which can be up to $\sim4$ orders magnitude higher than the case where no mutual coupling is considered. This leakage is, however, essentially confined at $k < 0.3$~$h$~Mpc$^{-1}$ for triads that include 29~m baselines. The leakage magnitude is more pronounced when bright foregrounds appear in the antenna sidelobes, as expected. Finally, we find that triads that include mutual coupling beams different from each other have power spectra similar to triads that include the same type of mutual coupling beam, indicating that beam-to-beam variation within triads (or visibility pairs) is not the major source of foreground leakage in the EoR window. △ Less

Submitted 19 May, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Measurements of the diffuse Galactic synchrotron spectral index and curvature from MeerKLASS pilot data

Authors: Melis O. Irfan, Philip Bull, Mario G. Santos, Jingying Wang, Keith Grainge, Yichao Li, Isabella P. Carucci, Marta Spinelli, Steven Cunnington

Abstract: 21cm intensity mapping experiments are bringing an influx of high spectral resolution observational data in the $\sim100$ MHz $- 1$ GHz regime. We use pilot $971-1075$ MHz data from MeerKAT in single-dish mode, recently used to test the calibration and data reduction scheme of the upcoming MeerKLASS survey, to probe the spectral index of diffuse synchrotron emission below 1 GHz within… ▽ More 21cm intensity mapping experiments are bringing an influx of high spectral resolution observational data in the $\sim100$ MHz $- 1$ GHz regime. We use pilot $971-1075$ MHz data from MeerKAT in single-dish mode, recently used to test the calibration and data reduction scheme of the upcoming MeerKLASS survey, to probe the spectral index of diffuse synchrotron emission below 1 GHz within $145^{\circ} < α< 180^{\circ}$, $-1^{\circ} < δ< 8^{\circ}$. Through comparisons with data from the OVRO Long Wavelength Array and the Maipu and MU surveys, we find an average spectral index of $-2.75 < β< -2.71$ between 45 and 1055 MHz. By fitting for spectral curvature with a spectral index of the form $β+ c \, {\rm{ln}}(ν/ 73~{\rm MHz})$, we measure $β= -2.55 \pm 0.13$ and $c = -0.12 \pm 0.05$ within our target field. Our results are in good agreement (within $1σ$) with existing measurements from experiments such as ARCADE2 and EDGES. These results show the calibration accuracy of current data and demonstrate that MeerKLASS will also be capable of achieving a secondary science goal of probing the interstellar medium. △ Less

Submitted 11 January, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS. Updated to match published paper (additional references and acknowledgements)

Self-similar abelian groups and their centralizers

Authors: Alex C. Dantas, Tulio M. G. Santos, Said N. Sidki

Abstract: We extend results on transitive self-similar abelian subgroups of the group of automorphisms $\mathcal{A}_m$ of an $m$-ary tree $\mathcal{T}_m$ in \cite{BS}, to the general case where the permutation group induced on the first level of the tree has $s\geq 1$ orbits. We prove that such a group $A$ embeds in a self-similar abelian group $A^*$ which is also a maximal abelian subgroup of… ▽ More We extend results on transitive self-similar abelian subgroups of the group of automorphisms $\mathcal{A}_m$ of an $m$-ary tree $\mathcal{T}_m$ in \cite{BS}, to the general case where the permutation group induced on the first level of the tree has $s\geq 1$ orbits. We prove that such a group $A$ embeds in a self-similar abelian group $A^*$ which is also a maximal abelian subgroup of $\mathcal{A}_m$. The construction of $A^*$ is based on the definition of a free monoid $Δ$ of rank $s$ of partial diagonal monomorphisms of $\mathcal{A}_m$, which is used to determine the structure of $C_{\mathcal{A}_m}(A)$, the centralizer of $A$ in $\mathcal{A}_m$. Indeed, we prove $A^*=C_{\mathcal{A}_m} (Δ(A))=\overline{ Δ({B(A)})}$, where $B(A)$ denotes the product of the projections of $A$ in its action on the different $s$ orbits of maximal subtrees of $\mathcal{T}_m$ and bar denotes the topological closure. When $A$ is a torsion self-similar abelian group, it is shown that it is necessarily of finite exponent. Moreover, we extend recent constructions of self-similar free abelian groups of infinite enumerable rank to examples of such groups which are also $Δ$-invariant for $s=2$. Finally, we focus on self-similar cyclic groups of automorphisms of $\mathcal{T}_m$ and compute their centralizers when $m=4.$ △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 25 pages

MIGHTEE: total intensity radio continuum imaging and the COSMOS / XMM-LSS Early Science fields

Authors: I. Heywood, M. J. Jarvis, C. L. Hale, I. H. Whittam, H. L. Bester, B. Hugo, J. S. Kenyon, M. Prescott, O. M. Smirnov, C. Tasse, J. M. Afonso, P. N. Best, J. D. Collier, R. P. Deane, B. S. Frank, M. J. Hardcastle, K. Knowles, N. Maddox, E. J. Murphy, I. Prandoni, S. M. Randriamampandry, M. G. Santos, S. Sekhar, F. Tabatabaei, A. R. Taylor , et al. (1 additional authors not shown)

Abstract: MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectro-polarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image $\sim$20 deg$^{2}$ over the COSMOS, E-CDFS, ELAIS-S1, and XMM-LSS extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength dataset… ▽ More MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectro-polarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image $\sim$20 deg$^{2}$ over the COSMOS, E-CDFS, ELAIS-S1, and XMM-LSS extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength datasets from numerous existing and forthcoming observational campaigns. Here we describe and validate the data processing strategy for the total intensity continuum aspect of MIGHTEE, using a single deep pointing in COSMOS (1.6 deg$^{2}$) and a three-pointing mosaic in XMM-LSS (3.5 deg$^{2}$). The processing includes the correction of direction-dependent effects, and results in thermal noise levels below 2~$\mathrmμ$Jy beam$^{-1}$ in both fields, limited in the central regions by classical confusion at $\sim$8$''$ angular resolution, and meeting the survey specifications. We also produce images at $\sim$5$''$ resolution that are $\sim$3 times shallower. The resulting image products form the basis of the Early Science continuum data release for MIGHTEE. From these images we extract catalogues containing 9,896 and 20,274 radio components in COSMOS and XMM-LSS respectively. We also process a close-packed mosaic of 14 additional pointings in COSMOS and use these in conjunction with the Early Science pointing to investigate methods for primary beam correction of broadband radio images, an analysis that is of relevance to all full-band MeerKAT continuum observations, and wide field interferometric imaging in general. A public release of the MIGHTEE Early Science continuum data products accompanies this article. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: 21 pages, 13 figures, 3 tables, to be published in MNRAS

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

Authors: Devin Crichton, Moumita Aich, Adam Amara, Kevin Bandura, Bruce A. Bassett, Carlos Bengaly, Pascale Berner, Shruti Bhatporia, Martin Bucher, Tzu-Ching Chang, H. Cynthia Chiang, Jean-Francois Cliche, Carolyn Crichton, Romeel Dave, Dirk I. L. de Villiers, Matt A. Dobbs, Aaron M. Ewall-Wice, Scott Eyono, Christopher Finlay, Sindhu Gaddam, Ken Ganga, Kevin G. Gayley, Kit Gerodias, Tim Gibbon, Austin Gumba , et al. (75 additional authors not shown)

Abstract: The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 40… ▽ More The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array. △ Less

Submitted 17 January, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 24 pages, 6 figures. Updated to match published version

Journal ref: J. of Astronomical Telescopes, Instruments, and Systems, 8(1), 011019 (2022)

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

Authors: Dara Storer, Joshua S. Dillon, Daniel C. Jacobs, Miguel F. Morales, Bryna J. Hazelton, Aaron Ewall-Wice, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Scott Dynes , et al. (53 additional authors not shown)

Abstract: We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for… ▽ More We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets. △ Less

Submitted 4 May, 2022; v1 submitted 26 September, 2021; originally announced September 2021.

Comments: 31 pages, 17 figures

Journal ref: Radio Science, vol. 57, no. 1, 2022

1/f Noise Analysis for FAST HI Intensity Mapping Drift-Scan Experiment

Authors: Wenkai Hu, Yichao Li, Yougang Wang, Fengquan Wu, Bo Zhang, Ming Zhu, Shifan Zuo, Guilaine Lagache, Yinzhe Ma, Mario G. Santos, Xuelei Chen

Abstract: We investigate the 1/f noise of the Five-hundred-meter Aperture Spherical Telescope (FAST) receiver system using drift-scan data from an intensity mapping pilot survey. All the 19 beams have 1/f fluctuations with similar structures. Both the temporal and the 2D power spectrum densities are estimated. The correlations directly seen in the time series data at low frequency $f$ are associated with th… ▽ More We investigate the 1/f noise of the Five-hundred-meter Aperture Spherical Telescope (FAST) receiver system using drift-scan data from an intensity mapping pilot survey. All the 19 beams have 1/f fluctuations with similar structures. Both the temporal and the 2D power spectrum densities are estimated. The correlations directly seen in the time series data at low frequency $f$ are associated with the sky signal, perhaps due to a coupling between the foreground and the system response. We use Singular Value Decomposition (SVD) to subtract the foreground. By removing the strongest components, the measured 1/f noise power can be reduced significantly. With 20 modes subtraction, the knee frequency of the 1/f noise in a 10 MHz band is reduced to $1.8 \times 10^{-3}\Hz$, well below the thermal noise over 500-seconds time scale. The 2D power spectra show that the 1/f-type variations are restricted to a small region in the time-frequency space and the correlations in frequency can be suppressed with SVD modes subtraction. The residual 1/f noise after the SVD mode subtraction is uncorrelated in frequency, and a simple noise diode frequency-independent calibration of the receiver gain at 8s interval does not affect the results. The 1/f noise can be important for HI intensity mapping, we estimate that the 1/f noise has a knee frequency $(f_{k}) \sim$ 6 $\times$ 10$^{-4}$Hz, and time and frequency correlation spectral indices $(α) \sim 0.65$, $(β) \sim 0.8$ after the SVD subtraction of 30 modes. This can bias the HI power spectrum measurement by 10 percent. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 13 pages,17 figures, MNRAS accepted

Measuring the baryonic Tully-Fisher relation below the detection threshold

Authors: Hengxing Pan, Matt J. Jarvis, Anastasia A. Ponomareva, Mario G. Santos, James R. Allison, Natasha Maddox, Bradley S. Frank

Abstract: We present a novel 2D flux density model for observed HI emission lines combined with a Bayesian stacking technique to measure the baryonic Tully-Fisher relation below the nominal detection threshold. We simulate a galaxy catalogue, which includes HI lines described either with Gaussian or busy function profiles, and HI data cubes with a range of noise and survey areas similar to the MeerKAT Inter… ▽ More We present a novel 2D flux density model for observed HI emission lines combined with a Bayesian stacking technique to measure the baryonic Tully-Fisher relation below the nominal detection threshold. We simulate a galaxy catalogue, which includes HI lines described either with Gaussian or busy function profiles, and HI data cubes with a range of noise and survey areas similar to the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. With prior knowledge of redshifts, stellar masses and inclinations of spiral galaxies, we find that our model can reconstruct the input baryonic Tully-Fisher parameters (slope and zero point) most accurately in a relatively broad redshift range from the local Universe to $z = 0.3$ for all the considered levels of noise and survey areas, and up to $z = 0.55$ for a nominal noise of $90\,μ$Jy/channel over 5 deg$^{2}$. Our model can also determine the $M_{\rm HI} - M_{\star}$ relation for spiral galaxies beyond the local Universe, and account for the detailed shape of the HI emission line, which is crucial for understanding the dynamics of spiral galaxies. Thus, we have developed a Bayesian stacking technique for measuring the baryonic Tully-Fisher relation for galaxies at low stellar and/or HI masses and/or those at high redshift, where the direct detection of HI requires prohibitive exposure times. △ Less

Submitted 9 September, 2021; originally announced September 2021.

Comments: 13 pages, 9 figures. Accepted for publication in MNRAS

Effects of model incompleteness on the drift-scan calibration of radio telescopes

Authors: Bharat K. Gehlot, Daniel C. Jacobs, Judd D. Bowman, Nivedita Mahesh, Steven G. Murray, Matthew Kolopanis, Adam P. Beardsley, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely , et al. (54 additional authors not shown)

Abstract: Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less at… ▽ More Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less attention. We investigate the sky based calibration of total power measurements with a HERA dish and an EDGES style antenna to understand the role of auto-correlations in the calibration of an interferometer and the role of sky in calibrating a total power instrument. Using simulations we study various scenarios such as time variable gain, incomplete sky calibration model, and primary beam model. We find that temporal gain drifts, sky model incompleteness, and beam inaccuracies cause biases in the receiver gain amplitude and the receiver temperature estimates. In some cases, these biases mix spectral structure between beam and sky resulting in spectrally variable gain errors. Applying the calibration method to the HERA and EDGES data, we find good agreement with calibration via the more standard methods. Although instrumental gains are consistent with beam and sky errors similar in scale to those simulated, the receiver temperatures show significant deviations from expected values. While we show that it is possible to partially mitigate biases due to model inaccuracies by incorporating a time-dependent gain model in calibration, the resulting errors on calibration products are larger and more correlated. Completely addressing these biases will require more accurate sky and primary beam models. △ Less

Submitted 15 July, 2021; v1 submitted 25 April, 2021; originally announced April 2021.

Comments: 16 pages, 13 figures, 1 table; accepted for publication in MNRAS main journal

Towards Automated Acceptance testing for industrial robots

Authors: Marcela G. dos Santos, Fabio Petrillo

Abstract: Industrial robots are important machines applied in numerous modern industries that execute repetitive tasks with high accuracy, replacing or supporting dangerous jobs. In this kind of system, with increased complexity in which cost is related to the time the system keeps working, the system must operate with a minimum number of failures. In other words, a quality aspect important in industry is r… ▽ More Industrial robots are important machines applied in numerous modern industries that execute repetitive tasks with high accuracy, replacing or supporting dangerous jobs. In this kind of system, with increased complexity in which cost is related to the time the system keeps working, the system must operate with a minimum number of failures. In other words, a quality aspect important in industry is reliability. We hypothesize that Automated Acceptance Testing improves reliability for industrial robot program. We present the research question, the motivation for this study, our hypothesis and future research efforts. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

Authors: James E. Aguirre, Steven G. Murray, Robert Pascua, Zachary E. Martinot, Jacob Burba, Joshua S. Dillon, Daniel C. Jacobs, Nicholas S. Kern, Piyanat Kittiwisit, Matthew Kolopanis, Adam Lanman, Adrian Liu, Lily Whitler, Zara Abdurashidova, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Steve Carey, Chris L. Carilli , et al. (51 additional authors not shown)

Abstract: We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the… ▽ More We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 32 pages, 20 figures. Submitted to the Astrophysical Journal

A Real Time Processing System for Big Data in Astronomy: Applications to HERA

Authors: Paul La Plante, Peter K. G. Williams, Matthew Kolopanis, Joshua S. Dillon, Adam P. Beardsley, Nicholas S. Kern, Michael Wilensky, Zaki S. Ali, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, John Ely , et al. (50 additional authors not shown)

Abstract: As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen… ▽ More As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen Epoch of Reionization Array (HERA), a radio interferometer endeavoring to provide the first detection of the highly redshifted 21 cm signal from Cosmic Dawn and the Epoch of Reionization by an interferometer. The RTP system consists of analysis routines run on raw data shortly after they are acquired, such as calibration and detection of radio-frequency interference (RFI) events. RTP works closely with the Librarian, the HERA data storage and transfer manager which automatically ingests data and transfers copies to other clusters for post-processing analysis. Both the RTP system and the Librarian are public and open source software, which allows for them to be modified for use in other scientific collaborations. When fully constructed, HERA is projected to generate over 50 terabytes (TB) of data each night, and the RTP system enables the successful scientific analysis of these data. △ Less

Submitted 30 September, 2021; v1 submitted 8 April, 2021; originally announced April 2021.

Comments: 15 pages, 3 figures, published in Astronomy and Computing

Journal ref: Astronomy and Computing 36 (2021) 100489

Software Engineering for Robotic Systems:a systematic mapping study

Authors: Marcela G. dos Santos, Fabio Petrillo

Abstract: Robots are being applied in a vast range of fields, leading researchers and practitioners to write tasks more complex than in the past. The robot software complexity increases the difficulty of engineering the robot's software components with quality requirements. Researchers and practitioners have applied software engineering (SE) approaches and robotic domains to address this issue in the last t… ▽ More Robots are being applied in a vast range of fields, leading researchers and practitioners to write tasks more complex than in the past. The robot software complexity increases the difficulty of engineering the robot's software components with quality requirements. Researchers and practitioners have applied software engineering (SE) approaches and robotic domains to address this issue in the last two decades. This study aims to identify, classify and evaluate the current state-of-the-art Software Engineering for Robotic Systems (SERS). We systematically selected and analyzed 50 primary studies extracted from an automated search on Scopus digital library and manual search on the two editions of the RoSE workshop. We present three main contributions. Firstly, we provide an analysis from three following perspectives: demographics of publication, SE areas applied in robotics domains, and RSE findings. Secondly, we show a catalogue of research studies that apply software engineering techniques in the robotic domain, classified with the SWEBOK guide. We have identified 5 of 15 software engineering areas from the SWEBOK guide applied explicitly in robotic domains. The majority of the studies focused on the development phase (design, models and methods and construction). Testing and quality software areas have little coverage in SERS. Finally, we identify research opportunities and gaps in software engineering for robotic systems for future studies. △ Less

Submitted 24 February, 2021; originally announced February 2021.

A deep radio view of the evolution of the cosmic star-formation rate density from a stellar-mass selected sample in VLA-COSMOS

Authors: Eliab D. Malefahlo, Matt J. Jarvis, Mario G. Santos, Sarah V. White, Nathan J. Adams, Rebecca A. A. Bowler

Abstract: We present the 1.4GHz radio luminosity functions (RLFs) of galaxies in the COSMOS field, measured above and below the $5σ$ detection threshold, using a Bayesian model-fitting technique. The radio flux-densities from VLA-COSMOS 3-GHz data, are extracted at the position of stellar mass-limited near-infrared (NIR) galaxies. We fit a local RLF model, which is a combination of active galactic nuclei (A… ▽ More We present the 1.4GHz radio luminosity functions (RLFs) of galaxies in the COSMOS field, measured above and below the $5σ$ detection threshold, using a Bayesian model-fitting technique. The radio flux-densities from VLA-COSMOS 3-GHz data, are extracted at the position of stellar mass-limited near-infrared (NIR) galaxies. We fit a local RLF model, which is a combination of active galactic nuclei (AGN) and star-forming galaxy (SFG), in 10 redshift bins with a pure luminosity evolution (PLE) model. We show that the evolution strength is similar to literature values up to $z\sim 1.6$. Beyond $z\sim 2$, we find that the SFG RLF exhibits a negative evolution ($L^*$ moves to lower luminosities) due to the decrease in low stellar-mass sources in our stellar mass-limited sample at high redshifts. From the RLF for SFGs, we determine the evolution in the cosmic star-formation-rate density (SFRD), which we find to be consistent with the established behaviour up to $z\sim 1$. Beyond $z\sim 1$ cosmic SFRD declines if one assumes an evolving infrared--radio correlation (IRRC), whereas it stays relatively higher if one adopts a constant IRRC. We find that the form of the relation between radio luminosity and SFR is therefore crucial in measuring the cosmic SFRD from radio data. We investigate the effects of stellar mass on the total RLF by splitting our sample into low ($10^{8.5} \leq M/\mathrm{M}_{\odot} \leq 10^{10}$) and high ($M>10^{10}\,\mathrm{M}_{\odot}$) stellar-mass subsets. We find that the SFRD is dominated by sources in the high stellar masses bin, at all redshifts. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: 28 pages, 16 figures, 6 tables

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.

Comments: 24 pages, 34 figures, accepted for publication in MNRAS

Hi intensity mapping with MeerKAT: Primary beam effects on foreground cleaning

Authors: Siyambonga D. Matshawule, Marta Spinelli, Mario G. Santos, Sibonelo Ngobese

Abstract: Upcoming and future neutral hydrogen Intensity Mapping surveys offer a great opportunity to constrain cosmology in the post-reionization Universe, provided a good accuracy is achieved in the separation between the strong foregrounds and the cosmological signal. Cleaning methods rely on the frequency smoothness of the foregrounds and are often applied under the assumption of a simplistic Gaussian p… ▽ More Upcoming and future neutral hydrogen Intensity Mapping surveys offer a great opportunity to constrain cosmology in the post-reionization Universe, provided a good accuracy is achieved in the separation between the strong foregrounds and the cosmological signal. Cleaning methods rely on the frequency smoothness of the foregrounds and are often applied under the assumption of a simplistic Gaussian primary beam. In this work, we test the cleaning in the presence of a realistic primary beam model with a non trivial frequency dependence. We focus on the Square Kilometre Array precursor MeerKAT telescope and simulate a single-dish wide area survey. We consider the main foreground components, including an accurate full sky point source catalogue. We find that the coupling between beam sidelobes and the foreground structure can complicate the cleaning. However, when the beam frequency dependence is smooth, we show that the cleaning is only problematic if the far sidelobes are unexpectedly large. Even in that case, a proper reconstruction is possible if the strongest point sources are removed and the cleaning is more aggressive. We then consider a non-trivial frequency dependence: a sinusoidal type feature in the beam width that is present in the MeerKAT beam and is expected in most dishes, including SKA1-MID. Such a feature, coupling with the foreground emission, biases the reconstruction of the signal across frequency, potentially impacting the cosmological analysis. Interestingly, such contamination is present at a lower level even when no point sources are included and the beam is Gaussian, showing that this frequency ripple can be problematic even within the main lobe. We show that this effect is constrained to a narrow region in $k_\parallel$ space and can be reduced if the maps are carefully re-smoothed to a common lower resolution. △ Less

Submitted 1 September, 2021; v1 submitted 21 November, 2020; originally announced November 2020.

Comments: Revised version accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 506, Issue 4, pp.5075-5092, 2021

MIGHTEE-HI: The HI emission project of the MeerKAT MIGHTEE survey

Authors: Natasha Maddox, Bradley S. Frank, A. A. Ponomareva, M. J. Jarvis, E. A. K. Adams, R. Davé, T. A. Oosterloo, M. G. Santos, S. L. Blyth, M. Glowacki, R. C. Kraan-Korteweg, W. Mulaudzi, B. Namumba, I. Prandoni, S. H. A. Rajohnson, K. Spekkens, N. J. Adams, R. A. A. Bowler, J. D. Collier, I. Heywood, S. Sekhar, A. R. Taylor

Abstract: We present the HI emission project within the MIGHTEE survey, currently being carried out with the newly commissioned MeerKAT radio telescope. This is one of the first deep, blind, medium-wide interferometric surveys for neutral hydrogen (HI) ever undertaken, extending our knowledge of HI emission to z=0.6. The science goals of this medium-deep, medium-wide survey are extensive, including the evol… ▽ More We present the HI emission project within the MIGHTEE survey, currently being carried out with the newly commissioned MeerKAT radio telescope. This is one of the first deep, blind, medium-wide interferometric surveys for neutral hydrogen (HI) ever undertaken, extending our knowledge of HI emission to z=0.6. The science goals of this medium-deep, medium-wide survey are extensive, including the evolution of the neutral gas content of galaxies over the past 5 billion years. Simulations predict nearly 3000 galaxies over 0<z<0.4 will be detected directly in HI, with statistical detections extending to z=0.6. The survey allows us to explore HI as a function of galaxy environment, with massive groups and galaxy clusters within the survey volume. Additionally, the area is large enough to contain as many as 50 local galaxies with HI mass $<10^8$ Msun, which allows us to study the low-mass galaxy population. The 20 deg$^2$ main survey area is centred on fields with exceptional multi-wavelength ancillary data, with photometry ranging from optical through far-infrared wavelengths, supplemented with multiple spectroscopic campaigns. We describe here the survey design and the key science goals. We also show first results from the Early Science observations, including kinematic modelling of individual sources, along with the redshift, HI, and stellar mass ranges of the sample to date. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: 14 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 646, A35 (2021)

HI intensity mapping with the MIGHTEE survey: power spectrum estimates

Authors: Sourabh Paul, Mario G. Santos, Junaid Townsend, Matt J. Jarvis, Natasha Maddox, Jordan D. Collier, Bradley S. Frank, Russ Taylor

Abstract: Intensity mapping (IM) with neutral hydrogen is a promising avenue to probe the large scale structure of the Universe. In this paper, we demonstrate that using the 64-dish MeerKAT radio telescope as a connected interferometer, it is possible to make a statistical detection of HI in the post-reionization Universe. With the MIGHTEE (MeerKAT International GHz Tiered Extragalactic Exploration) survey… ▽ More Intensity mapping (IM) with neutral hydrogen is a promising avenue to probe the large scale structure of the Universe. In this paper, we demonstrate that using the 64-dish MeerKAT radio telescope as a connected interferometer, it is possible to make a statistical detection of HI in the post-reionization Universe. With the MIGHTEE (MeerKAT International GHz Tiered Extragalactic Exploration) survey project observing in the L-band ($856 < ν< 1712$ MHz, $z < 0.66$), we can achieve the required sensitivity to measure the HI IM power spectrum on quasi-linear scales, which will provide an important complementarity to the single-dish IM MeerKAT observations. We present a purpose-built simulation pipeline that emulates the MIGHTEE observations and forecast the constraints that can be achieved on the HI power spectrum at $z = 0.27$ for $k > 0.3$ $\rm{Mpc}^{-1}$ using the foreground avoidance method. We present the power spectrum estimates with the current simulation on the COSMOS field that includes contributions from HI, noise and point source models constructed from the observed MIGHTEE data. The results from our \textit{visibility} based pipeline are in qualitative agreement to the already available MIGHTEE data. This paper demonstrates that MeerKAT can achieve very high sensitivity to detect HI with the full MIGHTEE survey on quasi-linear scales (signal-to-noise ratio $> 7$ at $k=0.49$ $\rm{Mpc}^{-1}$) which are instrumental in probing cosmological quantities such as the spectral index of fluctuation, constraints on warm dark matter, the quasi-linear redshift space distortions and the measurement of the HI content of the Universe up to $z\sim 0.5$. △ Less

Submitted 15 April, 2021; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: 13 pages, 8 figures, accepted for publication in MNRAS

Galaxy and Mass Assembly (GAMA): A $\textit{WISE}$ study of the activity of emission-line systems in G23

Authors: H. F. M. Yao, T. H. Jarrett, M. E. Cluver, L. Marchetti, Edward N. Taylor, M. G. Santos, Matt S. Owers, Angel R. Lopez-Sanchez, Y. A. Gordon, M. J. I. Brown, S. Brough, S. Phillipps, B. W. Holwerda, A. M. Hopkins, L. Wang

Abstract: We present a detailed study of emission-line systems in the GAMA G23 region, making use of $\textit{WISE}$ photometry that includes carefully measured resolved sources. After applying several cuts to the initial catalogue of $\sim$41,000 galaxies, we extract a sample of 9,809 galaxies. We then compare the spectral diagnostic (BPT) classification of 1154 emission-line galaxies (38$\%$ resolved in W… ▽ More We present a detailed study of emission-line systems in the GAMA G23 region, making use of $\textit{WISE}$ photometry that includes carefully measured resolved sources. After applying several cuts to the initial catalogue of $\sim$41,000 galaxies, we extract a sample of 9,809 galaxies. We then compare the spectral diagnostic (BPT) classification of 1154 emission-line galaxies (38$\%$ resolved in W1) to their location in the $\textit{WISE}$ colour-colour diagram, leading to the creation of a new zone for mid-infrared "warm" galaxies located 2$σ$ above the star-forming sequence, below the standard $\textit{WISE}$ AGN region. We find that the BPT and $\textit{WISE}$ diagrams agree on the classification for 85$\%$ and 8$\%$ of the galaxies as non-AGN (star forming = SF) and AGN, respectively, and disagree on $\sim$7$\%$ of the entire classified sample. 39$\%$ of the AGN (all types) are broad-line systems for which the [\ion{N}{ii}] and [H$α$] fluxes can barely be disentangled, giving in most cases spurious [\ion{N}{ii}]/[H$α$] flux ratios. However, several optical AGN appear to be completely consistent with SF in $\textit{WISE}$. We argue that these could be low power AGN, or systems whose hosts dominate the IR emission. Alternatively, given the sometimes high [\ion{O}{iii}] luminosity in these galaxies, the emission lines may be generated by shocks coming from super-winds associated with SF rather than the AGN activity. Based on our findings, we have created a new diagnostic: [W1-W2] vs [\ion{N}{ii}]/[H$α$], which has the virtue of separating SF from AGN and high-excitation sources. It classifies 3$\sim$5 times more galaxies than the classic BPT △ Less

Submitted 13 September, 2020; originally announced September 2020.

Comments: 43 pages, 32 figures, 4 tables

HI intensity mapping with MeerKAT: 1/f noise analysis

Authors: Yichao Li, Mario G. Santos, Keith Grainge, Stuart Harper, Jingying Wang

Abstract: The nature of the time correlated noise component (the 1/f noise) of single dish radio telescopes is critical to the detectability of the HI signal in intensity mapping experiments. In this paper, we present the 1/f noise properties of the MeerKAT receiver system using South Celestial Pole (SCP) tracking data. We estimate both the temporal power spectrum density and the 2D power spectrum density f… ▽ More The nature of the time correlated noise component (the 1/f noise) of single dish radio telescopes is critical to the detectability of the HI signal in intensity mapping experiments. In this paper, we present the 1/f noise properties of the MeerKAT receiver system using South Celestial Pole (SCP) tracking data. We estimate both the temporal power spectrum density and the 2D power spectrum density for each of the antennas and polarizations. We apply Singular Value Decomposition (SVD) to the dataset and show that, by removing the strongest components, the 1/f noise can be drastically reduced, indicating that it is highly correlated in frequency. Without SVD mode subtraction, the knee frequency over a $20\,$MHz integration is higher than $0.1\,\rm Hz$; with just $2$~mode subtraction, the knee frequency is reduced to $\sim 3\times 10^{-3}\,{\rm Hz}$, indicating that the system induced 1/f-type variations are well under the thermal noise fluctuations over a few hundred seconds time scales. The 2D power spectrum shows that the 1/f-type variations are restricted to a small region in the time-frequency space, either with long wavelength correlations in frequency or in time. This gives a wide range of cosmological scales where the 21cm signal can be measured without further need to calibrate the gain time fluctuations. Finally, we demonstrate that a simple power spectrum parameterization is sufficient to describe the data and provide fitting parameters for both the 1D and 2D power spectrum. △ Less

Submitted 3 July, 2020; originally announced July 2020.

Comments: 14 pages, 12 figures

Internet of Things Architectures: A Comparative Study

Authors: Marcela G. dos Santos, Darine Ameyed, Fabio Petrillo, Fehmi Jaafar, Mohamed Cheriet

Abstract: Over the past two decades, the Internet of Things (IoT) has become an underlying concept to a variety of solutions and technologies that it is now hardly possible to enumerate and describe all of them. The concept behind the Internet of Things is as powerful as it is complex, and for the components in the IoT solution tomesh together perfectly, they all have to be part of a well-thought-out struct… ▽ More Over the past two decades, the Internet of Things (IoT) has become an underlying concept to a variety of solutions and technologies that it is now hardly possible to enumerate and describe all of them. The concept behind the Internet of Things is as powerful as it is complex, and for the components in the IoT solution tomesh together perfectly, they all have to be part of a well-thought-out structure. That is where understanding the IoT architecture becomes paramount. Because of the vast domain of IoT, there is no single consensus on IoT architecture. Different researchers and organizations proposed different architectures under a variety of classifications, mainly: conceptual, standard and, industrial or commercial adoption. It is indispensable to make a systematic analysis of IoT architecture to be able to compare the industrial proposals and identify their similarities and their differences. In this work, we summarize information about seven IoT industrial architectures in order to propose an approach that makes possible a comparative analysis between different IoT architectures. This work presents two main contributions: (i) an approach for analyzing and comparing IoTarchitectures using Layer-Model; (ii) a comparative study of seven industrial IoT architectures. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Intransitive Self-similar Groups

Authors: Alex C. Dantas, Tulio M. G. Santos, Said N. Sidki

Abstract: A group is said to be self-similar provided it admits a faithful state-closed representation on some regular $m$-tree and the group is said to be transitive self-similar provided additionally it induces transitive action on the first level of the tree. A standard approach for constructing a transitive self-similar representation of a group has been by way of a single virtual endomorphism of \ the… ▽ More A group is said to be self-similar provided it admits a faithful state-closed representation on some regular $m$-tree and the group is said to be transitive self-similar provided additionally it induces transitive action on the first level of the tree. A standard approach for constructing a transitive self-similar representation of a group has been by way of a single virtual endomorphism of \ the group in question. Recently, it was shown that this approach when applied to the restricted wreath product $% \mathbb{Z}\wr \mathbb{Z}$ could not produce a faithful transitive self-similar representations for any $m\geq 2$ (see, \cite{DS}). In this work we study state-closed representations without assuming the transitivity condition. This general action is translated into a set of virtual endomorphisms corresponding to the different orbits of the action on the first level of the tree. In this manner, we produce faithful self-similar representations, some of which are also finite-state, for a number of groups such as $\mathbb{Z}^ω$, $\mathbb{Z}\wr \mathbb{Z}$ and $(\mathbb{Z} \wr \mathbb{Z}) \wr C_{2}$. △ Less

Submitted 21 April, 2020; v1 submitted 19 April, 2020; originally announced April 2020.

Foreground modelling via Gaussian process regression: an application to HERA data

Authors: Abhik Ghosh, Florent Mertens, Gianni Bernardi, Mário G. Santos, Nicholas S. Kern, Christopher L. Carilli, Trienko L. Grobler, Léon V. E. Koopmans, Daniel C. Jacobs, Adrian Liu, Aaron R. Parsons, Miguel F. Morales, James E. Aguirre, Joshua S. Dillon, Bryna J. Hazelton, Oleg M. Smirnov, Bharat K. Gehlot, Siyanda Matika, Paul Alexander, Zaki S. Ali, Adam P. Beardsley, Roshan K. Benefo, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (48 additional authors not shown)

Abstract: The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize… ▽ More The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz. △ Less

Submitted 12 May, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 15 pages, 15 figures, 1 table, Accepted to MNRAS

Redundant-Baseline Calibration of the Hydrogen Epoch of Reionization Array

Authors: Joshua S. Dillon, Max Lee, Zaki S. Ali, Aaron R. Parsons, Naomi Orosz, Chuneeta Devi Nunhokee, Paul La Plante, Adam P. Beardsley, Nicholas S. Kern, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo , et al. (54 additional authors not shown)

Abstract: In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of… ▽ More In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as "redundant-baseline calibration" resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions--both in data and in simulations--and present strategies for mitigating that structure. △ Less

Submitted 3 November, 2020; v1 submitted 18 March, 2020; originally announced March 2020.

Comments: 24 Pages, 19 Figures. Updated to match the accepted MNRAS version

Constraining the reionization history using deep learning from 21cm tomography with the Square Kilometre Array

Authors: Tumelo Mangena, Sultan Hassan, Mario G. Santos

Abstract: Upcoming 21cm surveys with the SKA1-LOW telescope will enable imaging of the neutral hydrogen distribution on cosmological scales in the early Universe. These surveys are expected to generate huge imaging datasets that will encode more information than the power spectrum. This provides an alternative unique way to constrain the reionization history, which might break the degeneracy in the power sp… ▽ More Upcoming 21cm surveys with the SKA1-LOW telescope will enable imaging of the neutral hydrogen distribution on cosmological scales in the early Universe. These surveys are expected to generate huge imaging datasets that will encode more information than the power spectrum. This provides an alternative unique way to constrain the reionization history, which might break the degeneracy in the power spectral analysis. Using Convolutional Neural Networks (CNN), we create a fast estimator of the neutral fraction from the 21cm maps that are produced by our large semi-numerical simulation. Our estimator is able to efficiently recover the neutral fraction ($x_{\rm HI}$) at several redshifts with a high accuracy of 99\% as quantified by the coefficient of determination $R^{2}$. Adding the instrumental effects from the SKA design slightly increases the loss function, but nevertheless we are still able to recover the neutral fraction with a similar high accuracy of 98\%, which is only 1 per cent less. While a weak dependence on redshift is observed, the accuracy increases rapidly with decreasing neutral fraction. This is due to the fact that the instrumental noise increases towards high redshift where the Universe is highly neutral. Our results show the promise of directly using 21cm-tomography to constrain the reionization history in a model independent way, complementing similar efforts, such as those of the optical depth measurements from the Cosmic Microwave Background (CMB) observations by {\it Planck}. △ Less

Submitted 21 April, 2020; v1 submitted 10 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS

A Full Treatment of Peculiar Velocities on the Reionization Light Cone

Authors: Emma Chapman, Mario G. Santos

Abstract: Accurate simulations of the 21-cm signal from the Epoch of Reionization require the generation of maps at specific frequencies given the values of astrophysical and cosmological parameters. The peculiar velocities of the hydrogen atoms producing the 21-cm radiation result in a shift in the observed frequency of the 21-cm radiation and changes the amplitude of the signal itself. This is not an effe… ▽ More Accurate simulations of the 21-cm signal from the Epoch of Reionization require the generation of maps at specific frequencies given the values of astrophysical and cosmological parameters. The peculiar velocities of the hydrogen atoms producing the 21-cm radiation result in a shift in the observed frequency of the 21-cm radiation and changes the amplitude of the signal itself. This is not an effect we can remove but instead needs to be accurately modelled to ensure we infer the correct physical parameters from an observation. We describe the full calculation of the distortion of the 21-cm signal, and propose a new code that integrates the 21-cm intensity along the line of sight for each individual light cone pixel to fully capture the intensity contributions from multiple redshifts. This algorithm naturally deals with the typical divergences found in standard approaches, allowing for large optical depths and 21-cm absorption events at high redshifts. We find the new method results in up to a 60% decrease in power on the largest scales on the sky, and an increase of over 80% on the smallest scales on the sky. We find that the new implementation of the light cone results in a longer tail of bright temperatures in the brightness temperature distribution, as a result of the successful circumventing of a previous cap that had to be implemented to avoid a divergence in the brightness temperature. We conclude that this full treatment of the evolution of the light cone pixel can be an important effect. △ Less

Submitted 21 July, 2020; v1 submitted 27 September, 2019; originally announced September 2019.

Comments: Accepted by MNRAS

The optically-selected 1.4-GHz quasar luminosity function below 1 mJy

Authors: Eliab Malefahlo, Mario G. Santos, Matt J. Jarvis, Sarah V. White, Jonathan T. L. Zwart

Abstract: We present the radio luminosity function (RLF) of optically-selected quasars below 1~mJy, constructed by applying a Bayesian-fitting stacking technique to objects well below the nominal radio flux-density limit. We test the technique using simulated data, confirming that we can reconstruct the RLF over three orders of magnitude below the typical $5σ$ detection threshold. We apply our method to 1.4… ▽ More We present the radio luminosity function (RLF) of optically-selected quasars below 1~mJy, constructed by applying a Bayesian-fitting stacking technique to objects well below the nominal radio flux-density limit. We test the technique using simulated data, confirming that we can reconstruct the RLF over three orders of magnitude below the typical $5σ$ detection threshold. We apply our method to 1.4-GHz flux-densities from the Faint Images of the Radio Sky at Twenty-cm survey (FIRST), extracted at the positions of optical quasars from the Sloan Digital Sky Survey (SDSS) over seven redshift bins up to $z=2.15$ {and measure the RLF down to two orders of magnitude below the FIRST detection threshold}. In the lowest redshift bin ($0.2<z<0.45$), we find that our measured RLF agrees well with deeper data from the literature. The RLF for the radio-loud quasars flattens below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 25.5$ and becomes steeper again below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 24.8$, where radio-quiet quasars start to emerge. The radio luminosity where radio-quiet quasars emerge coincides with the luminosity where star-forming galaxies are expected to start to dominate the radio source counts. This implies that there could be a significant contribution from star formation in the host galaxies, but additional data is required to investigate this further. The higher-redshift bins show a similar behaviour as for the lowest-$z$ bin, implying that the same physical process may be responsible. △ Less

Submitted 16 January, 2020; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: 20 pages, 14 figures, 5 tables

On the contamination of the global 21~cm signal from polarized foregrounds

Authors: Marta Spinelli, Gianni Bernardi, Mario G. Santos

Abstract: Global (i.e. sky-averaged) $21$~cm signal experiments can measure the evolution of the universe from the Cosmic Dawn to the Epoch of Reionization. These measurements are challenged by the presence of bright foreground emission that can be separated from the cosmological signal if its spectrum is smooth. This assumption fails in the case of single polarization antennas as they measure linearly pola… ▽ More Global (i.e. sky-averaged) $21$~cm signal experiments can measure the evolution of the universe from the Cosmic Dawn to the Epoch of Reionization. These measurements are challenged by the presence of bright foreground emission that can be separated from the cosmological signal if its spectrum is smooth. This assumption fails in the case of single polarization antennas as they measure linearly polarized foreground emission - which is inevitably Faraday rotated through the interstellar medium. We investigate the impact of Galactic polarized foregrounds on the extraction of the global 21~cm signal through realistic sky and dipole simulations both in a low frequency band from $50$ to $100$~MHz, where a 21~cm absorption profile is expected, and in a higher frequency band ($100-200$~MHz). We find that the presence of a polarized contaminant with complex frequency structure can bias the amplitude and the shape of the reconstructed signal parameters in both bands. We investigate if polarized foregrounds can explain the unexpected $21$~cm Cosmic Dawn signal recently reported by the EDGES collaboration. We find that unaccounted polarized foreground contamination can produce an enhanced and distorted $21$~cm absorption trough similar to the anomalous profile reported by Bowman et al. (2018), and whose amplitude is in mild tension with the assumed input Gaussian profile (at $\sim 1.5 σ$ level). Moreover, we note that, under the hypothesis of contamination from polarized foreground, the amplitude of the reconstructed EDGES signal can be overestimated by around $30\%$, mitigating the requirement for an explanation based on exotic physics. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Comments: 9 pages, 10 figures, comments welcome

Measuring the HI mass function below the detection threshold

Authors: Hengxing Pan, Matt J. Jarvis, James R. Allison, Ian Heywood, Mario G. Santos, Natasha Maddox, Bradley S. Frank, Xi Kang

Abstract: We present a Bayesian Stacking technique to directly measure the HI mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several sky areas given an assumed HIMF described by a Schechter function and simulate the HI emission lines with different levels of background noise to test the technique. We use Mul… ▽ More We present a Bayesian Stacking technique to directly measure the HI mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several sky areas given an assumed HIMF described by a Schechter function and simulate the HI emission lines with different levels of background noise to test the technique. We use Multinest to constrain the parameters of the HIMF in a broad redshift bin, demonstrating that the HIMF can be accurately reconstructed, using the simulated spectral cube far below the HI mass limit determined by the $5σ$ flux-density limit, i.e. down to $M_{\rm HI} = 10^{7.5}$ M$_{\odot}$ over the redshift range $0 < z < 0.55$ for this particular simulation, with a noise level similar to that expected for the MIGHTEE survey. We also find that the constraints on the parameters of the Schechter function, $φ_{\star}$, $M_\star$ and $α$ can be reliably fit, becoming tighter as the background noise decreases as expected, although the constraints on the redshift evolution are not significantly affected. All the parameters become better constrained as the survey area increases. In summary, we provide an optimal method for estimating the HI mass at cosmological distances that allows us to constrain the HI mass function below the detection threshold in forthcoming HI surveys. This study is a first step towards the measurement of the HIMF at high ($z>0.1$) redshifts. △ Less

Submitted 25 October, 2019; v1 submitted 20 July, 2019; originally announced July 2019.

Comments: 17 pages, 11 figures, revised manuscript accepted for publication in MNRAS

Prospects for cosmic magnification measurements using HI intensity mapping

Authors: Amadeus Witzemann, Alkistis Pourtsidou, Mario G. Santos

Abstract: We investigate the prospects of measuring the cosmic magnification effect by cross-correlating neutral hydrogen intensity mapping (HI IM) maps with background optical galaxies. We forecast the signal-to-noise ratio for HI IM data from SKA1-MID and HIRAX, combined with LSST photometric galaxy samples. We find that, thanks to their different resolutions, SKA1-MID and HIRAX are highly complementary i… ▽ More We investigate the prospects of measuring the cosmic magnification effect by cross-correlating neutral hydrogen intensity mapping (HI IM) maps with background optical galaxies. We forecast the signal-to-noise ratio for HI IM data from SKA1-MID and HIRAX, combined with LSST photometric galaxy samples. We find that, thanks to their different resolutions, SKA1-MID and HIRAX are highly complementary in such an analysis. We predict that SKA1-MID can achieve a detection with a signal-to-noise ratio of $\sim 15$ on a multipole range of $\ell \lesssim 200$, while HIRAX can reach a signal-to-noise ratio of $\sim 50$ on $200 < \ell < 2000$. We conclude that measurements of the cosmic magnification signal will be possible on a wide redshift range with foreground HI intensity maps up to $z \lesssim 2$, while optimal results are obtained when $0.6 \lesssim z \lesssim 1.3$. Finally, we perform a signal to noise analysis that shows how these measurements can constrain the HI parameters across a wide redshift range. △ Less

Submitted 3 June, 2020; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: Accepted version. 9 pages, 7 figures, 2 tables