-
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.
-
Aliasing from Galactic Plane Setting in Widefield Radio Interferometry
Authors:
N. Barry,
J. L. B. Line,
C. R. Lynch,
M. Kriele,
J. Cook
Abstract:
Measurements with widefield radio interferometers often include the near-infinite gradient between the sky and the horizon. This causes aliasing inherent to the measurement itself, and is purely a consequence of the Fourier basis. For this reason, the horizon is often attenuated by the instrumental beam down to levels deemed inconsequential. However, this effect is enhanced via our own Galactic pl…
▽ More
Measurements with widefield radio interferometers often include the near-infinite gradient between the sky and the horizon. This causes aliasing inherent to the measurement itself, and is purely a consequence of the Fourier basis. For this reason, the horizon is often attenuated by the instrumental beam down to levels deemed inconsequential. However, this effect is enhanced via our own Galactic plane as it sets over the course of a night. We show all-sky simulations of the Galactic plane setting in a radio interferometer in detail for the first time. We then apply these simulations to the Murchison Widefield Array to show that a beam attenuation of 0.1% is not sufficient in some precision science cases. We determine that the noise statistics of the residual data image are drastically more Gaussian with aliasing removal, and explore consequences in simulation for cataloging of extragalactic sources and 21 cm Epoch of Reionization detection via the power spectrum.
△ Less
Submitted 4 April, 2024; v1 submitted 12 December, 2023;
originally announced December 2023.
-
Multi-frequency angular power spectrum of the 21~cm signal from the Epoch of Reionisation using the Murchison Widefield Array
Authors:
Cathryn M. Trott,
Rajesh Mondal,
Garrelt Mellema,
Steven G. Murray,
Bradley Greig,
Jack L. B. Line,
Nichole Barry,
Miguel F. Morales
Abstract:
The Multi-frequency Angular Power Spectrum (MAPS) is an alternative to spherically-averaged power spectra, and computes local fluctuations in the angular power spectrum without need for line-of-sight spectral transform. To test different approaches to MAPS and treatment of the foreground contamination, and compare with the spherically-averaged power spectrum, and the single-frequency angular power…
▽ More
The Multi-frequency Angular Power Spectrum (MAPS) is an alternative to spherically-averaged power spectra, and computes local fluctuations in the angular power spectrum without need for line-of-sight spectral transform. To test different approaches to MAPS and treatment of the foreground contamination, and compare with the spherically-averaged power spectrum, and the single-frequency angular power spectrum. We apply the MAPS to 110~hours of data in $z=6.2-7.5$ obtained for the Murchison Widefield Array Epoch of Reionisation experiment to compute the statistical power of 21~cm brightness temperature fluctuations. In the presence of bright foregrounds, a filter is applied to remove large-scale modes prior to MAPS application, significantly reducing MAPS power due to systematics. The MAPS shows a contrast of 10$^2$--10$^3$ to a simulated 21~cm cosmological signal for spectral separations of 0--4~MHz after application of the filter, reflecting results for the spherically-averaged power spectrum. The single-frequency angular power spectrum is also computed. At $z=7.5$ and $l=200$, we find an angular power of 53~mK$^2$, exceeding a simulated cosmological signal power by a factor of one thousand. Residual spectral structure, inherent to the calibrated data, and not spectral leakage from large-scale modes, is the dominant source of systematic power bias. The single-frequency angular power spectrum yields slightly poorer results compared with the spherically-averaged power spectrum, having applied a spectral filter to reduce foregrounds. Exploration of other filters may improve this result, along with consideration of wider bandwidths.
△ Less
Submitted 11 August, 2022;
originally announced August 2022.
-
Optimising MWA EoR data processing for improved 21 cm power spectrum measurements -- fine-tuning ionospheric corrections
Authors:
J. Kariuki Chege,
C. H. Jordan,
C. Lynch,
C. M. Trott,
J. L. B. Line,
B. Pindor,
S. Yoshiura
Abstract:
The redshifted cosmological 21 cm signal emitted by neutral hydrogen during the first billion years of the universe is much fainter relative to other galactic and extragalactic radio emissions, posing a great challenge towards detection of the signal. Therefore, precise instrumental calibration is a vital prerequisite for the success of radio interferometers such as the Murchison Widefield Array (…
▽ More
The redshifted cosmological 21 cm signal emitted by neutral hydrogen during the first billion years of the universe is much fainter relative to other galactic and extragalactic radio emissions, posing a great challenge towards detection of the signal. Therefore, precise instrumental calibration is a vital prerequisite for the success of radio interferometers such as the Murchison Widefield Array (MWA), which aim for a 21 cm detection. Over the previous years, novel calibration techniques targeting the power spectrum paradigm of EoR science have been actively researched and where possible implemented. Using recently acquired computation resources for the MWA, we test the full capabilities of the state-of-the-art calibration techniques available for the MWA EoR project, with a focus on both direction dependent and direction independent calibration. Specifically, we investigate improvements that can be made in the vital calibration stages of sky modelling, ionospheric correction, and compact source foreground subtraction as applied in the hybrid foreground mitigation approach (one that combines both foreground subtraction and avoidance). Additionally, we investigate a method of ionospheric correction using interpolated ionospheric phase screens and assess its performance in the power spectrum space. Overall, we identify a refined RTS calibration configuration that leads to an at least 2 factor reduction of the EoR window power contamination at the $0.1 \; \text{hMpc}^{-1}$ scale. The improvement marks a step further towards detecting the 21 cm signal using the MWA and the forthcoming SKA low telescope.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
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.
-
The MWA Long Baseline Epoch of Reionisation Survey: I. Improved Source Catalogue for the EoR 0 field
Authors:
C. R. Lynch,
T. J. Galvin,
J. L. B. Line,
C. H. Jordan,
C. M. Trott,
J. K. Chege,
B. McKinley,
M. Johnston-Hollitt,
S. J. Tingay
Abstract:
One of the principal systematic constraints on the Epoch of Reionisation (EoR) experiment is the accuracy of the foreground calibration model. Recent results have shown that highly accurate models of extended foreground sources, and including models for sources in both the primary beam and its sidelobes, are necessary for reducing foreground power. To improve the accuracy of the source models for…
▽ More
One of the principal systematic constraints on the Epoch of Reionisation (EoR) experiment is the accuracy of the foreground calibration model. Recent results have shown that highly accurate models of extended foreground sources, and including models for sources in both the primary beam and its sidelobes, are necessary for reducing foreground power. To improve the accuracy of the source models for the EoR fields observed by the Murchison Widefield Array (MWA), we conducted the MWA Long Baseline Epoch of Reionisation Survey (LoBES). This survey consists of multi-frequency observations of the main MWA EoR fields and their eight neighbouring fields using the MWA Phase II extended array. We present the results of the first half of this survey centred on the MWA EoR0 observing field (centred at RA(J2000) 0 h, Dec(J2000) -27 deg). This half of the survey covers an area of 3069 degrees$^2$, with an average rms of 2.1 mJy beam$^{-1}$. The resulting catalogue contains a total of 80824 sources, with 16 separate spectral measurements between 100 and 230 MHz, and spectral modelling for 78$\%$ of these sources. Over this region we estimate that the catalogue is 90$\%$ complete at 32 mJy, and 70$\%$ complete at 10.5~mJy. The overall normalised source counts are found to be in good agreement with previous low-frequency surveys at similar sensitivities. Testing the performance of the new source models we measure lower residual rms values for peeled sources, particularly for extended sources, in a set of MWA Phase I data. The 2-dimensional power spectrum of these data residuals also show improvement on small angular scales -- consistent with the better angular resolution of the LoBES catalogue. It is clear that the LoBES sky models improve upon the current sky model used by the Australian MWA EoR group for the EoR0 field.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
Epoch of Reionization Power Spectrum Limits from Murchison Widefield Array Data Targeted at EoR1 Field
Authors:
M. Rahimi,
B. Pindor,
J. L. B. Line,
N. Barry,
C. M. Trott,
R. L. Webster,
C. H. Jordan,
M. Wilensky,
S. Yoshiura,
A. Beardsley,
J. Bowman,
R. Byrne,
A. Chokshi,
B. J. Hazelton,
K. Hasegawa,
E. Howard,
B. Greig,
D. Jacobs,
R. Joseph,
M. Kolopanis,
C. Lynch,
B. McKinley,
D. A. Mitchell,
S. Murray,
M. F. Morales
, et al. (6 additional authors not shown)
Abstract:
Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is ins…
▽ More
Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is instructive to explore different aspects of existing datasets. One such aspect is the choice of observing field. To date, MWA EoR observations have largely focused on the EoR0 field. In this work, we present a new detailed analysis of the EoR1 field. The EoR1 field is one of the coldest regions of the Southern radio sky, but contains the very bright radio galaxy Fornax-A. The presence of this bright extended source in the primary beam of the interferometer makes the calibration and analysis of EoR1 particularly challenging. We demonstrate the effectiveness of a recently developed shapelet model of Fornax-A in improving the results from this field. We also describe and apply a series of data quality metrics which identify and remove systematically contaminated data. With substantially improved source models, upgraded analysis algorithms and enhanced data quality metrics, we determine EoR power spectrum upper limits based on analysis of the best $\sim$14-hours data observed during 2015 and 2014 at redshifts 6.5, 6.8 and 7.1, with the lowest $2σ$ upper limit at z=6.5 of $Δ^2 \leq (73.78 ~\mathrm{mK)^2}$ at $k=0.13~\mathrm{h~ Mpc^{-1}}$, improving on previous EoR1 measurement results.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Constraining the 21cm brightness temperature of the IGM at $z$=6.6 around LAEs with the Murchison Widefield Array
Authors:
Cathryn M. Trott,
C. H. Jordan,
J. L. B. Line,
C. R. Lynch,
S. Yoshiura,
B. McKinley,
P. Dayal,
B. Pindor,
A. Hutter,
K. Takahashi,
R. B. Wayth,
N. Barry,
A. Beardsley,
J. Bowman,
R. Byrne,
A. Chokshi,
B. Greig,
K. Hasegawa,
B. J. Hazelton,
E. Howard,
D. Jacobs,
M. Kolopanis,
D. A. Mitchell,
M. F. Morales,
S. Murray
, et al. (7 additional authors not shown)
Abstract:
The locations of Ly-$α$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for t…
▽ More
The locations of Ly-$α$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for the spatial distribution of the gas, thereby providing prior information unavailable to power spectrum measurements. We use a test set of 12 hours of observations from the Murchison Widefield Array (MWA) in extended array configuration, to constrain the neutral hydrogen signature of 58 LAEs, detected with the Subaru Hypersuprime Cam in the \textit{Silverrush} survey, centred on $z$=6.58. We assume that detectable emitters reside in the centre of ionised HII bubbles during the end of reionization, and predict the redshifted neutral hydrogen signal corresponding to the remaining neutral regions using a set of different ionised bubble radii. A prewhitening matched filter detector is introduced to assess detectability. We demonstrate the ability to detect, or place limits upon, the amplitude of brightness temperature fluctuations, and the characteristic HII bubble size. With our limited data, we constrain the brightness temperature of neutral hydrogen to $Δ{\rm T}_B<$30 mK ($<$200 mK) at 95% (99%) confidence for lognormally-distributed bubbles of radii, $R_B =$ 15$\pm$2$h^{-1}$cMpc.
△ Less
Submitted 30 July, 2021;
originally announced July 2021.
-
A Map of Diffuse Radio Emission at 182 MHz to Enhance Epoch of Reionization Observations in the Southern Hemisphere
Authors:
Ruby Byrne,
Miguel F. Morales,
Bryna Hazelton,
Ian Sullivan,
Nichole Barry,
Christene Lynch,
Jack L. B. Line,
Daniel C. Jacobs
Abstract:
We present a broadband map of polarized diffuse emission at 167-198 MHz developed from data from the Murchison Widefield Array (MWA). The map is designed to improve visibility simulation and precision calibration for 21 cm Epoch of Reionization (EoR) experiments. It covers a large swath - 11,000 sq. deg. - of the Southern Hemisphere sky in all four Stokes parameters and captures emission on angula…
▽ More
We present a broadband map of polarized diffuse emission at 167-198 MHz developed from data from the Murchison Widefield Array (MWA). The map is designed to improve visibility simulation and precision calibration for 21 cm Epoch of Reionization (EoR) experiments. It covers a large swath - 11,000 sq. deg. - of the Southern Hemisphere sky in all four Stokes parameters and captures emission on angular scales of 1 to 9 degrees. The band-averaged diffuse structure is predominantly unpolarized but has significant linearly polarized structure near RA = 0 h. We evaluate the accuracy of the map by combining it with the GLEAM catalog and simulating an observation from the MWA, demonstrating that the accuracy of the short baselines (6.1-50 wavelengths) now approaches the accuracy of the longer baselines typically used for EoR calibration. We discuss how to use the map for visibility simulation for a variety of interferometric arrays. The map has potential to improve calibration accuracy for experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the forthcoming Square Kilometre Array (SKA) as well as the MWA.
△ Less
Submitted 15 November, 2021; v1 submitted 23 July, 2021;
originally announced July 2021.
-
A new MWA limit on the 21 cm Power Spectrum at Redshifts $\sim$ 13 $-$ 17
Authors:
S. Yoshiura,
B. Pindor,
J. L. B. Line,
N. Barry,
C. M. Trott,
A. Beardsley,
J. Bowman,
R. Byrne,
A. Chokshi,
B. J. Hazelton,
K. Hasegawa,
E. Howard,
B. Greig,
D. Jacobs,
C. H. Jordan,
R. Joseph,
M. Kolopanis,
C. Lynch,
B. McKinley,
D. A. Mitchell,
M. F. Morales,
S. G. Murray,
J. C. Pober,
M. Rahimi,
K. Takahashi
, et al. (7 additional authors not shown)
Abstract:
Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction indepen…
▽ More
Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration.
△ Less
Submitted 26 May, 2021;
originally announced May 2021.
-
Simulations of ionospheric refraction on radio interferometric data
Authors:
J. Kariuki Chege,
C. H. Jordan,
C. Lynch,
J. L. B. Line,
C. M. Trott
Abstract:
The Epoch of Reionisation (EoR) is the period within which the neutral universe transitioned to an ionised one. This period remains unobserved using low-frequency radio interferometers which target the 21 cm signal of neutral hydrogen emitted in this era. The Murchison Widefield Array (MWA) radio telescope was built with the detection of this signal as one of its major science goals. One of the mo…
▽ More
The Epoch of Reionisation (EoR) is the period within which the neutral universe transitioned to an ionised one. This period remains unobserved using low-frequency radio interferometers which target the 21 cm signal of neutral hydrogen emitted in this era. The Murchison Widefield Array (MWA) radio telescope was built with the detection of this signal as one of its major science goals. One of the most significant challenges towards a successful detection is that of calibration, especially in the presence of the Earth's ionosphere. By introducing refractive source shifts, distorting source shapes and scintillating flux densities, the ionosphere is a major nuisance in low-frequency radio astronomy. We introduce SIVIO, a software tool developed for simulating observations of the MWA through different ionospheric conditions estimated using thin screen approximation models and propagated into the visibilities. This enables us to directly assess the impact of the ionosphere on observed EoR data and the resulting power spectra. We show that the simulated data captures the dispersive behaviour of ionospheric effects. We show that the spatial structure of the simulated ionospheric media is accurately reconstructed either from the resultant source positional offsets or from parameters evaluated during the data calibration procedure. In turn, this will inform on the best strategies of identifying and efficiently eliminating ionospheric contamination in EoR data moving into the Square Kilometre Array era.
△ Less
Submitted 10 May, 2021;
originally announced May 2021.
-
Dual Polarization Measurements of MWA Beampatterns at 137 MHz
Authors:
A. Chokshi,
J. L. B. Line,
N. Barry,
D. Ung,
D. Kenney,
A. McPhail,
A. Williams,
R. L. Webster
Abstract:
The wide adoption of low-frequency radio interferometers as a tool for deeper and higher resolution astronomical observations has revolutionised radio astronomy. Despite their construction from static, relatively simple dipoles, the sheer number of distinct elements introduces new, complicated instrumental effects. Their necessary remote locations exacerbate failure rates, while electronic interac…
▽ More
The wide adoption of low-frequency radio interferometers as a tool for deeper and higher resolution astronomical observations has revolutionised radio astronomy. Despite their construction from static, relatively simple dipoles, the sheer number of distinct elements introduces new, complicated instrumental effects. Their necessary remote locations exacerbate failure rates, while electronic interactions between the many adjacent receiving elements can lead to non-trivial instrumental effects. The Murchison Widefield Array (MWA) employs phased array antenna elements (tiles), which improve collecting area at the expense of complex beam shapes. Advanced electromagnetic simulations have produced the Fully Embedded Element (FEE) simulated beam model which has been highly successful in describing the ideal beam response of MWA antennas. This work focuses on the relatively unexplored aspect of various in-situ, environmental perturbations to beam models and represents the first large-scale, in-situ, all-sky measurement of MWA beam shapes at multiple polarizations and pointings. Our satellite based beam measurement approach enables all-sky beam response measurements with a dynamic range of \sim 50 dB, at 137 MHz.
△ Less
Submitted 31 January, 2021;
originally announced February 2021.
-
The Impact of Tandem Redundant/Sky-Based Calibration in MWA Phase II Data Analysis
Authors:
Zheng Zhang,
Jonathan C. Pober,
Wenyang Li,
Bryna J. Hazelton,
Miguel F. Morales,
Cathryn M. Trott,
Christopher H. Jordan,
Ronniy C. Joseph,
Adam Beardsley,
Nichole Barry,
Ruby Byrne,
Steven J. Tingay,
Aman Chokshi,
Kenji Hasegawa,
Daniel C. Jacobs,
Adam Lanman,
Jack L. B. Line,
Christene Lynch,
Benjamin McKinley,
Daniel A. Mitchell,
Steven Murray,
Bart Pindor,
Mahsa Rahimi,
Keitaro Takahashi,
Randall B. Wayth
, et al. (4 additional authors not shown)
Abstract:
Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR)…
▽ More
Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in \cite{Li_2018} and \cite{Wenyang_2019} studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the power spectrum from tandem calibration are significant. To understand this result, we analyze both the calibration solutions themselves and the effects on the power spectrum over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model-incompleteness error.
△ Less
Submitted 19 September, 2020;
originally announced September 2020.
-
Murchison Widefield Array detection of steep-spectrum, diffuse, non-thermal radio emission within Abell 1127
Authors:
Stefan W. Duchesne,
Melanie Johnston-Hollitt,
Zhenghao Zhu,
Randall B. Wayth,
Jack L. B. Line
Abstract:
Diffuse, non-thermal emission in galaxy clusters is increasingly being detected in low-frequency radio surveys and images. We present a new diffuse, steep-spectrum, non-thermal radio source within the cluster Abell 1127 found in survey data from the Murchison Widefield Array (MWA). We perform follow-up observations with the 'extended' configuration MWA Phase II with improved resolution to better r…
▽ More
Diffuse, non-thermal emission in galaxy clusters is increasingly being detected in low-frequency radio surveys and images. We present a new diffuse, steep-spectrum, non-thermal radio source within the cluster Abell 1127 found in survey data from the Murchison Widefield Array (MWA). We perform follow-up observations with the 'extended' configuration MWA Phase II with improved resolution to better resolve the source and measure its low-frequency spectral properties. We use archival Very Large Array S-band data to remove the discrete source contribution from the MWA data, and from a power law model fit we find a spectral index of $-1.83\pm0.29$ broadly consistent with relic-type sources. The source is revealed by the Giant Metrewave Radio Telescope (GMRT) at 150 MHz to have an elongated morphology, with a projected linear size of 850 kpc as measured in the MWA data. Using Chandra observations we derive morphological estimators and confirm quantitatively that the cluster is in a disturbed dynamical state, consistent with the majority of phoenices and relics being hosted by merging clusters. We discuss the implications of relying on morphology and low-resolution imaging alone for the classification of such sources and highlight the usefulness of the MHz to GHz radio spectrum in classifying these types of emission. Finally, we discuss the benefits and limitations of using the MWA Phase II in conjunction with other instruments for detailed studies of diffuse, steep-spectrum, non-thermal radio emission within galaxy clusters.
△ Less
Submitted 29 July, 2020;
originally announced July 2020.
-
Modelling and peeling extended sources with shapelets: a Fornax A case study
Authors:
J. L. B. Line,
D. A. Mitchell,
B. Pindor,
J. L. Riding,
B. McKinley,
R. L. Webster,
C. M. Trott,
N. Hurley-Walker,
A. R. Offringa
Abstract:
To make a power spectrum (PS) detection of the 21 cm signal from the Epoch of Reionisation (EoR), one must avoid/subtract bright foreground sources. Sources such as Fornax A present a modelling challenge due to spatial structures spanning from arc seconds up to a degree. We compare modelling with multi-scale (MS) CLEAN components to 'shapelets', an alternative set of basis functions. We introduce…
▽ More
To make a power spectrum (PS) detection of the 21 cm signal from the Epoch of Reionisation (EoR), one must avoid/subtract bright foreground sources. Sources such as Fornax A present a modelling challenge due to spatial structures spanning from arc seconds up to a degree. We compare modelling with multi-scale (MS) CLEAN components to 'shapelets', an alternative set of basis functions. We introduce a new image-based shapelet modelling package, SHAMFI. We also introduce a new CUDA simulation code (WODEN) to generate point source, Gaussian, and shapelet components into visibilities. We test performance by modelling a simulation of Fornax A, peeling the model from simulated visibilities, and producing a residual PS. We find the shapelet method consistently subtracts large-angular-scale emission well, even when the angular-resolution of the data is changed. We find that when increasing the angular-resolution of the data, the MS CLEAN model worsens at large angular-scales. When testing on real MWA data, the expected improvement is not seen in real data because of the other dominating systematics still present. Through further simulation we find the expected differences to be lower than obtainable through current processing pipelines. We conclude shapelets are worthwhile for subtracting extended galaxies, and may prove essential for an EoR detection in the future, once other systematics have been addressed.
△ Less
Submitted 19 May, 2020;
originally announced May 2020.
-
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.
-
Measuring the global 21-cm signal with the MWA-I: improved measurements of the Galactic synchrotron background using lunar occultation
Authors:
B. McKinley,
G. Bernardi,
C. M. Trott,
J. L. B. Line,
R. B. Wayth,
A. R. Offringa,
B. Pindor,
C. H. Jordan,
M. Sokolowski,
S. J. Tingay,
E. Lenc,
N. Hurley-Walker,
J. D. Bowman,
F. Briggs,
R. L. Webster
Abstract:
We present early results from a project to measure the sky-averaged (global), redshifted $21\,$cm signal from the Epoch of Reionisation (EoR), using the Murchison Widefield Array (MWA) telescope. Because interferometers are not sensitive to a spatially-invariant global average, they cannot be used to detect this signal using standard techniques. However, lunar occultation of the radio sky imprints…
▽ More
We present early results from a project to measure the sky-averaged (global), redshifted $21\,$cm signal from the Epoch of Reionisation (EoR), using the Murchison Widefield Array (MWA) telescope. Because interferometers are not sensitive to a spatially-invariant global average, they cannot be used to detect this signal using standard techniques. However, lunar occultation of the radio sky imprints a spatial structure on the global signal, allowing us to measure the average brightness temperature of the patch of sky immediately surrounding the Moon. In this paper we present one night of Moon observations with the MWA between 72 - 230 MHz and verify our techniques to extract the background sky temperature from measurements of the Moon's flux density. We improve upon previous work using the lunar occultation technique by using a more sophisticated model for reflected `earthshine' and by employing image differencing to remove imaging artefacts. We leave the Moon's (constant) radio brightness temperature as a free parameter in our fit to the data and as a result, measure $T_{\rm{moon}} = 180 \pm 12 $ K and a Galactic synchrotron spectral index of $-2.64\pm0.14$, at the position of the Moon. Finally, we evaluate the prospects of the lunar occultation technique for a global EoR detection and map out a way forward for future work with the MWA.
△ Less
Submitted 5 September, 2018;
originally announced September 2018.
-
In situ measurement of MWA primary beam variation using ORBCOMM
Authors:
J. L. B. Line,
B. McKinley,
J. Rasti,
M. Bhardwaj,
R. B. Wayth,
R. L. Webster,
D. Ung,
D. Emrich,
L. Horsley,
A. Beardsley,
B. Crosse,
T. M. O. Franzen,
B. M. Gaensler,
M. Johnston-Hollitt,
D. L. Kaplan,
D. Kenney,
M. F. Morales,
D. Pallot,
K. Steele,
S. J. Tingay,
C. M. Trott,
M. Walker,
A. Williams,
C. Wu
Abstract:
We provide the first in situ measurements of antenna element (tile) beam shapes of the Murchison Widefield Array (MWA), a low radio-frequency interferometer and an SKA precursor. Most current MWA processing pipelines use an assumed beam shape, errors in which can cause absolute and relative flux density errors, as well as polarisation 'leakage'. This makes understanding the primary beam of paramou…
▽ More
We provide the first in situ measurements of antenna element (tile) beam shapes of the Murchison Widefield Array (MWA), a low radio-frequency interferometer and an SKA precursor. Most current MWA processing pipelines use an assumed beam shape, errors in which can cause absolute and relative flux density errors, as well as polarisation 'leakage'. This makes understanding the primary beam of paramount importance, especially for sensitive experiments such as a measurement of the 21 cm line from the epoch of reionisation (EoR). The calibration requirements for measuring the EoR 21 cm line are so extreme that tile to tile beam variations may affect our ability to make a detection. Measuring the primary beam shape from visibilities alone is challenging, as multiple instrumental, atmospheric, and astrophysical factors contribute to uncertainties in the data. Building on the methods of Neben et al. (2015), we tap directly into the receiving elements of the MWA before any digitisation or correlation of the signal. Using ORBCOMM satellite passes we are able to produce all-sky maps for 4 separate tiles in the XX polarisation. We find good agreement with the cutting-edge 'fully' embedded element (FEE) model of Sokolowski et al. (2017), and observe that the MWA beamformers consistently recreate beam shapes to within ~1dB in the reliable areas of our beam maps. We also clearly observe the effects of a missing dipole from a tile in one of our beam maps, and show that the FEE model is able to reproduce this modified beam shape. We end by motivating and outlining additional onsite experiments to further constrain the primary beam behaviour.
△ Less
Submitted 13 August, 2018;
originally announced August 2018.
-
Detectability of 21cm-signal during the Epoch of Reionization with 21cm-Lyman-α emitter cross-correlation. II. Foreground contamination
Authors:
S. Yoshiura,
J. L. B. Line,
K. Kubota,
K. Hasegawa,
K. Takahashi
Abstract:
Cross-correlation between the redshifted 21 cm signal and Lyman-α emitters (LAEs) is powerful tool to probe the Epoch of Reionization (EoR). Although the cross-power spectrum (PS) has an advantage of not correlating with foregrounds much brighter than the 21 cm signal, the galactic and extra-galactic foregrounds prevent detection since they contribute to the variance of the cross PS. Therefore, st…
▽ More
Cross-correlation between the redshifted 21 cm signal and Lyman-α emitters (LAEs) is powerful tool to probe the Epoch of Reionization (EoR). Although the cross-power spectrum (PS) has an advantage of not correlating with foregrounds much brighter than the 21 cm signal, the galactic and extra-galactic foregrounds prevent detection since they contribute to the variance of the cross PS. Therefore, strategies for mitigating foregrounds are required. In this work, we study the impact of foreground avoidance on the measurement of the 21 cm-LAE cross-correlation. We then simulate the 21 cm observation as observed by the Murchison Widefield Array (MWA). The point source foreground is modelled from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey catalogue, and the diffuse foreground is evaluated using a parametric model. For LAE observations, we assume a large survey of the Subaru Hyper Supreme-Cam (HSC), with spectroscopic observations of the Prime Focus Spectrograph (PFS). To predict the 21 cm signal, we employ a numerical simulation combining post processed radiative transfer and radiation hydrodynamics. Using these models, the signal-to-noise ratio of 2D PS shows the foreground contamination dominates the error of cross-PS even in the so-called `EoR window'. We find that at least 99% of the point source foreground and 80% of the galactic diffuse foreground must be removed to measure the EoR signal at large scales $k<0.5 h \rm Mpc^{-1}$. Additionally, a sensitivity 80 times larger than that of the MWA operating with 128 tiles and 99% of the point source foreground removal are required for a detection at small scales.
△ Less
Submitted 4 July, 2018; v1 submitted 13 September, 2017;
originally announced September 2017.
-
Extragalactic Peaked-Spectrum Radio Sources at Low Frequencies
Authors:
J. R. Callingham,
R. D. Ekers,
B. M. Gaensler,
J. L. B. Line,
N. Hurley-Walker,
E. M. Sadler,
S. J. Tingay,
P. J. Hancock,
M. E. Bell,
K. S. Dwarakanath,
B. -Q. For,
T. M. O. Franzen,
L. Hindson,
M. Johnston-Hollitt,
A. D. Kapinska,
E. Lenc,
B. McKinley,
J. Morgan,
A. R. Offringa,
P. Procopio,
L. Staveley-Smith,
R. B. Wayth,
C. Wu,
Q. Zheng
Abstract:
We present a sample of 1,483 sources that display spectral peaks between 72 MHz and 1.4 GHz, selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The GLEAM survey is the widest fractional bandwidth all-sky survey to date, ideal for identifying peaked-spectrum sources at low radio frequencies. Our peaked-spectrum sources are the low frequency analogues of g…
▽ More
We present a sample of 1,483 sources that display spectral peaks between 72 MHz and 1.4 GHz, selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The GLEAM survey is the widest fractional bandwidth all-sky survey to date, ideal for identifying peaked-spectrum sources at low radio frequencies. Our peaked-spectrum sources are the low frequency analogues of gigahertz-peaked spectrum (GPS) and compact-steep spectrum (CSS) sources, which have been hypothesized to be the precursors to massive radio galaxies. Our sample more than doubles the number of known peaked-spectrum candidates, and 95% of our sample have a newly characterized spectral peak. We highlight that some GPS sources peaking above 5 GHz have had multiple epochs of nuclear activity, and demonstrate the possibility of identifying high redshift ($z > 2$) galaxies via steep optically thin spectral indices and low observed peak frequencies. The distribution of the optically thick spectral indices of our sample is consistent with past GPS/CSS samples but with a large dispersion, suggesting that the spectral peak is a product of an inhomogeneous environment that is individualistic. We find no dependence of observed peak frequency with redshift, consistent with the peaked-spectrum sample comprising both local CSS sources and high-redshift GPS sources. The 5 GHz luminosity distribution lacks the brightest GPS and CSS sources of previous samples, implying that a convolution of source evolution and redshift influences the type of peaked-spectrum sources identified below 1 GHz. Finally, we discuss sources with optically thick spectral indices that exceed the synchrotron self-absorption limit.
△ Less
Submitted 10 January, 2017;
originally announced January 2017.
-
PUMA: The Positional Update and Matching Algorithm
Authors:
J. L. B. Line,
R. L. Webster,
B. Pindor,
D. A. Mitchell,
C. M. Trott
Abstract:
We present new software to cross-match low-frequency radio catalogues: the Positional Update and Matching Algorithm (PUMA). PUMA combines a positional Bayesian probabilistic approach with spectral matching criteria, allowing for confusing sources in the matching process. We go on to create a radio sky model using PUMA based on the Murchison Widefield Array Commissioning Survey, and are able to aut…
▽ More
We present new software to cross-match low-frequency radio catalogues: the Positional Update and Matching Algorithm (PUMA). PUMA combines a positional Bayesian probabilistic approach with spectral matching criteria, allowing for confusing sources in the matching process. We go on to create a radio sky model using PUMA based on the Murchison Widefield Array Commissioning Survey, and are able to automatically cross-match 98.5% of sources. Using the characteristics of this sky model, we create simple simulated mock catalogues on which to test PUMA, and find that PUMA can reliably find the correct spectral indices of sources, along with being able to recover ionospheric offsets. Finally, we use this sky model to calibrate and remove foreground sources from simulated interferometric data, generated using OSKAR (the Oxford University visibility generator). We demonstrate that there is a substantial improvement in foreground source removal when using higher frequency and higher resolution source positions, even when correcting positions by an average of 0.3 given a synthesized beam-width of 2.3.
△ Less
Submitted 16 November, 2016;
originally announced November 2016.
