-
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.
-
All Sky Modelling Requirements for Bayesian 21 cm Power Spectrum Estimation with BayesEoR
Authors:
Jacob Burba,
Peter H. Sims,
Jonathan C. Pober
Abstract:
We present a comprehensive simulation-based study of the BayesEoR code for 21 cm power spectrum recovery when analytically marginalizing over foreground parameters. To account for covariance between the 21 cm signal and contaminating foreground emission, BayesEoR jointly constructs models for both signals within a Bayesian framework. Due to computational constraints, the forward model is construct…
▽ More
We present a comprehensive simulation-based study of the BayesEoR code for 21 cm power spectrum recovery when analytically marginalizing over foreground parameters. To account for covariance between the 21 cm signal and contaminating foreground emission, BayesEoR jointly constructs models for both signals within a Bayesian framework. Due to computational constraints, the forward model is constructed using a restricted field-of-view (FoV) in the image domain. When the only EoR contaminants are noise and foregrounds, we demonstrate that BayesEoR can accurately recover the 21 cm power spectrum when the component of sky emission outside this forward-modelled region is downweighted by the beam at the level of the dynamic range between the foreground and 21 cm signals. However, when all-sky foreground emission is included along with a realistic instrument primary beam with sidelobes above this threshold extending to the horizon, the recovered power spectrum is contaminated by unmodelled sky emission outside the restricted FoV model. Expanding the combined cosmological and foreground model to cover the whole sky is computationally prohibitive. To address this, we present a modified version of BayesEoR that allows for an all-sky foreground model, while the modelled 21 cm signal remains only within the primary FoV of the telescope. With this modification, it will be feasible to run an all-sky BayesEoR analysis on a sizeable compute cluster. We also discuss several future directions for further reducing the need to model all-sky foregrounds, including wide-field foreground subtraction, an image-domain likelihood utilizing a tapering function, and instrument primary beam design.
△ Less
Submitted 8 February, 2023;
originally announced February 2023.
-
Analytic approximations of scattering effects on beam chromaticity in 21-cm global experiments
Authors:
Alan E. E. Rogers,
John P. Barrett,
Judd D. Bowman,
Rigel Cappallo,
Colin J. Lonsdale,
Nivedita Mahesh,
Raul A. Monsalve,
Steven G. Murray,
Peter H. Sims
Abstract:
Scattering from objects near an antenna produce correlated signals from strong compact radio sources in a manner similar to those used by the Sea Interferometer to measure the radio source positions using the fine frequency structure in the total power spectrum of a single antenna. These fringes or ripples due to correlated signal interference are present at a low level in the spectrum of any sing…
▽ More
Scattering from objects near an antenna produce correlated signals from strong compact radio sources in a manner similar to those used by the Sea Interferometer to measure the radio source positions using the fine frequency structure in the total power spectrum of a single antenna. These fringes or ripples due to correlated signal interference are present at a low level in the spectrum of any single antenna and are a major source of systematics in systems used to measure the global redshifted 21-cm signal from the early universe. In the Sea Interferometer a single antenna on a cliff above the sea is used to add the signal from the direct path to the signal from the path reflected from the sea thereby forming an interferometer. This was used for mapping radio sources with a single antenna by Bolton and Slee in the 1950s. In this paper we derive analytic expressions to determine the level of these ripples and compare these results in a few simple cases with electromagnetic modeling software to verify that the analytic calculations are sufficient to obtain the magnitude of the scattering effects on the measurements of the global 21-cm signal. These analytic calculations are needed to evaluate the magnitude of the effects in cases that are either too complex or take too much time to be modeled using software.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
A Bayesian approach to modelling spectrometer data chromaticity corrected using beam factors -- I. Mathematical formalism
Authors:
Peter H. Sims,
Judd D. Bowman,
Nivedita Mahesh,
Steven G. Murray,
John P. Barrett,
Rigel Cappallo,
Raul A. Monsalve,
Alan E. E. Rogers,
Titu Samson,
Akshatha K. Vydula
Abstract:
Accurately accounting for spectral structure in spectrometer data induced by instrumental chromaticity on scales relevant for detection of the 21-cm signal is among the most significant challenges in global 21-cm signal analysis. In the publicly available EDGES low-band data set, this complicating structure is suppressed using beam-factor based chromaticity correction (BFCC), which works by dividi…
▽ More
Accurately accounting for spectral structure in spectrometer data induced by instrumental chromaticity on scales relevant for detection of the 21-cm signal is among the most significant challenges in global 21-cm signal analysis. In the publicly available EDGES low-band data set, this complicating structure is suppressed using beam-factor based chromaticity correction (BFCC), which works by dividing the data by a sky-map-weighted model of the spectral structure of the instrument beam. Several analyses of this data have employed models that start with the assumption that this correction is complete. However, while BFCC mitigates the impact of instrumental chromaticity on the data, given realistic assumptions regarding the spectral structure of the foregrounds, the correction is only partial. This complicates the interpretation of fits to the data with intrinsic sky models (models that assume no instrumental contribution to the spectral structure of the data). In this paper, we derive a BFCC data model from an analytic treatment of BFCC and demonstrate using simulated observations that, in contrast to using an intrinsic sky model for the data, the BFCC data model enables unbiased recovery of a simulated global 21-cm signal from beam-factor chromaticity corrected data in the limit that the data is corrected with an error-free beam-factor model.
△ Less
Submitted 28 March, 2023; v1 submitted 7 December, 2022;
originally announced December 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.
-
A Bayesian Calibration Framework for EDGES
Authors:
Steven G. Murray,
Judd D. Bowman,
Peter H. Sims,
Nivedita Mahesh,
Alan E. E. Rogers,
Raul A. Monsalve,
Titu Samson,
Akshatha Konakondula Vydula
Abstract:
We develop a Bayesian model that jointly constrains receiver calibration, foregrounds and cosmic 21cm signal for the EDGES global 21\,cm experiment. This model simultaneously describes calibration data taken in the lab along with sky-data taken with the EDGES low-band antenna. We apply our model to the same data (both sky and calibration) used to report evidence for the first star formation in 201…
▽ More
We develop a Bayesian model that jointly constrains receiver calibration, foregrounds and cosmic 21cm signal for the EDGES global 21\,cm experiment. This model simultaneously describes calibration data taken in the lab along with sky-data taken with the EDGES low-band antenna. We apply our model to the same data (both sky and calibration) used to report evidence for the first star formation in 2018. We find that receiver calibration does not contribute a significant uncertainty to the inferred cosmic signal (<1%), though our joint model is able to more robustly estimate the cosmic signal for foreground models that are otherwise too inflexible to describe the sky data. We identify the presence of a significant systematic in the calibration data, which is largely avoided in our analysis, but must be examined more closely in future work. Our likelihood provides a foundation for future analyses in which other instrumental systematics, such as beam corrections and reflection parameters, may be added in a modular manner.
△ Less
Submitted 7 September, 2022;
originally announced September 2022.
-
A Bayesian approach to high fidelity interferometric calibration I: mathematical formalism
Authors:
Peter H. Sims,
Jonathan C. Pober,
Jonathan L. Sievers
Abstract:
High fidelity radio interferometric data calibration that minimises spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been show…
▽ More
High fidelity radio interferometric data calibration that minimises spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been shown to impart spurious spectral structure to the calibrated data if the sky model used to calibrate the data is incomplete. In this paper, we introduce BayesCal: a novel solution to the sky-model incompleteness problem in interferometric calibration, designed to enable high fidelity data calibration. The BayesCal data model supplements the a priori known component of the forward model of the sky with a statistical model for the missing and uncertain flux contribution to the data, constrained by a prior on the power in the model. We demonstrate how the parameters of this model can be marginalised out analytically, reducing the dimensionality of the parameter space to be sampled from and allowing one to sample directly from the posterior probability distribution of the calibration parameters. Additionally, we show how physically motivated priors derived from theoretical and measurement-based constraints on the spectral smoothness of the instrumental gains can be used to constrain the calibration solutions. In a companion paper, we apply this algorithm to simulated observations with a HERA-like array and demonstrate that it enables up to four orders of magnitude suppression of power in spurious spectral fluctuations relative to standard calibration approaches.
△ Less
Submitted 27 June, 2022;
originally announced June 2022.
-
A Bayesian approach to high fidelity interferometric calibration II: demonstration with simulated data
Authors:
Peter H. Sims,
Jonathan C. Pober,
Jonathan L. Sievers
Abstract:
In a companion paper, we presented BayesCal, a mathematical formalism for mitigating sky-model incompleteness in interferometric calibration. In this paper, we demonstrate the use of BayesCal to calibrate the degenerate gain parameters of full-Stokes simulated observations with a HERA-like hexagonal close-packed redundant array, for three assumed levels of completeness of the a priori known compon…
▽ More
In a companion paper, we presented BayesCal, a mathematical formalism for mitigating sky-model incompleteness in interferometric calibration. In this paper, we demonstrate the use of BayesCal to calibrate the degenerate gain parameters of full-Stokes simulated observations with a HERA-like hexagonal close-packed redundant array, for three assumed levels of completeness of the a priori known component of the calibration sky model. We compare the BayesCal calibration solutions to those recovered by calibrating the degenerate gain parameters with only the a priori known component of the calibration sky model both with and without imposing physically motivated priors on the gain amplitude solutions and for two choices of baseline length range over which to calibrate. We find that BayesCal provides calibration solutions with up to four orders of magnitude lower power in spurious gain amplitude fluctuations than the calibration solutions derived for the same data set with the alternate approaches, and between $\sim10^7$ and $\sim10^{10}$ times smaller than in the mean degenerate gain amplitude on the full range of spectral scales accessible in the data. Additionally, we find that in the scenarios modelled only BayesCal has sufficiently high fidelity calibration solutions for unbiased recovery of the 21 cm power spectrum on large spectral scales ($k_\parallel \lesssim 0.15~h\mathrm{Mpc}^{-1}$). In all other cases, in the completeness regimes studied, those scales are contaminated.
△ Less
Submitted 27 June, 2022;
originally announced June 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.
-
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.
-
Testing for calibration systematics in the EDGES low-band data using Bayesian model selection
Authors:
Peter H. Sims,
Jonathan C. Pober
Abstract:
Cosmic Dawn, when the first stars and proto-galaxies began to form, is commonly expected to be accompanied by an absorption signature at radio frequencies. This feature arises as Lyman-$α$ photons emitted by these first luminous objects couple the 21 cm excitation temperature of intergalactic hydrogen gas to its kinetic temperature, driving it into absorption relative to the CMB. The detailed prop…
▽ More
Cosmic Dawn, when the first stars and proto-galaxies began to form, is commonly expected to be accompanied by an absorption signature at radio frequencies. This feature arises as Lyman-$α$ photons emitted by these first luminous objects couple the 21 cm excitation temperature of intergalactic hydrogen gas to its kinetic temperature, driving it into absorption relative to the CMB. The detailed properties of this absorption profile encode powerful information about the physics of Cosmic Dawn. Recently, Bowman et al. analysed data from the EDGES low-band radio antenna and found an unexpectedly deep absorption profile centred at 78 MHz, which could be a detection of this signature. Their specific analysis fit their measurements using a polynomial foreground model, a flattened Gaussian absorption profile and a white noise model; we argue that a more accurate model, that includes a detailed noise model and accounting for the effects of plausible calibration errors, is essential for describing the EDGES data set. We perform a Bayesian evidence-based comparison of models of the EDGES low-band data set and find that those incorporating these additional components are decisively preferred. The subset of the best fitting models of the data that include a global signal favour an amplitude consistent with standard cosmological assumptions (A < 209 mK). However, there is not strong evidence to favour models of the data including a global 21 cm signal over those without one. Ultimately, we find that the derivation of robust constraints on astrophysics from the data is limited by the presence of systematics.
△ Less
Submitted 7 October, 2019;
originally announced October 2019.
-
A Roadmap for Astrophysics and Cosmology with High-Redshift 21 cm Intensity Mapping
Authors:
The Hydrogen Epoch of Reionization Array,
Collaboration,
James E. Aguirre,
Adam P. Beardsley,
Gianni Bernardi,
Judd D. Bowman,
Philip Bull,
Chris L. Carilli,
Wei-Ming Dai,
David R. DeBoer,
Joshua S. Dillon,
Aaron Ewall-Wice,
Steve R. Furlanetto,
Bharat K. Gehlot,
Deepthi Gorthi,
Bradley Greig,
Bryna J. Hazelton,
Jacqueline N. Hewitt,
Daniel C. Jacobs,
Nicholas S. Kern,
Piyanat Kittiwisit,
Matthew Kolopanis,
Paul La Plante,
Adrian Liu,
Yin-Zhe Ma
, et al. (7 additional authors not shown)
Abstract:
In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the syste…
▽ More
In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the systematics confronting 21 cm cosmology into a plan for overcoming them, enabling next-generation, mid-scale 21 cm arrays to be proposed late in the decade. Submitted for consideration by the Astro2020 Decadal Survey Program Panel for Radio, Millimeter, and Submillimeter Observations from the Ground as a Medium-Sized Project.
△ Less
Submitted 15 July, 2019;
originally announced July 2019.
-
Joint estimation of the Epoch of Reionization power spectrum and foregrounds
Authors:
Peter H. Sims,
Jonathan C. Pober
Abstract:
The power spectrum of redshifted 21 cm emission brightness temperature fluctuations is a powerful probe of the Epoch of Reionization (EoR). However, bright foreground emission presents a significant impediment to its unbiased recovery from interferometric data. We build on the Bayesian power spectral estimation methodology introduced in Sims et al. 2016 and demonstrate that incorporating a priori…
▽ More
The power spectrum of redshifted 21 cm emission brightness temperature fluctuations is a powerful probe of the Epoch of Reionization (EoR). However, bright foreground emission presents a significant impediment to its unbiased recovery from interferometric data. We build on the Bayesian power spectral estimation methodology introduced in Sims et al. 2016 and demonstrate that incorporating a priori knowledge of the spectral structure of foregrounds in the large spectral scale component of the data model enables significantly improved modelling of the foregrounds without increasing the model complexity. We explore two astrophysically motivated parametrisations of the large spectral scale model: (i) a constant plus power law model of the form $q_{0}+q_{1}(ν/ν_{0})^{b_{1}}$ for two values of $b_{1}$: $b_{1} = <β>_\mathrm{GDSE}$ and $b_{1} = <β>_\mathrm{EGS}$, the mean spectral indices of the Galactic diffuse synchrotron emission and extragalactic source foreground emission, respectively, and (ii) a constant plus double power law model of the form $q_{0}+q_{1}(ν/ν_{0})^{b_{1}}+q_{2}(ν/ν_{0})^{b_{2}}$ with $b_{1} = <β>_\mathrm{GDSE}$ and $b_{2} = <β>_\mathrm{EGS}$. We estimate the EoR power spectrum from simulated interferometric data consisting of an EoR signal, Galactic diffuse synchrotron emission, extragalactic sources and diffuse free-free emission from the Galaxy. We show that, by jointly estimating a model of the EoR signal with the constant plus double power law parametrisation of the large spectral scale model, unbiased estimates of the EoR power spectrum are recoverable on all spatial scales accessible in the data set, including on the large spatial scales that were found to be contaminated in earlier work.
△ Less
Submitted 4 July, 2019;
originally announced July 2019.
-
Bayesian power spectrum estimation at the Epoch of Reionization
Authors:
Peter H. Sims,
Lindley Lentati,
Jonathan C. Pober,
Chris Carilli,
Michael P. Hobson,
Paul Alexander,
Paul Sutter
Abstract:
We introduce a new method for performing robust Bayesian estimation of the three-dimensional spatial power spectrum at the Epoch of Reionization (EoR), from interferometric observations. The versatility of this technique allows us to present two approaches. First, when the observations span only a small number of independent spatial frequencies ($k$-modes) we sample directly from the spherical pow…
▽ More
We introduce a new method for performing robust Bayesian estimation of the three-dimensional spatial power spectrum at the Epoch of Reionization (EoR), from interferometric observations. The versatility of this technique allows us to present two approaches. First, when the observations span only a small number of independent spatial frequencies ($k$-modes) we sample directly from the spherical power spectrum coefficients that describe the EoR signal realisation. Second, when the number of $k$-modes to be included in the model becomes large, we sample from the joint probability density of the spherical power spectrum and the signal coefficients, using Hamiltonian Monte Carlo methods to explore this high dimensional ($\sim$ 20000) space efficiently. This approach has been successfully applied to simulated observations that include astrophysically realistic foregrounds in a companion publication (Sims et al. 2016). Here we focus on explaining the methodology in detail, and use simple foreground models to both demonstrate its efficacy, and highlight salient features. In particular, we show that including an arbitrary flat spectrum continuum foreground that is $10^8$ times greater in power than the EoR signal has no detectable impact on our parameter estimates of the EoR power spectrum recovered from the data.
△ Less
Submitted 3 November, 2018; v1 submitted 12 January, 2017;
originally announced January 2017.
-
Contamination of the Epoch of Reionization power spectrum in the presence of foregrounds
Authors:
Peter H. Sims,
Lindley Lentati,
Paul Alexander,
Chris L. Carilli
Abstract:
We construct foreground simulations comprising spatially correlated extragalactic and diffuse Galactic emission components and calculate the `intrinsic' (instrument-free) two-dimensional spatial power spectrum and the cylindrically and spherically averaged three-dimensional k-space power spectra of the Epoch of Reionization (EoR) and our foreground simulations using a Bayesian power spectral estim…
▽ More
We construct foreground simulations comprising spatially correlated extragalactic and diffuse Galactic emission components and calculate the `intrinsic' (instrument-free) two-dimensional spatial power spectrum and the cylindrically and spherically averaged three-dimensional k-space power spectra of the Epoch of Reionization (EoR) and our foreground simulations using a Bayesian power spectral estimation framework. This leads us to identify a model dependent region of optimal signal estimation for our foreground and EoR models, within which the spatial power in the EoR signal relative to foregrounds is maximised. We identify a target field dependent region, in k-space, of intrinsic foreground power spectral contamination at low k_perp and k_parallel and a transition to a relatively foreground-free intrinsic EoR window in the complement to this region. The contaminated region of k-space demonstrates that simultaneous estimation of the EoR and foregrounds is important for obtaining statistically robust estimates of the EoR power spectrum; biased results will be obtained from methodologies that ignore their covariance. Using simulated observations with frequency dependent uv-coverage and primary beam, with the former derived for HERA in 37-antenna and 331-antenna configuration, we recover instrumental power spectra consistent with their intrinsic counterparts. We discuss the implications of these results for optimal strategies for unbiased estimation of the EoR power spectrum.
△ Less
Submitted 26 July, 2016;
originally announced July 2016.