-
Impact of extragalactic point sources on the low-frequency sky spectrum and cosmic dawn global 21-cm measurements
Authors:
Shikhar Mittal,
Girish Kulkarni,
Dominic Anstey,
Eloy de Lera Acedo
Abstract:
Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements o…
▽ More
Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements of the luminosity function and the angular correlation function of extragalactic point sources, we develop a model for the contribution of these sources to the low-frequency sky spectrum. Using this model, we find that the contribution of sources with flux density $>10^{-6}\,$Jy to the sky-averaged spectrum is smooth and of the order of a few kelvins at 50 - $200\,$MHz. We combine this model with measurements of the galactic foreground spectrum and convolve the result with the beam of the conical log-spiral antenna planned as part of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) project. We find that the contribution of point sources to resultant spectrum is $\sim0.4\%$ of the total foregrounds, but still larger by at least an order of magnitude than the standard predictions for the cosmological 21-cm signal. As a result, not accounting for the point-source contribution leads to a systematic bias in 21-cm signal recovery. We show, however, that in the REACH case, this reconstruction bias can be removed by modelling the point-source contribution as a power law with a running spectral index. We make our code publicly available as a Python package labelled epspy.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
The Spectrometer Development of CosmoCube, Lunar Orbiting Satellite to Detect 21-cm Hydrogen Signal from Cosmic Dark Ages
Authors:
Kaan Artuc,
Eloy de Lera Acedo
Abstract:
The cosmic Dark Ages represent a pivotal epoch in the evolution of the Universe, marked by the emergence of the first cosmic structures under the influence of dark matter. The 21-cm hydrogen line, emanating from the hyperfine transition of neutral hydrogen, serves as a critical probe into this era. We describe the development and implementation of the spectrometer for CosmoCube, a novel lunar orbi…
▽ More
The cosmic Dark Ages represent a pivotal epoch in the evolution of the Universe, marked by the emergence of the first cosmic structures under the influence of dark matter. The 21-cm hydrogen line, emanating from the hyperfine transition of neutral hydrogen, serves as a critical probe into this era. We describe the development and implementation of the spectrometer for CosmoCube, a novel lunar orbiting CubeSat designed to detect the redshifted 21-cm signal within the redshift range of 13 to 150. Our instrumentation utilizes a Xilinx RFSoC, which integrates both Analog-to-Digital Converters (ADCs) and Digital-to-Analog Converters (DACs), tailored for the spectrometer component of the radiometer. This system is characterized by a 4096 FFT length at 62.5 kHz steps using a Polyphase Filter Bank (PFB), achieving an average Effective Number of Bits (ENOB) of 11.5 bits throughout the frequency of interest, from 10 MHz to 100 MHz. The spectrometer design is further refined through loopback tests involving both DAC and ADC of the RFSoC, with DAC outputs varying between high (+1 dBm) and low (-3 dBm) power modes to characterize system performance. The power consumption was optimized to 5.45 W using three ADCs and one DAC for the radiometer. Additionally, the stability of the ADC noise floor was investigated in a thermal chamber with environmental temperatures ranging from 5°C to 40°C. A consistent noise floor of approximately -152.5 dBFS/Hz was measured, with a variation of $\pm$0.2 dB, ensuring robust performance under varying thermal conditions.
△ Less
Submitted 14 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.
-
Calibration Error in 21-centimeter Global Spectrum Experiments
Authors:
Shijie Sun,
Eloy de Lera Acedo,
Fengquan Wu,
Bin Yue,
Jiacong Zhu,
Xuelei Chen
Abstract:
The redshifted 21 cm line signal is a powerful probe of the cosmic dawn and the epoch of reionization. The global spectrum can potentially be detected with a single antenna and spectrometer. However, this measurement requires an extremely accurate calibration of the instrument to facilitate the separation of the 21 cm signal from the much brighter foregrounds and possible variations in the instrum…
▽ More
The redshifted 21 cm line signal is a powerful probe of the cosmic dawn and the epoch of reionization. The global spectrum can potentially be detected with a single antenna and spectrometer. However, this measurement requires an extremely accurate calibration of the instrument to facilitate the separation of the 21 cm signal from the much brighter foregrounds and possible variations in the instrument response. Understanding how the measurement errors propagate in a realistic instrument system and affect system calibration is the focus of this work. We simulate a 21 cm global spectrum observation based on the noise wave calibration scheme. We focus on how measurement errors in reflection coefficients affect the noise temperature and how typical errors impact the recovery of the 21 cm signal, especially in the frequency domain. Results show that for our example set up, a typical vector network analyzer (VNA) measurement error in the magnitude of the reflection coefficients of the antenna, receiver, and open cable, which are 0.001, 0.001, and 0.002 (linear), respectively, would result in a 200 mK deviation on the detected signal, and a typical measurement error of 0.48 degree, 0.78 degree, or 0.15 degree in the respective phases would cause a 40 mK deviation. The VNA measurement error can greatly affect the result of a 21 cm global spectrum experiment using this calibration technique, and such a feature could be mistaken for or be combined with the 21 cm signal
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
Authors:
Anchal Saxena,
P. Daniel Meerburg,
Christoph Weniger,
Eloy de Lera Acedo,
Will Handley
Abstract:
The redshifted 21-cm signal from the Cosmic Dawn and Epoch of Reionization carries invaluable information about the cosmology and astrophysics of the early Universe. Analyzing the data from a sky-averaged 21-cm signal experiment typically involves navigating through an intricate parameter space to accurately address various factors such as foregrounds, beam uncertainties, ionospheric distortions,…
▽ More
The redshifted 21-cm signal from the Cosmic Dawn and Epoch of Reionization carries invaluable information about the cosmology and astrophysics of the early Universe. Analyzing the data from a sky-averaged 21-cm signal experiment typically involves navigating through an intricate parameter space to accurately address various factors such as foregrounds, beam uncertainties, ionospheric distortions, and receiver noise for the search of the cosmological 21-cm signal. The traditional likelihood-based sampling methods for modeling these effects could become computationally demanding for such highly complex models, which makes it infeasible to include physically motivated 21-cm signal models in the analysis. Moreover, the inference with these traditional methods is driven by the assumed functional form of the likelihood function. This work demonstrates how Simulation-Based Inference through Truncated Marginal Neural Ratio Estimation (TMNRE) can naturally handle these issues at a significantly reduced computational cost than the likelihood-based methods. We estimate the posterior distribution on our model parameters with TMNRE for simulated mock observations, composed of beam-weighted foregrounds, physically motivated 21-cm signal, and radiometric noise. We find that maximizing the information content by simultaneously analyzing the data from multiple time slices and antennas significantly improves the parameter constraints and leads to a better exploration of the cosmological signal. We discuss the application of TMNRE for the current configuration of the REACH experiment and demonstrate how it can be utilized to explore potential avenues for REACH. The method presented here can be easily generalized for any sky-averaged 21-cm signal experiment.
△ Less
Submitted 21 March, 2024;
originally announced March 2024.
-
High-dimensional inference of radio interferometer beam patterns I: Parametric model of the HERA beams
Authors:
Michael J. Wilensky,
Jacob Burba,
Philip Bull,
Hugh Garsden,
Katrine A. Glasscock,
Nicolas Fagnoni,
Eloy de Lera Acedo,
David R. DeBoer,
Nima Razavi-Ghods
Abstract:
Accurate modelling of the primary beam is an important but difficult task in radio astronomy. For high dynamic range problems such as 21cm intensity mapping, small modelling errors in the sidelobes and spectral structure of the beams can translate into significant systematic errors. Realistic beams exhibit complex spatial and spectral structure, presenting a major challenge for beam measurement an…
▽ More
Accurate modelling of the primary beam is an important but difficult task in radio astronomy. For high dynamic range problems such as 21cm intensity mapping, small modelling errors in the sidelobes and spectral structure of the beams can translate into significant systematic errors. Realistic beams exhibit complex spatial and spectral structure, presenting a major challenge for beam measurement and calibration methods. In this paper series, we present a Bayesian framework to infer per-element beam patterns from the interferometric visibilities for large arrays with complex beam structure, assuming a particular (but potentially uncertain) sky model and calibration solution. In this first paper, we develop a compact basis for the beam so that the Bayesian computation is tractable with high-dimensional sampling methods. We use the Hydrogen Epoch of Reionization Array (HERA) as an example, verifying that the basis is capable of describing its single-element E-field beam (i.e. without considering array effects like mutual coupling) with a relatively small number of coefficients. We find that 32 coefficients per feed, incident polarization, and frequency, are sufficient to give percent-level and $\sim$10\% errors in the mainlobe and sidelobes respectively for the current HERA Vivaldi feeds, improving to $\sim 0.1\%$ and $\sim 1\%$ for 128 coefficients.
△ Less
Submitted 24 July, 2024; v1 submitted 20 March, 2024;
originally announced March 2024.
-
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.
-
Understanding spectral artefacts in SKA-LOW 21-cm cosmology experiments: the impact of cable reflections
Authors:
Oscar Sage David O'Hara,
Fred Dulwich,
Eloy de Lera Acedo,
Jiten Dhandha,
Thomas Gessey-Jones,
Dominic Anstey,
Anastasia Fialkov
Abstract:
The Cosmic Dawn marks the first star formations and preceded the Epoch-of-Reionization, when the Universe underwent a fundamental transformation propelled by the radiation from these first stars and galaxies. Interferometric 21-cm experiments aim to probe redshifted neutral hydrogen signals from these periods, constraining the conditions of the early Universe. The SKA-LOW instrument of the Square…
▽ More
The Cosmic Dawn marks the first star formations and preceded the Epoch-of-Reionization, when the Universe underwent a fundamental transformation propelled by the radiation from these first stars and galaxies. Interferometric 21-cm experiments aim to probe redshifted neutral hydrogen signals from these periods, constraining the conditions of the early Universe. The SKA-LOW instrument of the Square Kilometre Array telescope is envisaged to be the largest and most sensitive radio telescope at m and cm wavelengths. The latest Aperture Array Verification Systems feature 7m coaxial transmission lines connecting the Low Noise Amplifiers to optical transmitters at the front of the analogue-receiving chain. An impedance mismatch between these components results in a partially reflected electromagnetic signal, which introduces chromatic aberrations in the instrument bandpass. This causes power from the foreground signals to appear at higher delays, potentially contaminating the EoR window, a region at which the 21-cm signal should be detectable. We present an end-to-end simulation pipeline for SKA-LOW using a composite sky model combining radio foregrounds from The GLEAM Survey, Haslam $408$MHz, and a $1.5$cGpc 21-cm brightness temperature cube generated with the 21cmSPACE simulator. Iterating a parametric approach, we derive a model for the scattering parameters of a coaxial transmission line in terms of its specifications and bulk material properties. Assuming identical cables of length $\leq 15.0$m with impedance mismatch $\leq 10Ω$ confines the reflection to k-modes below the EoR window. However, we demonstrate that even a $0.1$% length tolerance introduces contamination with an RMSE of $\sim 10$% across all accessible k-modes.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
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.
-
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.
-
On the Constraints on Superconducting Cosmic Strings from 21-cm Cosmology
Authors:
T. Gessey-Jones,
S. Pochinda,
H. T. J. Bevins,
A. Fialkov,
W. J. Handley,
E. de Lera Acedo,
S. Singh,
R. Barkana
Abstract:
Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainti…
▽ More
Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainties in foregrounds and high redshift astrophysics. We include the latest publicly available 21-cm power spectrum upper limits from HERA, 21-cm global signal data from SARAS 3, and the synergistic probe of the unresolved X-ray background in our final analysis. This paper thus constitutes the first attempt to use 21-cm power spectrum data to probe cosmic strings. In contrast to previous works, we find no strong constraints can be placed on superconducting cosmic strings from current 21-cm measurements. This is because of uncertainties in the X-ray emission efficiency of the first galaxies, with X-ray emissivities greater than $3 \times 10^{40}$erg s$^{-1}$ M$_{\odot}^{-1}$ yr able to mask the presence of cosmic strings in the 21-cm signal. We conclude by discussing the prospects for future constraints from definitive 21-cm signal measurements and argue that the recently proposed soft photon heating should be cause for optimism due to its potential to break degeneracies that would have otherwise made the signatures of cosmic strings difficult to distinguish from those of astrophysical origin.
△ Less
Submitted 15 February, 2024; v1 submitted 14 December, 2023;
originally announced December 2023.
-
Constraining the properties of Population III galaxies with multi-wavelength observations
Authors:
S. Pochinda,
T. Gessey-Jones,
H. T. J. Bevins,
A. Fialkov,
S. Heimersheim,
I. Abril-Cabezas,
E. de Lera Acedo,
S. Singh,
S. Sikder,
R. Barkana
Abstract:
The early Universe, spanning 400,000 to 400 million years after the Big Bang ($z\approx1100-11$), has been left largely unexplored as the light from luminous objects is too faint to be observed directly. While new experiments are pushing the redshift limit of direct observations, measurements in the low-frequency radio band promise to probe early star and black hole formation via observations of t…
▽ More
The early Universe, spanning 400,000 to 400 million years after the Big Bang ($z\approx1100-11$), has been left largely unexplored as the light from luminous objects is too faint to be observed directly. While new experiments are pushing the redshift limit of direct observations, measurements in the low-frequency radio band promise to probe early star and black hole formation via observations of the hydrogen 21-cm line. In this work we explore synergies between 21-cm data from the HERA and SARAS 3 experiments and observations of the unresolved radio and X-ray backgrounds using multi-wavelength Bayesian analysis. We use the combined data set to constrain properties of Population II and Population III stars as well as early X-ray and radio sources. The joint fit reveals a 68 percentile disfavouring of Population III star formation efficiencies $\gtrsim5.7\%$. We also show how the 21-cm and the X-ray background data synergistically constrain opposite ends of the X-ray efficiency prior distribution to produce a peak in the 1D posterior of the X-ray luminosity per star formation rate. We find (at 68\% confidence) that early galaxies were likely 0.3 to 318 times as X-ray efficient as present-day starburst galaxies. We also show that the functional posteriors from our joint fit rule out global 21-cm signals deeper than $\lesssim-203\ \mathrm{mK}$ and power spectrum amplitudes at $k=0.34\ h\mathrm{Mpc^{-1}}$ greater than $Δ_{21}^2 \gtrsim 946\ \mathrm{mK}^2$ with $3σ$ confidence.
△ Less
Submitted 1 May, 2024; v1 submitted 13 December, 2023;
originally announced December 2023.
-
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.
-
FlexKnot as a Generalised Model of the Sky-averaged 21-cm Signal at z ~ 6-30 in the Presence of Systematics
Authors:
Emma Shen,
Dominic Anstey,
Eloy de Lera Acedo,
Anastasia Fialkov
Abstract:
Global 21-cm experiments are built to study the evolution of the Universe between the cosmic dawn and the epoch of reionisation. FlexKnot is a function parameterised by freely moving knots stringed together by splines. Adopting the FlexKnot function as the signal model has the potential to separate the global 21-cm signal from the foregrounds and systematics while being capable of recovering the c…
▽ More
Global 21-cm experiments are built to study the evolution of the Universe between the cosmic dawn and the epoch of reionisation. FlexKnot is a function parameterised by freely moving knots stringed together by splines. Adopting the FlexKnot function as the signal model has the potential to separate the global 21-cm signal from the foregrounds and systematics while being capable of recovering the crucial features given by theoretical predictions. In this paper, we implement the FlexKnot method by integrating twice over a function of freely moving knots interpolated linearly. The function is also constrained at the lower frequencies corresponding to the dark ages by theoretical values. The FlexKnot model is tested in the framework of the realistic data analysis pipeline of the REACH global signal experiment using simulated antenna temperature data. We demonstrate that the FlexKnot model performs better than existing signal models, e.g. the Gaussian signal model, at reconstructing the shape of the true signals present in the simulated REACH data, especially for injected signals with complex structures. The capabilities of the FlexKnot signal model is also tested by introducing various systematics and simulated global signals of different types. These tests show that four to five knots are sufficient to recover the general shape of most realistic injected signals, with or without sinusoidal systematics. We show that true signals whose absorption trough of amplitude between 120 to 450 mK can be well recovered with systematics up to about 50 mK.
△ Less
Submitted 24 November, 2023;
originally announced November 2023.
-
Direct Optimal Mapping Image Power Spectrum and its Window Functions
Authors:
Zhilei Xu,
Honggeun Kim,
Jacqueline N. Hewitt,
Kai-Feng Chen,
Nicholas S. Kern,
Eleanor Rath,
Ruby Byrne,
Adélie Gorce,
Robert Pascua,
Zachary E. Martinot,
Joshua S. Dillon,
Bryna J. Hazelton,
Adrian Liu,
Miguel F. Morales,
Zara 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
, et al. (57 additional authors not shown)
Abstract:
The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based…
▽ More
The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show $<10^{-11}$ of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR.
△ Less
Submitted 5 July, 2024; v1 submitted 17 November, 2023;
originally announced November 2023.
-
The effects of the antenna power pattern uncertainty within a global 21 cm experiment
Authors:
John Cumner,
Carla Pieterse,
Dirk De Villiers,
Eloy de Lera Acedo
Abstract:
Experimental 21 cm cosmology aims to detect the formation of the first stars during the cosmic dawn and the subsequent epoch of reionization by utilizing the 21 cm hydrogen line transition. While several experiments have published results that begin to constrain the shape of this signal, a definitive detection has yet to be achieved. In this paper, we investigate the influence of uncertain antenna…
▽ More
Experimental 21 cm cosmology aims to detect the formation of the first stars during the cosmic dawn and the subsequent epoch of reionization by utilizing the 21 cm hydrogen line transition. While several experiments have published results that begin to constrain the shape of this signal, a definitive detection has yet to be achieved. In this paper, we investigate the influence of uncertain antenna-sky interactions on the possibility of detecting the signal. This paper aims to define the level of accuracy to which a simulated antenna beam pattern is required to agree with the actual observing beam pattern of the antenna to allow for a confident detection of the global 21 cm signal. By utilising singular value decomposition, we construct a set of antenna power patterns that incorporate minor, physically motivated variations. We take the absolute mean averaged difference between the original beam and the perturbed beam averaged over frequency ($ΔD$) to quantifying this difference, identifying the correlation between $ΔD$ and antenna temperature. To analyse the impact of $ΔD$ on making a confident detection, we utilize the REACH Bayesian analysis pipeline and compare the Bayesian evidence $\log \mathcal{Z}$ and root-mean-square error for antenna beams of different $ΔD$ values. Our calculations suggest that achieving an agreement between the original and perturbed antenna power pattern with $ΔD$ better than -35 dB is necessary for confident detection of the global 21 cm signal. Furthermore, we discuss potential methods to achieve the required high level of accuracy within a global 21~cm experiment.
△ Less
Submitted 13 November, 2023;
originally announced November 2023.
-
A Bayesian Method to Mitigate the Effects of Unmodelled Time-Varying Systematics for 21-cm Cosmology Experiments
Authors:
Christian J. Kirkham,
Dominic J. Anstey,
Eloy de Lera Acedo
Abstract:
Radio observations of the neutral hydrogen signal from the Cosmic Dawn and Epoch of Reionisation have helped to provide constraints on the properties of the first stars and galaxies. Since this global 21-cm cosmological signal from the Cosmic Dawn is effectively constant on observing timescales and since effects resulting from systematics will vary with time, the effects of these systematics can b…
▽ More
Radio observations of the neutral hydrogen signal from the Cosmic Dawn and Epoch of Reionisation have helped to provide constraints on the properties of the first stars and galaxies. Since this global 21-cm cosmological signal from the Cosmic Dawn is effectively constant on observing timescales and since effects resulting from systematics will vary with time, the effects of these systematics can be mitigated without the need for a model of the systematic. We present a method to account for unmodelled time-varying systematics in 21-cm radio cosmology experiments using a squared-exponential Gaussian process kernel to account for correlations between time bins in a fully Bayesian way. We find by varying the model parameters of a simulated systematic that the Gaussian process method improves our ability to recover the signal parameters by widening the posterior in the presence of a systematic and reducing the bias in the mean fit parameters. When varying the amplitude of a model sinusoidal systematic between 0.25 and 2.00 times the 21-cm signal amplitude and the period between 0.5 and 4.0 times the signal width, we find on average a 5% improvement in the root mean squared error of the fitted signal. We can use the fitted Gaussian process hyperparameters to identify the presence of a systematic in the data, demonstrating the method's utility as a diagnostic tool. Furthermore, we can use Gaussian process regression to calculate a mean fit to the residuals over time, providing a basis for producing a model of the time-varying systematic.
△ Less
Submitted 27 October, 2023;
originally announced October 2023.
-
The Impact of Beam Variations on Power Spectrum Estimation for 21 cm Cosmology II: Mitigation of Foreground Systematics for HERA
Authors:
Honggeun Kim,
Nicholas S. Kern,
Jacqueline N. Hewitt,
Bang D. Nhan,
Joshua S. Dillon,
Eloy de Lera Acedo,
Scott B. C. Dynes,
Nivedita Mahesh,
Nicolas Fagnoni,
David R. DeBoer
Abstract:
One key challenge in detecting 21 cm cosmological signal at z > 6 is to separate the cosmological signal from foreground emission. This can be studied in a power spectrum space where the foreground is confined to low delay modes whereas the cosmological signal can spread out to high delay modes. When there is a calibration error, however, chromaticity of gain errors propagates to the power spectru…
▽ More
One key challenge in detecting 21 cm cosmological signal at z > 6 is to separate the cosmological signal from foreground emission. This can be studied in a power spectrum space where the foreground is confined to low delay modes whereas the cosmological signal can spread out to high delay modes. When there is a calibration error, however, chromaticity of gain errors propagates to the power spectrum estimate and contaminates the modes for cosmological detection. The Hydrogen Epoch of Reionization Array (HERA) employs a high-precision calibration scheme using redundancy in measurements. In this study, we focus on the gain errors induced by nonredundancies arising from feed offset relative to the HERA's 14 meter parabolic dish element, and investigate how to mitigate the chromatic gain errors using three different methods: restricting baseline lengths for calibration, smoothing the antenna gains, and applying a temporal filter prior to calibration. With 2 cm/2 degree perturbations for translation/tilting motions, a level achievable under normal HERA operating conditions, the combination of the baseline cut and temporal filtering indicates that the spurious gain feature due to nonredundancies is significantly reduced, and the power spectrum recovers the clean foreground-free region. We found that the mitigation technique works even for large feed motions but in order to keep a stable calibration process, the feed positions need to be constrained to 2 cm for translation motions and 2 degree for tilting offset relative to the dish's vertex.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Modelling a Hot Horizon in Global 21 cm Experimental Foregrounds
Authors:
Joe H. N. Pattison,
Dominic J. Anstey,
Eloy de Lera Acedo
Abstract:
The 21 cm signal from cosmic hydrogen is one of the most propitious probes of the early Universe. The detection of this signal would reveal key information about the first stars, the nature of dark matter, and early structure formation. We explore the impact of an emissive and reflective, or `hot', horizon on the recovery of this signal for global 21 cm experiments. It is demonstrated that using p…
▽ More
The 21 cm signal from cosmic hydrogen is one of the most propitious probes of the early Universe. The detection of this signal would reveal key information about the first stars, the nature of dark matter, and early structure formation. We explore the impact of an emissive and reflective, or `hot', horizon on the recovery of this signal for global 21 cm experiments. It is demonstrated that using physically motivated foreground models to recover the sky-averaged 21 cm signal one must accurately describe the horizon around the radiometer. We show that not accounting for the horizon will lead to a signal recovery with residuals an order of magnitude larger than the injected signal, with a log Bayesian evidence of almost 1600 lower than when one does account for the horizon. It is shown that signal recovery is sensitive to incorrect values of soil temperature and reflection coefficient in describing the horizon, with even a 10% error in reflectance causing twofold increases in the RMSE of a given fit. We also show these parameters may be fitted using Bayesian inference to mitigate for these issues without overfitting and mischaracterising a non-detection. We further demonstrate that signal recovery is sensitive to errors in measurements of the horizon projection onto the sky, but fitting for soil temperature and reflection coefficients with priors that extend beyond physical expectation can resolve these problems. We show that using an expanded prior range can reliably recover the signal even when the height of the horizon is mismeasured by up to 20%, decreasing the RMSE from the model that does not perform this fitting by a factor of 9.
△ Less
Submitted 6 July, 2023;
originally announced July 2023.
-
Receiver design for the REACH global 21-cm signal experiment
Authors:
Nima Razavi-Ghods,
Ian L. V. Roque,
Steven H. Carey,
John A. Ely,
Will Handley,
Alessio Magro,
Riccardo Chiello,
Tian Huang,
P. Alexander,
D. Anstey,
G. Bernardi,
H. T. J. Bevins,
J. Cavillot,
W. Croukamp,
J. Cumner,
E. de Lera Acedo,
D. I. L. de Villiers,
A. Fialkov,
T. Gessey-Jones,
Q. Gueuning,
A. T. Josaitis,
G. Kulkarni,
S. A. K. Leeney,
R. Maiolino,
P. D. Meerburg
, et al. (13 additional authors not shown)
Abstract:
We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrai…
▽ More
We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrained least squares approach demonstrating a calibration RMSE of 80 mK for five hours of integration on a custom-made source with comparable impedance to that of the antenna used in the field. This paper therefore documents the state of the calibrator and data analysis in December 2022 in Cambridge before shipping to South Africa.
△ Less
Submitted 14 July, 2023; v1 submitted 30 June, 2023;
originally announced July 2023.
-
Signatures of Cosmic Ray Heating in 21-cm Observables
Authors:
T. Gessey-Jones,
A. Fialkov,
E. de Lera Acedo,
W. J. Handley,
R. Barkana
Abstract:
Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D semi-numerical 21-cm signal simulations and con…
▽ More
Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D semi-numerical 21-cm signal simulations and conduct the first investigations of the signatures it imprints on the 21-cm power spectrum and tomographic maps. We find that cosmic ray heating of the IGM is short-ranged, leading to heating clustered around star-forming sites, and a sharp contrast between heated regions of 21-cm emission and unheated regions of absorption. This contrast results in greater small-scale power for cosmic ray heated scenarios compared to what is found for X-ray heating, thus suggesting a way to test the nature of IGM heating with future 21-cm observations. Finally, we find an unexpectedly rich thermal history in models where cosmic rays can only escape efficiently from low-mass halos, such as in scenarios where these energetic particles originate from population III star supernovae remnants. The interplay of heating and the Lyman-Werner feedback in these models can produce a local peak in the IGM kinetic temperature and, for a limited parameter range, a flattened absorption trough in the global 21-cm signal.
△ Less
Submitted 3 October, 2023; v1 submitted 14 April, 2023;
originally announced April 2023.
-
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.
-
On the use of temporal filtering for mitigating galactic synchrotron calibration bias in 21 cm reionization observations
Authors:
Ntsikelelo Charles,
Nicholas Kern,
Gianni Bernardi,
Landman Bester,
Oleg Smirnov,
Nicolas Fagnoni,
Eloy de Lera Acedo
Abstract:
Precision antenna calibration is required for mitigating the impact of foreground contamination in 21 cm cosmological radio surveys. One widely studied source of error is the effect of missing point sources in the calibration sky model; however, poorly understood diffuse galactic emission also creates a calibration bias that can complicate the clean separation of foregrounds from the 21 cm signal.…
▽ More
Precision antenna calibration is required for mitigating the impact of foreground contamination in 21 cm cosmological radio surveys. One widely studied source of error is the effect of missing point sources in the calibration sky model; however, poorly understood diffuse galactic emission also creates a calibration bias that can complicate the clean separation of foregrounds from the 21 cm signal. In this work, we present a technique for suppressing this bias with temporal filtering of radio interferometric visibilities observed in a drift-scan mode. We demonstrate this technique on mock simulations of the Hydrogen Epoch of Reionization Array (HERA) experiment. Inspecting the recovered calibration solutions, we find that our technique reduces spurious errors by over an order of magnitude. This improved accuracy approaches the required accuracy needed to make a fiducial detection of the 21 cm signal with HERA, but is dependent on a number of external factors that we discuss. We also explore different types of temporal filtering techniques and discuss their relative performance and tradeoffs.
△ Less
Submitted 4 April, 2023; v1 submitted 1 February, 2023;
originally announced February 2023.
-
Joint analysis constraints on the physics of the first galaxies with low frequency radio astronomy data
Authors:
Harry T. J. Bevins,
Stefan Heimersheim,
Irene Abril-Cabezas,
Anastasia Fialkov,
Eloy de Lera Acedo,
William Handley,
Saurabh Singh,
Rennan Barkana
Abstract:
The first billion years of cosmic history remains largely unobserved. We demonstrate, using a novel machine learning technique, how combining upper limits on the spatial fluctuations in the 21-cm signal with observations of the sky-averaged 21-cm signal from neutral hydrogen can improve our understanding of this epoch. By jointly analysing data from SARAS3 (redshift $z\approx15-25$) and limits fro…
▽ More
The first billion years of cosmic history remains largely unobserved. We demonstrate, using a novel machine learning technique, how combining upper limits on the spatial fluctuations in the 21-cm signal with observations of the sky-averaged 21-cm signal from neutral hydrogen can improve our understanding of this epoch. By jointly analysing data from SARAS3 (redshift $z\approx15-25$) and limits from HERA ($z\approx8$ and $10$), we show that such a synergetic analysis provides tighter constraints on the astrophysics of galaxies 200 million years after the Big Bang than can be achieved with the individual data sets. Although our constraints are weak, this is the first time data from a sky-averaged 21-cm experiment and power spectrum experiment have been analysed together. In synergy, the two experiments leave only $64.9^{+0.3}_{-0.1}$% of the explored broad theoretical parameter space to be consistent with the joint data set, in comparison to $92.3^{+0.3}_{-0.1}$% for SARAS3 and $79.0^{+0.5}_{-0.2}$% for HERA alone. We use the joint analysis to constrain star formation efficiency, minimum halo mass for star formation, X-ray luminosity of early emitters and the radio luminosity of early galaxies. The joint analysis disfavours at 68% confidence a combination of galaxies with X-ray emission that is $\lesssim 33$ and radio emission that is $\gtrsim 32$ times as efficient as present day galaxies. We disfavour at $95$% confidence scenarios in which power spectra are $\geq126$ mK$^{2}$ at $z=25$ and the sky-averaged signals are $\leq-277$ mK.
△ Less
Submitted 18 December, 2023; v1 submitted 9 January, 2023;
originally announced January 2023.
-
Sky-averaged 21-cm signal extraction using multiple antennas with an SVD framework: the REACH case
Authors:
Anchal Saxena,
P. Daniel Meerburg,
Eloy de Lera Acedo,
Will Handley,
Léon V. E. Koopmans
Abstract:
In a sky-averaged 21-cm signal experiment, the uncertainty on the extracted signal depends mainly on the covariance between the foreground and 21-cm signal models. In this paper, we construct these models using the modes of variation obtained from the Singular Value Decomposition of a set of simulated foreground and 21-cm signals. We present a strategy to reduce this overlap between the 21-cm and…
▽ More
In a sky-averaged 21-cm signal experiment, the uncertainty on the extracted signal depends mainly on the covariance between the foreground and 21-cm signal models. In this paper, we construct these models using the modes of variation obtained from the Singular Value Decomposition of a set of simulated foreground and 21-cm signals. We present a strategy to reduce this overlap between the 21-cm and foreground modes by simultaneously fitting the spectra from multiple different antennas, which can be used in combination with the method of utilizing the time dependence of foregrounds while fitting multiple drift scan spectra. To demonstrate this idea, we consider two different foreground models (i) a simple foreground model, where we assume a constant spectral index over the sky, and (ii) a more realistic foreground model, with a spatial variation of the spectral index. For the simple foreground model, with just a single antenna design, we are able to extract the signal with good accuracy if we simultaneously fit the data from multiple time slices. The 21-cm signal extraction is further improved when we simultaneously fit the data from different antennas as well. This improvement becomes more pronounced while using the more realistic mock observations generated from the detailed foreground model. We find that even if we fit multiple time slices, the recovered signal is biased and inaccurate for a single antenna. However, simultaneously fitting the data from different antennas reduces the bias and the uncertainty by a factor of 2-3 on the extracted 21-cm signal.
△ Less
Submitted 6 April, 2023; v1 submitted 14 December, 2022;
originally announced December 2022.
-
Astrophysical Constraints from the SARAS3 non-detection of the Cosmic Dawn Sky-Averaged 21-cm Signal
Authors:
H. T. J. Bevins,
A. Fialkov,
E. de Lera Acedo,
W. J. Handley,
S. Singh,
R. Subrahmanyan,
R. Barkana
Abstract:
Observations of the redshifted 21-cm line of atomic hydrogen have resulted in several upper limits on the 21-cm power spectrum and a tentative detection of the sky-averaged signal at $z\sim17$. Made with the EDGES Low-Band antenna, this claim was recently disputed by the SARAS3 experiment, which reported a non-detection and is the only available upper limit strong enough to constrain cosmic dawn a…
▽ More
Observations of the redshifted 21-cm line of atomic hydrogen have resulted in several upper limits on the 21-cm power spectrum and a tentative detection of the sky-averaged signal at $z\sim17$. Made with the EDGES Low-Band antenna, this claim was recently disputed by the SARAS3 experiment, which reported a non-detection and is the only available upper limit strong enough to constrain cosmic dawn astrophysics. We use these data to constrain a population of radio-luminous galaxies $\sim 200$ million years after the Big Bang ($z\approx 20$). We find, using Bayesian data analysis, that the data disfavours (at 68% confidence) radio-luminous galaxies in dark matter halos with masses of $4.4\times10^{5}$ M$_\odot \lesssim M \lesssim 1.1\times10^{7}$M$_\odot$ (where $M_\odot$ is the mass of the Sun) at $z = 20$ and galaxies in which $>5$% of the gas is converted into stars. The data disfavour galaxies with radio luminosity per star formation rate of $L_\mathrm{r}/\mathrm{SFR} \gtrsim 1.549 \times 10^{25}$ W Hz$^{-1}$M$_\odot^{-1}$ yr at 150 MHz, a thousand times brighter than today, and, separately, a synchrotron radio background in excess of the CMB by $\gtrsim 6%$ at 1.42 GHz.
△ Less
Submitted 1 December, 2022;
originally announced December 2022.
-
Bayesian approach to radio frequency interference mitigation
Authors:
S. A. K. Leeney,
W. J. Handley,
E. de Lera Acedo
Abstract:
Interfering signals such as Radio Frequency Interference from ubiquitous satellite constellations are becoming an endemic problem in fields involving physical observations of the electromagnetic spectrum. To address this we propose a novel data cleaning methodology. Contamination is simultaneously flagged and managed at the likelihood level. It is modeled in a Bayesian fashion through a piecewise…
▽ More
Interfering signals such as Radio Frequency Interference from ubiquitous satellite constellations are becoming an endemic problem in fields involving physical observations of the electromagnetic spectrum. To address this we propose a novel data cleaning methodology. Contamination is simultaneously flagged and managed at the likelihood level. It is modeled in a Bayesian fashion through a piecewise likelihood that is constrained by a Bernoulli prior distribution. The techniques described in this paper can be implemented with just a few lines of code.
△ Less
Submitted 19 February, 2024; v1 submitted 28 November, 2022;
originally announced November 2022.
-
A General Bayesian Framework to Account for Foreground Map Errors in Global 21-cm Experiments
Authors:
Michael Pagano,
Peter Sims,
Adrian Liu,
Dominic Anstey,
Will Handley,
Eloy De Lera Acedo
Abstract:
Measurement of the global 21-cm signal during Cosmic Dawn (CD) and the Epoch of Reionization (EoR) is made difficult by bright foreground emission which is 2-5 orders of magnitude larger than the expected signal. Fitting for a physics-motivated parametric forward model of the data within a Bayesian framework provides a robust means to separate the signal from the foregrounds, given sufficient info…
▽ More
Measurement of the global 21-cm signal during Cosmic Dawn (CD) and the Epoch of Reionization (EoR) is made difficult by bright foreground emission which is 2-5 orders of magnitude larger than the expected signal. Fitting for a physics-motivated parametric forward model of the data within a Bayesian framework provides a robust means to separate the signal from the foregrounds, given sufficient information about the instrument and sky. It has previously been demonstrated that, within such a modelling framework, a foreground model of sufficient fidelity can be generated by dividing the sky into $N$ regions and scaling a base map assuming a distinct uniform spectral index in each region. Using the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) as our fiducial instrument, we show that, if unaccounted-for, amplitude errors in low-frequency radio maps used for our base map model will prevent recovery of the 21-cm signal within this framework, and that the level of bias in the recovered 21-cm signal is proportional to the amplitude and the correlation length of the base-map errors in the region. We introduce an updated foreground model that is capable of accounting for these measurement errors by fitting for a monopole offset and a set of spatially-dependent scale factors describing the ratio of the true and model sky temperatures, with the size of the set determined by Bayesian evidence-based model comparison. We show that our model is flexible enough to account for multiple foreground error scenarios allowing the 21-cm sky-averaged signal to be detected without bias from simulated observations with a smooth conical log spiral antenna.
△ Less
Submitted 1 November, 2023; v1 submitted 18 November, 2022;
originally announced November 2022.
-
The Impact of Beam Variations on Power Spectrum Estimation for 21-cm Cosmology I: Simulations of Foreground Contamination for HERA
Authors:
Honggeun Kim,
Bang D. Nhan,
Jacqueline N. Hewitt,
Nicholas S. Kern,
Joshua S. Dillon,
Eloy de Lera Acedo,
Scott Dynes,
Nivedita Mahesh,
Nicolas Fagnoni,
David R. DeBoer
Abstract:
Detecting cosmological signals from the Epoch of Reionization (EoR) requires high-precision calibration to isolate the cosmological signals from foreground emission. In radio interferometery, perturbed primary beams of antenna elements can disrupt the precise calibration, which results in contaminating the foreground-free region, or the EoR window, in the cylindrically averaged power spectrum. For…
▽ More
Detecting cosmological signals from the Epoch of Reionization (EoR) requires high-precision calibration to isolate the cosmological signals from foreground emission. In radio interferometery, perturbed primary beams of antenna elements can disrupt the precise calibration, which results in contaminating the foreground-free region, or the EoR window, in the cylindrically averaged power spectrum. For Hydrogen Epoch of Reionization Array (HERA), we simulate and characterize the perturbed primary beams induced by feed motions such as axial, lateral, and tilting motions, above the 14-meter dish. To understand the effect of the perturbed beams, visibility measurements are modeled with two different foreground components, point sources and diffuse sources, and we find different feed motions present a different reaction to each type of sky source. HERA's redundant-baseline calibration in the presence of non-redundant antenna beams due to feed motions introduces chromatic errors in gain solutions, which produces foreground power leakage into the EoR window. The observed leakage from vertical feed motions comes predominately from point sources around zenith. Furthermore, the observed leakage from horizontal and tilting feed motion comes predominately from the diffuse components near the horizon. Mitigation of chromatic gain errors will be necessary for robust detection of the EoR signals with minimal foreground bias, and this will be discussed in the subsequent paper.
△ Less
Submitted 28 October, 2022;
originally announced October 2022.
-
Characterization Of Inpaint Residuals In Interferometric Measurements of the Epoch Of Reionization
Authors:
Michael Pagano,
Jing Liu,
Adrian Liu,
Nicholas S. Kern,
Aaron Ewall-Wice,
Philip Bull,
Robert Pascua,
Siamak Ravanbakhsh,
Zara 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,
Jacob Burba,
Steven Carey,
Chris L. Carilli,
Carina Cheng,
David R. DeBoer
, et al. (53 additional authors not shown)
Abstract:
Radio Frequency Interference (RFI) is one of the systematic challenges preventing 21cm interferometric instruments from detecting the Epoch of Reionization. To mitigate the effects of RFI on data analysis pipelines, numerous inpaint techniques have been developed to restore RFI corrupted data. We examine the qualitative and quantitative errors introduced into the visibilities and power spectrum du…
▽ More
Radio Frequency Interference (RFI) is one of the systematic challenges preventing 21cm interferometric instruments from detecting the Epoch of Reionization. To mitigate the effects of RFI on data analysis pipelines, numerous inpaint techniques have been developed to restore RFI corrupted data. We examine the qualitative and quantitative errors introduced into the visibilities and power spectrum due to inpainting. We perform our analysis on simulated data as well as real data from the Hydrogen Epoch of Reionization Array (HERA) Phase 1 upper limits. We also introduce a convolutional neural network that capable of inpainting RFI corrupted data in interferometric instruments. We train our network on simulated data and show that our network is capable at inpainting real data without requiring to be retrained. We find that techniques that incorporate high wavenumbers in delay space in their modeling are best suited for inpainting over narrowband RFI. We also show that with our fiducial parameters Discrete Prolate Spheroidal Sequences (DPSS) and CLEAN provide the best performance for intermittent ``narrowband'' RFI while Gaussian Progress Regression (GPR) and Least Squares Spectral Analysis (LSSA) provide the best performance for larger RFI gaps. However we caution that these qualitative conclusions are sensitive to the chosen hyperparameters of each inpainting technique. We find these results to be consistent in both simulated and real visibilities. We show that all inpainting techniques reliably reproduce foreground dominated modes in the power spectrum. Since the inpainting techniques should not be capable of reproducing noise realizations, we find that the largest errors occur in the noise dominated delay modes. We show that in the future, as the noise level of the data comes down, CLEAN and DPSS are most capable of reproducing the fine frequency structure in the visibilities of HERA data.
△ Less
Submitted 20 February, 2023; v1 submitted 26 October, 2022;
originally announced October 2022.
-
The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28
Authors:
E. de Lera Acedo,
D. I. L. de Villiers,
N. Razavi-Ghods,
W. Handley,
A. Fialkov,
A. Magro,
D. Anstey,
H. T. J. Bevins,
R. Chiello,
J. Cumner,
A. T. Josaitis,
I. L. V. Roque,
P. H. Sims,
K. H. Scheutwinkel,
P. Alexander,
G. Bernardi,
S. Carey,
J. Cavillot,
W. Croukamp,
J. A. Ely,
T. Gessey-Jones,
Q. Gueuning,
R. Hills,
G. Kulkarni,
R. Maiolino
, et al. (9 additional authors not shown)
Abstract:
Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The fea…
▽ More
Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The feature is deeper than expected, and, if confirmed, would call for new physics. However, different groups have re-analyzed the EDGES data and questioned the reliability of the signal. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is a sky-averaged 21-cm experiment aiming at improving the current observations by tackling the issues faced by current instruments related to residual systematic signals in the data. The novel experimental approach focuses on detecting and jointly explaining these systematics together with the foregrounds and the cosmological signal using Bayesian statistics. To achieve this, REACH features simultaneous observations with two different antennas, an ultra wideband system (redshift range 7.5 to 28), and a receiver calibrator based on in-field measurements. Simulated observations forecast percent-level constraints on astrophysical parameters, potentially opening up a new window to the infant Universe.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
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.
-
Use of Time Dependent Data in Bayesian Global 21cm Foreground and Signal Modelling
Authors:
Dominic Anstey,
Eloy de Lera Acedo,
Will Handley
Abstract:
Global 21cm cosmology aims to investigate the cosmic dawn and epoch of reionisation by measuring the sky averaged HI absorption signal, which requires, accurate modelling of, or correction for, the bright radio foregrounds and distortions arising from chromaticity of the antenna beam. We investigate the effect of improving foreground modelling by fitting data sets from many observation times simul…
▽ More
Global 21cm cosmology aims to investigate the cosmic dawn and epoch of reionisation by measuring the sky averaged HI absorption signal, which requires, accurate modelling of, or correction for, the bright radio foregrounds and distortions arising from chromaticity of the antenna beam. We investigate the effect of improving foreground modelling by fitting data sets from many observation times simultaneously in a single Bayesian analysis, fitting for the same parameter set by performing these fits on simulated data. We find that for a hexagonal dipole antenna, this simultaneous fitting produces a significant improvement in the accuracy of the recovered 21cm signal, relative to fitting a time average of the data. Furthermore, the recovered models of the foreground are also seen to become more accurate by up to a factor of $\sim$2-3 relative to time averaged fitting. For a less chromatic log spiral antenna, no significant improvement in signal recovery was found by this process. However, the modelling of the foregrounds was still significantly improved. We also investigate extending this technique to fit multiple data sets from different antennae simultaneously for the same parameters. This is also found to improve both 21cm signal and foreground modelling, to a higher degree than fitting data set from multiple times from the same antenna.
△ Less
Submitted 16 January, 2023; v1 submitted 10 October, 2022;
originally announced October 2022.
-
Impact of instrument and data characteristics in the interferometric reconstruction of the 21 cm power spectrum
Authors:
Adélie Gorce,
Samskruthi Ganjam,
Adrian Liu,
Steven G. Murray,
Zara 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,
Carina Cheng,
David R. DeBoer,
Eloy de Lera Acedo,
Matt Dexter,
Joshua S. Dillon
, et al. (53 additional authors not shown)
Abstract:
Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand t…
▽ More
Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on these window functions. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that long-baseline observations correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage, whilst an informed choice of bandwidth and frequency taper can reduce said tails. With simple test cases and realistic simulations, we show that, apart from tracing mode mixing, the window functions help accurately reconstruct the power spectrum estimator of simulated visibilities. The window functions depend strongly on the beam chromaticity, and less on its spatial structure - a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented here correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and will be used in future analyses to confront theoretical models and data directly in cylindrical space.
△ Less
Submitted 11 January, 2023; v1 submitted 7 October, 2022;
originally announced October 2022.
-
Marginal Bayesian Statistics Using Masked Autoregressive Flows and Kernel Density Estimators with Examples in Cosmology
Authors:
Harry Bevins,
Will Handley,
Pablo Lemos,
Peter Sims,
Eloy de Lera Acedo,
Anastasia Fialkov
Abstract:
Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments.…
▽ More
Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments. However, these experiments are typically plagued by differing systematics, instrumental effects and contaminating signals, which we collectively refer to as `nuisance' components, that have to be modelled alongside target signals of interest. This leads to high dimensional parameter spaces, especially when combining data sets, with > 20 dimensions of which only around 5 correspond to key physical quantities. We present a means by which to combine constraints from different data sets in a computationally efficient manner by generating rapid, reusable and reliable marginal probability density estimators, giving us access to nuisance-free likelihoods. This is possible through the unique combination of nested sampling, which gives us access to Bayesian evidences, and the marginal Bayesian statistics code MARGARINE. Our method is lossless in the signal parameters, resulting in the same posterior distributions as would be found from a full nested sampling run over all nuisance parameters, and typically quicker than evaluating full likelihoods. We demonstrate our approach by applying it to the combination of posteriors from the Dark Energy Survey and Planck.
△ Less
Submitted 24 November, 2022; v1 submitted 23 July, 2022;
originally announced July 2022.
-
Marginal Post Processing of Bayesian Inference Products with Normalizing Flows and Kernel Density Estimators
Authors:
Harry T. J. Bevins,
William J. Handley,
Pablo Lemos,
Peter H. Sims,
Eloy de Lera Acedo,
Anastasia Fialkov,
Justin Alsing
Abstract:
Bayesian analysis has become an indispensable tool across many different cosmological fields including the study of gravitational waves, the Cosmic Microwave Background and the 21-cm signal from the Cosmic Dawn among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrument…
▽ More
Bayesian analysis has become an indispensable tool across many different cosmological fields including the study of gravitational waves, the Cosmic Microwave Background and the 21-cm signal from the Cosmic Dawn among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrumental effects modelled with `nuisance parameters'. In this paper, we summarise a method that uses Masked Autoregressive Flows and Kernel Density Estimators to learn marginal posterior densities corresponding to core science parameters. We find that the marginal or 'nuisance-free' posteriors and the associated likelihoods have an abundance of applications including; the calculation of previously intractable marginal Kullback-Leibler divergences and marginal Bayesian Model Dimensionalities, likelihood emulation and prior emulation. We demonstrate each application using toy examples, examples from the field of 21-cm cosmology and samples from the Dark Energy Survey. We discuss how marginal summary statistics like the Kullback-Leibler divergences and Bayesian Model Dimensionalities can be used to examine the constraining power of different experiments and how we can perform efficient joint analysis by taking advantage of marginal prior and likelihood emulators. We package our multipurpose code up in the pip-installable code margarine for use in the wider scientific community.
△ Less
Submitted 18 December, 2023; v1 submitted 25 May, 2022;
originally announced May 2022.
-
Bayesian Data Analysis for Sky-averaged 21-cm Experiments in the Presence of Ionospheric Effects
Authors:
Emma Shen,
Dominic Anstey,
Eloy de Lera Acedo,
Anastasia Fialkov
Abstract:
The ionosphere introduces chromatic distortions on low frequency radio waves, and thus poses a hurdle for 21-cm cosmology. In this paper we introduce time-varying chromatic ionospheric effects on simulated antenna temperature data of a global 21-cm data analysis pipeline, and try to detect the injected global signal. We demonstrate that given turbulent ionospheric conditions, more than 5\% error i…
▽ More
The ionosphere introduces chromatic distortions on low frequency radio waves, and thus poses a hurdle for 21-cm cosmology. In this paper we introduce time-varying chromatic ionospheric effects on simulated antenna temperature data of a global 21-cm data analysis pipeline, and try to detect the injected global signal. We demonstrate that given turbulent ionospheric conditions, more than 5\% error in our knowledge of the ionospheric parameters could lead to comparatively low evidence and high root-mean-square error (RMSE), suggesting a false or null detection. When using a constant antenna beam for cases that include data at different times, the significance of the detection lowers as the number of time samples increases. It is also shown that for observations that include data at different times, readjusting beam configurations according to the time-varying ionospheric conditions should greatly improve the significance of a detection, yielding higher evidences and lower RMSE, and that it is a necessary procedure for a successful detection when the ionospheric conditions are not ideal.
△ Less
Submitted 22 April, 2022;
originally announced April 2022.
-
Direct Optimal Mapping for 21cm Cosmology: A Demonstration with the Hydrogen Epoch of Reionization Array
Authors:
Zhilei Xu,
Jacqueline N. Hewitt,
Kai-Feng Chen,
Honggeun Kim,
Joshua S. Dillon,
Nicholas S. Kern,
Miguel F. Morales,
Bryna J. Hazelton,
Ruby Byrne,
Nicolas Fagnoni,
Eloy de Lera Acedo,
Zara 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
, et al. (56 additional authors not shown)
Abstract:
Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipe…
▽ More
Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays.
△ Less
Submitted 26 October, 2022; v1 submitted 12 April, 2022;
originally announced April 2022.
-
Bayesian evidence-driven likelihood selection for sky-averaged 21-cm signal extraction
Authors:
K. H. Scheutwinkel,
W. Handley,
E. de Lera Acedo
Abstract:
We demonstrate that the Bayesian evidence can be used to find a good approximation of the ground truth likelihood function of a dataset, a goal of the likelihood-free inference (LFI) paradigm. As a concrete example, we use forward modelled sky-averaged 21-cm signal antenna temperature datasets where we artificially inject noise structures of various physically motivated forms. We find that the Gau…
▽ More
We demonstrate that the Bayesian evidence can be used to find a good approximation of the ground truth likelihood function of a dataset, a goal of the likelihood-free inference (LFI) paradigm. As a concrete example, we use forward modelled sky-averaged 21-cm signal antenna temperature datasets where we artificially inject noise structures of various physically motivated forms. We find that the Gaussian likelihood performs poorly when the noise distribution deviates from the Gaussian case e.g. heteroscedastic radiometric or heavy-tailed noise. For these non-Gaussian noise structures, we show that the generalised normal likelihood is on a similar Bayesian evidence scale with comparable sky-averaged 21-cm signal recovery as the ground truth likelihood function of our injected noise. We therefore propose the generalised normal likelihood function as a good approximation of the true likelihood function if the noise structure is a priori unknown.
△ Less
Submitted 20 March, 2023; v1 submitted 9 April, 2022;
originally announced April 2022.
-
Bayesian evidence-driven diagnosis of instrumental systematics for sky-averaged 21-cm cosmology experiments
Authors:
K. H. Scheutwinkel,
E. de Lera Acedo,
W. Handley
Abstract:
We demonstrate the effectiveness of a Bayesian evidence-based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include…
▽ More
We demonstrate the effectiveness of a Bayesian evidence-based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include sky-averaged 21-cm posterior estimates resembling a very deep or wide signal. However, when including parameterised models of the systematic, the signal recovery is dramatically improved in performance. Most importantly, a Bayesian evidence-based model comparison is capable of determining whether or not such a systematic model is needed as the true underlying generative model of an experimental dataset is in principle unknown. We, therefore, advocate a pipeline capable of testing a variety of potential systematic errors with the Bayesian evidence acting as the mechanism for detecting their presence.
△ Less
Submitted 23 September, 2022; v1 submitted 9 April, 2022;
originally announced April 2022.
-
Impact of the Primordial Stellar Initial Mass Function on the 21-cm Signal
Authors:
T. Gessey-Jones,
N. S. Sartorio,
A. Fialkov,
G. M. Mirouh,
M. Magg,
R. G. Izzard,
E. de Lera Acedo,
W. J. Handley,
R. Barkana
Abstract:
Properties of the first generation of stars (Pop III), such as their initial mass function (IMF), are poorly constrained by observations and have yet to converge between simulations. The cosmological 21-cm signal of neutral hydrogen is predicted to be sensitive to Lyman-band photons produced by these stars, thus providing a unique way to probe the first stellar population. In this paper, we invest…
▽ More
Properties of the first generation of stars (Pop III), such as their initial mass function (IMF), are poorly constrained by observations and have yet to converge between simulations. The cosmological 21-cm signal of neutral hydrogen is predicted to be sensitive to Lyman-band photons produced by these stars, thus providing a unique way to probe the first stellar population. In this paper, we investigate the impacts of the Pop III IMF on the cosmic dawn 21-cm signal via the Wouthuysen-Field effect, Lyman-Werner feedback, Ly-alpha heating, and CMB heating. We calculate the emission spectra of star-forming halos for different IMFs by integrating over individual metal-free stellar spectra, computed from a set of stellar evolution histories and stellar atmospheres, and taking into account variability of the spectra with stellar age. Through this study, we therefore relax two common assumptions: that the zero-age main sequence emission rate of a Pop III star is representative of its lifetime mean emission rate, and that Pop III emission can be treated as instantaneous. Exploring a bottom-heavy, a top-heavy, and intermediate IMFs, we show that variations in the 21-cm signal are driven by stars lighter than 20 solar masses. For the explored models we find maximum relative differences of 59% in the cosmic dawn global 21-cm signal, and 131% between power spectra. Although this impact is modest, precise modelling of the first stars and their evolution is necessary for accurate prediction and interpretation of the 21-cm signal.
△ Less
Submitted 21 July, 2022; v1 submitted 4 February, 2022;
originally announced February 2022.
-
A Comprehensive Bayesian re-analysis of the SARAS2 data from the Epoch of Reionization
Authors:
H. T. J. Bevins,
E. de Lera Acedo,
A. Fialkov,
W. J. Handley,
S. Singh,
R. Subrahmanyan,
R. Barkana
Abstract:
We present a Bayesian re-analysis of the sky-averaged 21-cm experimental data from SARAS2 using nested sampling implemented with polychord, spectrally smooth foreground modelling implemented with maxsmooth, detailed systematic modelling and rapid signal emulation with globalemu. Our analysis differs from previous analysis of the SARAS2 data through the use of a full Bayesian framework and separate…
▽ More
We present a Bayesian re-analysis of the sky-averaged 21-cm experimental data from SARAS2 using nested sampling implemented with polychord, spectrally smooth foreground modelling implemented with maxsmooth, detailed systematic modelling and rapid signal emulation with globalemu. Our analysis differs from previous analysis of the SARAS2 data through the use of a full Bayesian framework and separate modelling of the foreground and non-smooth systematics. We use the most up-to-date signal models including Lyman-$α$ and CMB heating parameterised by astrophysical parameters such as star formation efficiency, X-ray heating efficiency, minimal virial circular velocity of star forming galaxies, CMB optical depth and the low energy cutoff of the X-ray spectral energy distribution. We consider models with an excess radio background above the CMB produced via radio emission from early galaxies and parameterised by a radio production efficiency. A non-smooth systematic is identified and modelled as both a frequency damped sinusoid introduced by the electronics and separately from the sky. The latter is modulated by the total efficiency of the antenna and marginally favoured by the data. We consider three different models for the noise in the data. The SARAS2 constraints on individual astrophysical parameters are extremely weak however we identify classes of disfavoured signals. We weakly disfavour standard astrophysical models with high Lyman-$α$ fluxes and weak heating and more confidently disfavour exotic models with high Lyman-$α$ fluxes, low X-ray efficiencies and high radio production efficiencies in early galaxies.
△ Less
Submitted 20 May, 2022; v1 submitted 27 January, 2022;
originally announced January 2022.
-
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.
-
Array Element Coupling in Radio Interferometry I: A Semi-Analytic Approach
Authors:
Alec T. Josaitis,
Aaron Ewall-Wice,
Nicolas Fagnoni,
Eloy de Lera Acedo
Abstract:
We derive a general formalism for interferometric visibilities, which considers first-order antenna-antenna coupling and assumes steady-state, incident radiation. We simulate such coupling features for non-polarized skies on a compact, redundantly-spaced array and present a phenomenological analysis of the coupling features. Contrary to previous studies, we find mutual coupling features manifest t…
▽ More
We derive a general formalism for interferometric visibilities, which considers first-order antenna-antenna coupling and assumes steady-state, incident radiation. We simulate such coupling features for non-polarized skies on a compact, redundantly-spaced array and present a phenomenological analysis of the coupling features. Contrary to previous studies, we find mutual coupling features manifest themselves at nonzero fringe rates. We compare power spectrum results for both coupled and non-coupled (noiseless, simulated) data and find coupling effects to be highly dependent on LST, baseline length, and baseline orientation. For all LSTs, lengths, and orientations, coupling features appear at delays which are outside the foreground 'wedge', which has been studied extensively and contains non-coupled astrophysical foreground features. Further, we find that first-order coupling effects threaten our ability to average data from baselines with identical length and orientation. Two filtering strategies are proposed which may mitigate such coupling systematics. The semi-analytic coupling model herein presented may be used to study mutual coupling systematics as a function of LST, baseline length, and baseline orientation. Such a model is not only helpful to the field of 21cm cosmology, but any study involving interferometric measurements, where coupling effects at the level of at least 1 part in 10^4 could corrupt the scientific result. Our model may be used to mitigate coupling systematics in existing radio interferometers and to design future arrays where the configuration of array elements inherently mitigates coupling effects at desired LSTs and angular resolutions.
△ Less
Submitted 4 April, 2022; v1 submitted 20 October, 2021;
originally announced October 2021.
-
The Aperture Array Verification System 1: System overview and early commissioning results
Authors:
P. Benthem,
R. Wayth,
E. de Lera Acedo,
K. Zarb Adami,
M. Alderighi,
C. Belli,
P. Bolli,
T. Booler,
J. Borg,
J. W. Broderick,
S. Chiarucci,
R. Chiello,
L. Ciani,
G. Comoretto,
B. Crosse,
D. Davidson,
A. DeMarco,
D. Emrich,
A. van Es,
D. Fierro,
A. Faulkner,
M. Gerbers,
N. Razavi-Ghods,
P. Hall,
L. Horsley
, et al. (29 additional authors not shown)
Abstract:
The design and development process for the Square Kilometre Array (SKA) radio telescope, the Low Frequency Aperture Array component, was progressed during the SKA pre-construction phase by an international consortium, with the goal of meeting requirements for a critical design review. As part of the development process a full-sized prototype SKA Low station was deployed, the Aperture Array Verific…
▽ More
The design and development process for the Square Kilometre Array (SKA) radio telescope, the Low Frequency Aperture Array component, was progressed during the SKA pre-construction phase by an international consortium, with the goal of meeting requirements for a critical design review. As part of the development process a full-sized prototype SKA Low station was deployed, the Aperture Array Verification System 1 (AAVS1). We provide a system overview and describe the commissioning results of AAVS1, which is a low frequency radio telescope with 256 dual-polarisation log-periodic dipole antennas working as a phased array. A detailed system description is provided, including an in-depth overview of relevant sub-systems, ranging from hardware, firmware, software, calibration,and control sub-systems. Early commissioning results cover initial bootstrapping, array calibration, stability testing, beam-forming,and on-sky sensitivity validation. Lessons learned are presented, along with future developments.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
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.
-
Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case
Authors:
J. Cumner,
E. De Lera Acedo,
D. I. L. de Villiers,
D. Anstey,
C. I. Kolitsidas,
B. Gurdon,
N. Fagnoni,
P. Alexander,
G. Bernardi,
H. T. J. Bevins,
S. Carey,
J. Cavillot,
R. Chiello,
C. Craeye,
W. Croukamp,
J. A. Ely,
A. Fialkov,
T. Gessey-Jones,
Q. Gueuning,
W. Handley,
R. Hills,
A. T. Josaitis,
G. Kulkarni,
A. Magro,
R. Maiolino
, et al. (13 additional authors not shown)
Abstract:
Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experim…
▽ More
Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis.
△ Less
Submitted 12 January, 2023; v1 submitted 21 September, 2021;
originally announced September 2021.
-
HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization
Authors:
The HERA Collaboration,
Zara Abdurashidova,
James E. Aguirre,
Paul Alexander,
Zaki Ali,
Yanga Balfour,
Rennan Barkana,
Adam Beardsley,
Gianni Bernardi,
Tashalee Billings,
Judd Bowman,
Richard Bradley,
Phillip Bull,
Jacob Burba,
Steven Carey,
Christopher Carilli,
Carina Cheng,
David DeBoer,
Matthew Dexter,
Eloy de Lera Acedo,
Joshua Dillon,
John Ely,
Aaron Ewall-Wice,
Nicolas Fagnoni,
Anastasia Fialkov
, et al. (59 additional authors not shown)
Abstract:
Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated…
▽ More
Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,ν}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results.
△ Less
Submitted 20 December, 2022; v1 submitted 16 August, 2021;
originally announced August 2021.
-
First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum
Authors:
The HERA Collaboration,
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,
Steve Carey,
Chris L. Carilli,
Carina Cheng,
David R. DeBoer,
Matt Dexter,
Eloy de Lera Acedo,
Taylor Dibblee-Barkman,
Joshua S. Dillon,
John Ely,
Aaron Ewall-Wice,
Nicolas Fagnoni,
Randall Fritz
, et al. (52 additional authors not shown)
Abstract:
We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground…
▽ More
We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.
△ Less
Submitted 4 August, 2021;
originally announced August 2021.
-
Informing antenna design for sky-averaged 21-cm experiments using a simulated Bayesian data analysis pipeline
Authors:
Dominic Anstey,
John Cumner,
Eloy de Lera Acedo,
Will Handley
Abstract:
Global 21cm experiments aim to measure the sky averaged HI absorption signal from cosmic dawn and the epoch of reionisation. However, antenna chromaticity coupling to bright foregrounds can introduce distortions into the observational data of such experiments. We demonstrate a method for guiding the antenna design of a global experiment through data analysis simulations. This is done by performing…
▽ More
Global 21cm experiments aim to measure the sky averaged HI absorption signal from cosmic dawn and the epoch of reionisation. However, antenna chromaticity coupling to bright foregrounds can introduce distortions into the observational data of such experiments. We demonstrate a method for guiding the antenna design of a global experiment through data analysis simulations. This is done by performing simulated observations for a range of inserted 21cm signals, then attempting to identify the signals with a data analysis pipeline. We demonstrate this method on five antennae that were considered as potential designs for the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH); a conical log spiral antenna, an inverted conical sinuous antenna and polygonal-, square- and elliptical-bladed dipoles. We find that the log spiral performs significantly better than the other antennae tested, able to correctly and confidently identify every inserted 21cm signal. In second place is the polygonal dipole antenna, which was only unable to detect signals with both very low amplitudes of 0.05K and low centre frequency of 80MHz. The conical sinuous antenna was found to perform least accurately, only able to detect the highest amplitude 21cm signals, and even then with biases. We also demonstrate that, due to the non-trivial nature of chromatic distortion and the processes of correcting for it, these are not the results that could have been expected superficially from the extent of chromatic variation in each antenna.
△ Less
Submitted 11 November, 2021; v1 submitted 18 June, 2021;
originally announced June 2021.