Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBayesian power spectrum estimation at the Epoch of Reionization
345
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic 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 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.
rate research
Read More
Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $Delta^2 leq 3.9 times 10^3$ mK$^2$ at $k=0.20$ $textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude.
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}( u/ u_{0})^{b_{1}}$ for two values of $b_{1}$: $b_{1} = <beta>_mathrm{GDSE}$ and $b_{1} = <beta>_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}( u/ u_{0})^{b_{1}}+q_{2}( u/ u_{0})^{b_{2}}$ with $b_{1} = <beta>_mathrm{GDSE}$ and $b_{2} = <beta>_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.
21 cm Epoch of Reionization observations promise to transform our understanding of galaxy formation, but these observations are impossible without unprecedented levels of instrument calibration. We present end-to-end simulations of a full EoR power spectrum analysis including all of the major components of a real data processing pipeline: models of astrophysical foregrounds and EoR signal, frequency-dependent instrument effects, sky-based antenna calibration, and the full PS analysis. This study reveals that traditional sky-based per-frequency antenna calibration can only be implemented in EoR measurement analyses if the calibration model is unrealistically accurate. For reasonable levels of catalog completeness, the calibration introduces contamination in otherwise foreground-free power spectrum modes, precluding a PS measurement. We explore the origin of this contamination and potential mitigation techniques. We show that there is a strong joint constraint on the precision of the calibration catalog and the inherent spectral smoothness of antennae, and that this has significant implications for the instrumental design of the SKA and other future EoR observatories.
Epoch of Reionization data analysis requires unprecedented levels of accuracy in radio interferometer pipelines. We have developed an imaging power spectrum analysis to meet these requirements and generate robust 21 cm EoR measurements. In this work, we build a signal path framework to mathematically describe each step in the analysis, from data reduction in the FHD package to power spectrum generation in the $varepsilon$ppsilon package. In particular, we focus on the distinguishing characteristics of FHD/$varepsilon$ppsilon: highly accurate spectral calibration, extensive data verification products, and end-to-end error propagation. We present our key data analysis products in detail to facilitate understanding of the prominent systematics in image-based power spectrum analyses. As a verification to our analysis, we also highlight a full-pipeline analysis simulation to demonstrate signal preservation and lack of signal loss. This careful treatment ensures that the FHD/$varepsilon$ppsilon power spectrum pipeline can reduce radio interferometric data to produce credible 21 cm EoR measurements.
We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERAs redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
