-
Verifying the Australian MWA EoR pipeline I: 21-cm sky model and correlated measurement density
Authors:
J. L. B. Line,
C. M. Trott,
J. H. Cook,
B. Greig,
N. Barry,
C. H. Jordan
Abstract:
We present the first of two papers dedicated to verifying the Australian Epoch of Reionisation pipeline (AusEoRPipe) through simulation. The AusEoRPipe aims to disentangle 21-cm radiation emitted by gas surrounding the very first stars from contaminating foreground astrophysical sources, and has been in development for close to a decade. In this paper, we build an accurate 21-cm sky model that can…
▽ More
We present the first of two papers dedicated to verifying the Australian Epoch of Reionisation pipeline (AusEoRPipe) through simulation. The AusEoRPipe aims to disentangle 21-cm radiation emitted by gas surrounding the very first stars from contaminating foreground astrophysical sources, and has been in development for close to a decade. In this paper, we build an accurate 21-cm sky model that can be used by the WODEN simulation software to create visibilities containing a predictable 21-cm signal. We verify that the power spectrum estimator CHIPS can recover this signal in the absence of foregrounds. We also investigate how measurements in Fourier-space are correlated, and how their gridded density affects the power spectrum. We measure and fit for this effect using Gaussian-noise simulations of the MWA phase I layout. We find a gridding density correction factor of 2.651 appropriate for integrations equal to or greater than 30 minutes of data, which contain observations with multiple primary beam pointings and LSTs. Paper II of this series will use the results of this paper to test the AusEoRPipe in the presence of foregrounds and instrumental effects.
△ Less
Submitted 7 April, 2024;
originally announced April 2024.
-
The Impact of ionization Morphology and X-ray Heating on the Cosmological 21cm Skew Spectrum
Authors:
J. H. Cook,
S. Balu,
B. Greig,
C. M. Trott,
J. L. B. Line,
Y. Qin,
J. S. B. Wyithe
Abstract:
The cosmological 21cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially-dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal however is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heatin…
▽ More
The cosmological 21cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially-dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal however is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heating and reionization, by calculating the skew spectrum (SS) of the 21cm signal using MERAXES, which couples a semi-analytic galaxy population with semi-numerical reionization simulations. The SS is the cross-spectrum of the quadratic temperature brightness field with itself. We generate a set of seven simulations from $z = 30$ to $z = 5$, varying the halo mass threshold for hosting star-formation, the X-ray luminosity per star-formation rate, and the minimum X-ray energy escaping host galaxies. We find the SS is predominantly negative as a function of redshift, transitioning to positive towards the start of reionization, and peaking during the midpoint of reionization. We do not see a negative dip in the SS during reionization, likely due to the specifics of modelling ionization sources. We normalise the SS by the PS during reionization isolating the non-Gaussianities. We find a trough ($k\sim\,0.1\,\textrm{Mpc}^{-1}$) and peak ($k\sim\,0.4-1\,\textrm{Mpc}^{-1}$) in the normalised SS during the mid to late periods of reionization. These correlate to the ionization topology, and neutral islands in the IGM. We calculate the cosmic variance of the normalised SS, and find these features are detectable in the absence of foregrounds with the SKA_LOW.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Investigating the Contribution of Extended Radio Sources to the Epoch of Reionisation Power Spectrum
Authors:
J H. Cook,
C. M. Trott,
J. L. B. Line
Abstract:
We investigate the contribution of extended radio sources such as Centaurus A, and Galactic supernova remnants (SNRs) to our ability to detect the statistical $21\,\rm{cm}$ signal from the Epoch of Reionisation (EoR) with the Murchison Widefield Array (MWA). These sources are typically ignored because they are in highly attenuated parts of the MWA primary beam, however in aggregate these sources h…
▽ More
We investigate the contribution of extended radio sources such as Centaurus A, and Galactic supernova remnants (SNRs) to our ability to detect the statistical $21\,\rm{cm}$ signal from the Epoch of Reionisation (EoR) with the Murchison Widefield Array (MWA). These sources are typically ignored because they are in highly attenuated parts of the MWA primary beam, however in aggregate these sources have apparent flux densities of $10\,\rm{Jy}$ on angular scales we expect to detect the $21\,\rm{cm}$ signal. We create bespoke multi-component 2D Gaussian models for Galactic SNRs and for Centaurus A, and simulate the visibilities for two MWA snapshot observations. We grid those visibilities and then Fourier transform them with respect to frequency, averaging them both spherically and cylindrically to produce the 1D and 2D power spectra. We compare the simulated 1D power spectra to the expected $21\,\rm{cm}$ power spectrum. We find that although these extended sources are in highly attenuated parts of the MWA primary beam pattern, collectively they have enough power ($\sim10^4-10^5\,\rm{mK^2}\,\it{h^{-3}} \,\rm{Mpc^{3}}$) on EoR significant modes ($|\mathbf{k}| \leq 0.1 h \rm{Mpc}^{-1}$) to prohibit detection of the $21\,\rm{cm}$ signal ($10^4\,\rm{mK^2}\,\it{h^{-3}} \,\rm{Mpc^{3}}$). We find that $50-90\%$ of sources must be removed in order to reduce leakage to a level of $10-20\%$ of the $21\,\rm{cm}$ power spectrum on EoR significant modes. The effects of widefield extended sources will have implications on the detectability of the $21\,\rm{cm}$ signal for the MWA and with the future Square Kilometre Array (SKA).
△ Less
Submitted 9 May, 2022;
originally announced May 2022.
-
Imaging at 300 MHz With the MWA
Authors:
J. H. Cook,
N. Seymour,
M. Sokolowski
Abstract:
At relatively high frequencies, highly sensitive grating sidelobes occur in the primary beam patterns of low frequency aperture arrays (LFAA) such as the Murchison Widefield Array (MWA). This occurs when the observing wavelength becomes comparable to the dipole separation for LFAA tiles, which for the MWA occurs at approximately 300 MHz. The presence of these grating sidelobes has made calibration…
▽ More
At relatively high frequencies, highly sensitive grating sidelobes occur in the primary beam patterns of low frequency aperture arrays (LFAA) such as the Murchison Widefield Array (MWA). This occurs when the observing wavelength becomes comparable to the dipole separation for LFAA tiles, which for the MWA occurs at approximately 300 MHz. The presence of these grating sidelobes has made calibration and image processing for 300 MHz MWA observations difficult. This work presents a new calibration and imaging strategy which employs existing techniques to process two example 300 MHz MWA observations. Observations are initially calibrated using a new 300 MHz sky-model which has been interpolated from low frequency and high frequency all-sky surveys. Using this 300 MHz model in conjunction with the accurate MWA tile primary beam model, we perform sky-model calibration for the two example observations. After initial calibration a self-calibration loop is performed by all-sky imaging each observation with WSCLEAN. Using the output all-sky image we mask the main lobe of the image. Using this masked image we perform a sky-subtraction by estimating the masked image visibilities using WSCLEAN. We then image the main lobe of the observations with WSCLEAN. This results in high dynamic range images of the two example observation main lobes. These images have a resolution of 2.4 arcminutes, with a maximum sensitivity of 31 mJy/beam. The calibration and imaging strategy demonstrated in this work, opens the door to performing science at 300 MHz with the MWA, which was previously an inaccessible domain. With this paper we release the code described below and the cross-matched catalogue along with the code to produce a sky-model in the range 70-1400 MHz.
△ Less
Submitted 2 April, 2021;
originally announced April 2021.
-
Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations
Authors:
Cathryn M. Trott,
C. H. Jordan,
S. Midgley,
N. Barry,
B. Greig,
B. Pindor,
J. H. Cook,
G. Sleap,
S. J. Tingay,
D. Ung,
P. Hancock,
A. Williams,
J. Bowman,
R. Byrne,
A. Chokshi,
B. J. Hazelton,
K. Hasegawa,
D. Jacobs,
R. C. Joseph,
W. Li,
J. L. B Line,
C. Lynch,
B. McKinley,
D. A. Mitchell,
M. F. Morales
, et al. (11 additional authors not shown)
Abstract:
We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8…
▽ More
We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$σ$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.
△ Less
Submitted 6 February, 2020;
originally announced February 2020.
