Showing 1–50 of 104 results for author: Trott, M

MOSEL survey: Unwrapping the Epoch of Reionization through mimic galaxies at Cosmic Noon

Authors: Ravi Jaiswar, Anshu Gupta, Elisabete da Cunha, Cathryn M. Trott, Anishya Harshan, Andrew Battisti, Ben Forrest

Abstract: The nature of the first galaxies that reionized the universe during the Epoch of Reionization (EoR) remains unclear. Attempts to directly determine spectral properties of these early galaxies are affected by both limited photometric constraints across the spectrum and by the opacity of the intergalactic medium (IGM) to the Lyman Continuum (LyC) at high redshift. We approach this by analysing prope… ▽ More The nature of the first galaxies that reionized the universe during the Epoch of Reionization (EoR) remains unclear. Attempts to directly determine spectral properties of these early galaxies are affected by both limited photometric constraints across the spectrum and by the opacity of the intergalactic medium (IGM) to the Lyman Continuum (LyC) at high redshift. We approach this by analysing properties of analogous extreme emission line galaxies (EELGs, [OIII]+Hbeta EW $>400$ Angstrom) at $2.5<z<4$ from the ZFOURGE survey using the Multi-wavelength Analysis of Galaxy Physical Properties (MAGPHYS) SED fitting code. We compare these to galaxies at $z>5.5$ observed with the James Webb Space Telesope (JWST) with self-consistent spectral energy distribution fitting methodology. This work focuses on the comparison of their UV slopes ($β_P$), ionizing photon production efficiencies $ξ_{ion}$, star formation rates and dust properties to determine the effectiveness of this analogue selection technique. We report the median ionizing photon production efficiencies as log$_{10}(ξ_{ion}/(Hz\ {\rm erg}^{-1}))=$$25.14^{+0.06}_{-0.04}$,$25.16^{+0.06}_{-0.05}$,$25.16^{+0.04}_{-0.05}$,$25.18^{+0.06}_{-0.07}$ for our ZFOURGE control, ZFOURGE EELG, JADES and CEERS samples respectively. ZFOURGE EELGs are 0.57 dex lower in stellar mass and have half the dust extinction, compared to their ZFOURGE control counterparts. They also have a similar specific star formation rates and $β_P$ to the $z>5.5$ samples. We find that EELGs at low redshift ($2.5<z<4$) are analogous to EoR galaxies in their dust attenuation and specific star formation rates. Their extensive photometric coverage and the accessibility of their LyC region opens pathways to infer stellar population properties in the EoR. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 12 pages, 8 figures, 2 tables, published in PASA

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.

Comments: 13 pages, 11 Figures, 1 table, accepted in the Publications of the Astronomical Society of Australia (PASA)

MOSEL survey: Spatially offset Lyman-continuum emission in a new emitter at z=3.088

Authors: Anshu Gupta, Cathryn M. Trott, Ravi Jaiswar, E. V. Ryan-Weber, Andrew J. Bunker, Ayan Acharyya, Alex J. Cameron, Ben Forrest, Glenn G. Kacprzak, Themiya Nanayakkara, Kim-Vy Tran, Aman Chokshi

Abstract: We present the discovery of a unique Lyman-continuum (LyC) emitter at z=3.088. The LyC emission were detected using the Hubble Space Telescope (HST) WFC3/UVIS F336W filter, covering a rest-frame wavelength range of 760-900 Angstrom. The peak signal-to-noise ratio (SNR) of LyC emission is 3.9 in a r=0.24'' aperture and is spatially offset by 0.29''+/-0.04'' (~ 2.2+/-0.3 kpc) from the rest-UV emissi… ▽ More We present the discovery of a unique Lyman-continuum (LyC) emitter at z=3.088. The LyC emission were detected using the Hubble Space Telescope (HST) WFC3/UVIS F336W filter, covering a rest-frame wavelength range of 760-900 Angstrom. The peak signal-to-noise ratio (SNR) of LyC emission is 3.9 in a r=0.24'' aperture and is spatially offset by 0.29''+/-0.04'' (~ 2.2+/-0.3 kpc) from the rest-UV emission peak (F606W). By combining imaging and spectroscopic data from the James Webb Space Telescope (JWST) JADES, FRESCO and JEMS surveys, along with VLT/MUSE data from the MXDF survey, we estimate that the probability of random alignment with an interloper galaxy causing the LyC emission is less than 6x10^-5. The interstellar medium (ISM) conditions in the galaxy are similar to other LyC emitters at high redshift (12+log(O/H)=7.79+/-0.06, logU =-3.27+/-0.14, O32 = 3.65+/-0.22), although the single-peaked Lyman-alpha profile and lack of rest-UV emission lines suggest an optically thick ISM. We think that LyC photons are leaking through a narrow cone of optically thin neutral ISM, most likely created by a past merger (as evidenced by medium-band F210M and F182M images). Using the escape fraction constraints from individual leakers and a simple model, we estimate that the opening half-angle of ionization cones can be as low as 16^deg (2% ionised fraction) to reproduce some of the theoretical constraints on the average escape fraction for galaxies. The narrow opening angle required can explain the low number density of confirmed LyC leakers. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 12 pages, 8 figures, and 1 table. Submitted to The Astrophysical Journal

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.

Comments: 19 pages, 14 figures

Evidence of Ultra-faint Radio Frequency Interference in Deep 21~cm Epoch of Reionization Power Spectra with the Murchison Widefield Array

Authors: Michael J. Wilensky, Miguel F. Morales, Bryna J. Hazelton, Pyxie L. Star, Nichole Barry, Ruby Byrne, C. H. Jordan, Daniel C. Jacobs, Jonathan C. Pober, C. M. Trott

Abstract: We present deep upper limits from the 2014 Murchison Widefield Array (MWA) Phase I observing season, with a particular emphasis on identifying the spectral fingerprints of extremely faint radio frequency interference (RFI) contamination in the 21~cm power spectra (PS). After meticulous RFI excision involving a combination of the \textsc{SSINS} RFI flagger and a series of PS-based jackknife tests,… ▽ More We present deep upper limits from the 2014 Murchison Widefield Array (MWA) Phase I observing season, with a particular emphasis on identifying the spectral fingerprints of extremely faint radio frequency interference (RFI) contamination in the 21~cm power spectra (PS). After meticulous RFI excision involving a combination of the \textsc{SSINS} RFI flagger and a series of PS-based jackknife tests, our lowest upper limit on the Epoch of Reionization (EoR) 21~cm PS signal is $Δ^2 \leq 1.61\cdot10^4 \text{ mK}^2$ at $k=0.258\text{ h Mpc}^{-1}$ at a redshift of 7.1 using 14.7 hours of data. By leveraging our understanding of how even fainter RFI is likely to contaminate the EoR PS, we are able to identify ultra-faint RFI signals in the cylindrical PS. Surprisingly this signature is most obvious in PS formed with less than an hour of data, but is potentially subdominant to other systematics in multiple-hour integrations. Since the total RFI budget in a PS detection is quite strict, this nontrivial integration behavior suggests a need to more realistically model coherently integrated ultra-faint RFI in PS measurements so that its potential contribution to a future detection can be diagnosed. △ Less

Submitted 7 November, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: Update acknowledgements, author metadata, and attach journal doi

Journal ref: The Astrophysical Journal, 2023, Volume 957, Number 2, p. 78

First constraints on Helium $^{+}{\rm He}^3$ evolution in $z=3-4$ using the 8.67GHz hyperfine transition

Authors: Cathryn M. Trott, Randall B. Wayth

Abstract: We present the first constraints on the cross-correlation power spectrum of HeII ($^{+}{\rm He}^3$) signal strength using the redshifted 8.67GHz hyperfine transition between $z=2.9$ and $z=4.1$ and with interferometric data obtained from the public archive of the Australia Telescope Compact Array. 210 hours of observations of the primary calibrator source B1934-638 were extracted from data obtaine… ▽ More We present the first constraints on the cross-correlation power spectrum of HeII ($^{+}{\rm He}^3$) signal strength using the redshifted 8.67GHz hyperfine transition between $z=2.9$ and $z=4.1$ and with interferometric data obtained from the public archive of the Australia Telescope Compact Array. 210 hours of observations of the primary calibrator source B1934-638 were extracted from data obtained with the telescope from 2014--2021, and coherently combined in a power spectrum pipeline to measure the HeII power across a range of spatial scales, and at three redshifts that span the period of Helium reionization. Our best limit places the brightness temperature fluctuation to be less than 557$μ$K on spatial scales of 30 arcmin at $z=2.91$, and less than 755$μ$K on scales of 30 arcmin at $z=4.14$ (2-sigma noise-limited). We measure a temperature of 489$μ$K at $z=2.91$. ATCA's few antennas and persistent remaining RFI in the data prevent deeper integrations improving the results. This work is a proof of principle to demonstrate how this type of experiment can be undertaken to reach the 0.01--1$μ$K level expected for the Helium signal at $z \sim 4$. △ Less

Submitted 20 November, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 10 pages, 4 figures, accepted in ApJ

Measuring the global 21-cm signal with the MWA-II: improved characterisation of lunar-reflected radio frequency interference

Authors: Himanshu Tiwari, Benjamin McKinley, Cathryn M. Trott, Nithyanandan Thyagarajan

Abstract: Radio interferometers can potentially detect the sky-averaged signal from the Cosmic Dawn (CD) and the Epoch of Reionisation (EoR) by studying the Moon as a thermal block to the foreground sky. The first step is to mitigate the Earth-based RFI reflections (Earthshine) from the Moon, which significantly contaminate the FM band $\approx 88-110$ MHz, crucial to CD-EoR science. We analysed MWA phase-I… ▽ More Radio interferometers can potentially detect the sky-averaged signal from the Cosmic Dawn (CD) and the Epoch of Reionisation (EoR) by studying the Moon as a thermal block to the foreground sky. The first step is to mitigate the Earth-based RFI reflections (Earthshine) from the Moon, which significantly contaminate the FM band $\approx 88-110$ MHz, crucial to CD-EoR science. We analysed MWA phase-I data from $72-180$ MHz at $40$ kHz resolution to understand the nature of Earthshine over three observing nights. We took two approaches to correct the Earthshine component from the Moon. In the first method, we mitigated the Earthshine using the flux density of the two components from the data, while in the second method, we used simulated flux density based on an FM catalogue to mitigate the Earthshine. Using these methods, we were able to recover the expected Galactic foreground temperature of the patch of sky obscured by the Moon. We performed a joint analysis of the Galactic foregrounds and the Moon's intrinsic temperature $(T_{\rm Moon})$ while assuming that the Moon has a constant thermal temperature throughout three epochs. We found $T_{\rm Moon}$ to be at $184.4\pm{2.6}\rm ~K$ and $173.8\pm{2.5}\rm ~K$ using the first and the second methods, respectively, and the best-fit values of the Galactic spectral index $(α)$ to be within the $5\%$ uncertainty level when compared with the global sky models. Compared with our previous work, these results improved constraints on the Galactic spectral index and the Moon's intrinsic temperature. We also simulated the Earthshine at MWA between November and December 2023 to find suitable observing times less affected by the Earthshine. Such observing windows act as Earthshine avoidance and can be used to perform future global CD-EoR experiments using the Moon with the MWA. △ Less

Submitted 1 November, 2023; v1 submitted 31 May, 2023; originally announced June 2023.

Comments: 17 pages, 14 figures and 8 tables, accepted for publication in PASA

Two-Screen Scattering in CRAFT FRBs

Authors: Mawson W. Sammons, Adam T. Deller, Marcin Glowacki, Kelly Gourdji, C. W. James, J. Xavier Prochaska, Hao Qiu, Danica R. Scott, R. M. Shannon, C. M. Trott

Abstract: Temporal broadening is a commonly observed property of fast radio bursts (FRBs), associated with turbulent media which cause radiowave scattering. Similarly to dispersion, scattering is an important probe of the media along the line of sight to an FRB source, such as the circum-burst or circum-galactic mediums (CGM). Measurements of characteristic scattering times alone are insufficient to constra… ▽ More Temporal broadening is a commonly observed property of fast radio bursts (FRBs), associated with turbulent media which cause radiowave scattering. Similarly to dispersion, scattering is an important probe of the media along the line of sight to an FRB source, such as the circum-burst or circum-galactic mediums (CGM). Measurements of characteristic scattering times alone are insufficient to constrain the position of the dominant scattering media along the line of sight. However, where more than one scattering screen exists, Galactic scintillation can be leveraged to form strong constraints. We quantify the scattering and scintillation in 10 FRBs with 1) known host galaxies and redshifts and 2) captured voltage data enabling high-time resolution analysis. We find strong evidence for two screens in three cases. For FRBs 20190608B and 20210320C, we find evidence for scattering screens less than approximately 16.7 and 3000 kpc respectively, from their sources, consistent with the scattering occurring in the circum-burst environment, the host ISM (inter-stellar medium) or the CGM. For FRB 20201124A we find a low modulation index that evolves over the burst's scattering tail, indicating the presence of a scattering screen $\approx9$ kpc from the host, and excluding the circum-burst environment from potential scattering sites. By assuming that pulse broadening is contributed by the host galaxy ISM or circum-burst environment, the lack of observed scintillation in four FRBs in our sample suggests that existing models may be poor estimators of scattering times associated with the Milky Way's ISM, similar to the anomalously low scattering observed for FRB 20201124A. △ Less

Submitted 25 August, 2023; v1 submitted 19 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures, 1 table. Accepted in MNRAS

Singular spectrum analysis of time series data from low frequency radiometers, with an application to SITARA data

Authors: Jishnu N. Thekkeppattu, Cathryn M. Trott, Benjamin McKinley

Abstract: Understanding the temporal characteristics of data from low frequency radio telescopes is of importance in devising suitable calibration strategies. Application of time series analysis techniques to data from radio telescopes can reveal a wealth of information that can aid in calibration. In this paper, we investigate singular spectrum analysis (SSA) as an analysis tool for radio data. We show the… ▽ More Understanding the temporal characteristics of data from low frequency radio telescopes is of importance in devising suitable calibration strategies. Application of time series analysis techniques to data from radio telescopes can reveal a wealth of information that can aid in calibration. In this paper, we investigate singular spectrum analysis (SSA) as an analysis tool for radio data. We show the intimate connection between SSA and Fourier techniques. We develop the relevant mathematics starting with an idealised periodic dataset and proceeding to include various non-ideal behaviours. We propose a novel technique to obtain long-term gain changes in data, leveraging the periodicity arising from sky drift through the antenna beams. We also simulate several plausible scenarios and apply the techniques to a 30-day time series data collected during June 2021 from SITARA - a short-spacing two element interferometer for global 21-cm detection. Applying the techniques to real data, we find that the first reconstructed component - the trend - has a strong anti-correlation with the local temperature suggesting temperature fluctuations as the most likely origin for the observed variations in the data. We also study the limitations of the calibration in the presence of diurnal gain variations and find that such variations are the likely impediment to calibrating SITARA data with SSA. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS

The Effect of Gravitational Lensing on Fast Transient Event Rates

Authors: Mawson W. Sammons, Clancy W. James, Cathryn M. Trott, Mark Walker

Abstract: Fast cosmological transients such as fast radio bursts (FRBs) and gamma-ray bursts (GRBs) represent a class of sources more compact than any other cosmological object. As such they are sensitive to significant magnification via gravitational lensing from a class of lenses which are not well-constrained by observations today. Low-mass primordial black holes are one such candidate which may constitu… ▽ More Fast cosmological transients such as fast radio bursts (FRBs) and gamma-ray bursts (GRBs) represent a class of sources more compact than any other cosmological object. As such they are sensitive to significant magnification via gravitational lensing from a class of lenses which are not well-constrained by observations today. Low-mass primordial black holes are one such candidate which may constitute a significant fraction of the Universe's dark matter. Current observations only constrain their density in the nearby Universe, giving fast transients from cosmological distances the potential to form complementary constraints. Motivated by this, we calculate the effect that gravitational lensing from a cosmological distribution of compact objects would have on the observed rates of FRBs and GRBs. For static lensing geometries, we rule out the prospect that all FRBs are gravitationally lensed for a range of lens masses and show that lens masses greater than $10^{-5}M_\odot$ can be constrained with 8000 un-localised high fluence FRBs at 1.4GHz, as might be detected by the next generation of FRB-finding telescopes. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 10 pages (+7 Supplementary). Accepted to MNRAS

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.

Comments: 9 pages, 9 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 666, A106 (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.

Comments: Accepted for publication by Publications of the Astronomical Society of Australia (PASA)

Generating extremely large-volume reionisation simulations

Authors: Bradley Greig, J. Stuart B. Wyithe, Steven G. Murray, Simon J. Mutch, Cathryn M. Trott

Abstract: Preparing for the first detection of the cosmic 21-cm signal from large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. To validate that these pipelines do not erroneously remove or add features that can mimic the cosmic signal (e.g. from side-lobes or large-scale power leakage), we require reionisation simulations larger than the experiment… ▽ More Preparing for the first detection of the cosmic 21-cm signal from large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. To validate that these pipelines do not erroneously remove or add features that can mimic the cosmic signal (e.g. from side-lobes or large-scale power leakage), we require reionisation simulations larger than the experiments primary field of view. For an experiment such as the MWA, with a field of view of $\sim25^{2}$ deg.$^{2}$, this would require a simulation of several Gpcs, which is currently infeasible. To overcome this, we developed a simplified version of the semi-numerical reionisation simulation code 21CMFAST preferencing large volumes over some physical accuracy by assuming linear theory for structure formation. With this, we constructed a 7.5 Gpc comoving volume with voxel resolution of $\sim1.17$ cMpc tailored specifically to the binned spectral resolution of the MWA. This simulation was used for validating the pipelines for the 2020 MWA 21-cm power spectrum (PS) upper limits (Trott et al.). We then use this large-volume simulation to explore: (i) whether smaller volume simulations are biased by the missing large-scale modes, (ii) non-Gaussianity in estimates of the cosmic variance, (iii) biases in the recovered 21-cm PS following foreground wedge removal and (iv) the impact of tiling smaller volume simulations to achieve extremely large volumes. In summary, we find: (i) no biases from missing large-scale power, (ii) significant contribution from non-Gaussianity in the cosmic variance as expected following Mondal et al. (iii) an over-estimate of the 21-cm PS of 10-20 per cent following wedge mode excision for our particular model and (iv) tiling smaller volume simulations under-estimates the large-scale power and also the estimated cosmic variance. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 13 pages, 9 figures and 1 table. Submitted to MNRAS, comments welcome

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.

Comments: 17 pages, 8 png figures

System design and calibration of SITARA -- a global 21 cm short spacing interferometer prototype

Authors: Jishnu N. Thekkeppattu, Benjamin McKinley, Cathryn M. Trott, Jake Jones, Daniel C. X. Ung

Abstract: Global 21-cm experiments require exquisitely precise calibration of the measurement systems in order to separate the weak 21-cm signal from Galactic and extragalactic foregrounds as well as instrumental systematics. Hitherto, experiments aiming to make this measurement have concentrated on measuring this signal using the single element approach. However, an alternative approach based on interferom… ▽ More Global 21-cm experiments require exquisitely precise calibration of the measurement systems in order to separate the weak 21-cm signal from Galactic and extragalactic foregrounds as well as instrumental systematics. Hitherto, experiments aiming to make this measurement have concentrated on measuring this signal using the single element approach. However, an alternative approach based on interferometers with short baselines is expected to alleviate some of the difficulties associated with a single element approach such as precision modelling of the receiver noise spectrum. Short spacing Interferometer Telescope probing cosmic dAwn and epoch of ReionisAtion (SITARA) is a short spacing interferometer deployed at the Murchison Radio-astronomy Observatory (MRO). It is intended to be a prototype or a test-bed to gain a better understanding of interferometry at short baselines, and develop tools to perform observations and data calibration. In this paper, we provide a description of the SITARA system and its deployment at the MRO, and discuss strategies developed to calibrate SITARA. We touch upon certain systematics seen in SITARA data and their modelling. We find that SITARA has sensitivity to all sky signals as well as non-negligible noise coupling between the antennas. It is seen that the coupled receiver noise has a spectral shape that broadly matches the theoretical calculations reported in prior works. We also find that when appropriately modified antenna radiation patterns taking into account the effects of mutual coupling are used, the measured data are well modelled by the standard visibility equation. △ Less

Submitted 20 March, 2022; originally announced March 2022.

Comments: 27 pages, 17 figures, Accepted for publication in PASA

Imaging the Southern Sky at 159MHz using Spherical Harmonics with the Engineering Development Array 2

Authors: Michael A. Kriele, Randall B. Wayth, Mark J. Bentum, Budi Juswardy, Cathryn M. Trott

Abstract: One of the major priorities of international radio astronomy is to study the early universe through the detection of the 21 cm HI line from the epoch of reionisation (EoR). Due to the weak nature of the 21 cm signal, an important part in the detection of the EoR is removing contaminating foregrounds from our observations as they are multiple orders of magnitude brighter. In order to achieve this,… ▽ More One of the major priorities of international radio astronomy is to study the early universe through the detection of the 21 cm HI line from the epoch of reionisation (EoR). Due to the weak nature of the 21 cm signal, an important part in the detection of the EoR is removing contaminating foregrounds from our observations as they are multiple orders of magnitude brighter. In order to achieve this, sky maps spanning a wide range of frequencies and angular scales are required for calibration and foreground subtraction. Complementing the existing low-frequency sky maps, we have constructed a Southern Sky map through spherical harmonic transit interferometry utilising the engineering development array 2 (EDA2), a square kilometre array (SKA) low-frequency array prototype system. We use the m-mode formalism to create an all-sky map at 159MHz with an angular resolution of 3 degrees, with data from the (EDA2) providing information over +60 degrees to -90 degrees in declination. We also introduce a new method for visualising and quantifying how the baseline distribution of an interferometer maps to the spherical harmonics, and discuss how prior information can be used to constrain spherical harmonic components that the interferometer is not sensitive to. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 29 pages, 35 figures, accepted in PASA (11-01-2022)

The impact of station far sidelobes on EoR/CD power spectra

Authors: Cathryn M. Trott

Abstract: Stations of dipole antennas for SKA1-Low will comprise 256 elements spread over an area with a diameter of 38~m. We consider the effect of residual unsubtracted sources well outside of the main beam for differing numbers of unique station configurations, in the Epoch of Reionisation (EoR) and the Cosmic Dawn (CD), both in simulation and with theoretical considerations. We find that beam sidelobes… ▽ More Stations of dipole antennas for SKA1-Low will comprise 256 elements spread over an area with a diameter of 38~m. We consider the effect of residual unsubtracted sources well outside of the main beam for differing numbers of unique station configurations, in the Epoch of Reionisation (EoR) and the Cosmic Dawn (CD), both in simulation and with theoretical considerations. We find that beam sidelobes imprint power that renders the cosmological signal unobservable over a range of scales unless compact sources are subtracted beyond θ_Z=30 degrees from zenith, and that station apodization will likely be required to control far sidelobes. An array with N_b=4 unique station configurations is sufficient to reduce the contamination, with an increase to N_b=8 showing little improvement. Comparison with an achromatic Airy disk beam model shows that beam sidelobe level is the main contributor to excess power in the EoR window, and beam chromaticity is less relevant. In the EoR, z=8.5, subtracting sources above 200~mJy out to θ_Z=45 degrees, will be required to access relevant modes of the power spectrum. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 20 pages, 11 figures, accepted for publication in SPIE JATIS SKA Special Edition

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.

Comments: 25 pages, 14 figures, Accepted for publication in PASA

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.

Comments: 18 pages, 17 figures, Accepted for publication in MNRAS

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.

Comments: 9 pages, 6 figures, accepted for publication in MNRAS

A broadband radio view of transient jet ejecta in the black hole candidate X-ray binary MAXI J1535-571

Authors: J. Chauhan, J. C. A. Miller-Jones, G. E. Anderson, A. Paduano, M. Sokolowski, C. Flynn, P. J. Hancock, N. Hurley-Walker, D. L. Kaplan, T. D. Russell, A. Bahramian, S. W. Duchesne, D. Altamirano, S. Croft, H. A. Krimm, G. R. Sivakoff, R. Soria, C. M. Trott, R. B. Wayth, V. Gupta, M. Johnston-Hollitt, S. J. Tingay

Abstract: We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535-571, which underwent a prolonged outburst beginning on 2 September 2017. We monitored MAXI J1535-571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 20 September to 14 October 2017. The source was quasi-simultaneously observed over… ▽ More We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535-571, which underwent a prolonged outburst beginning on 2 September 2017. We monitored MAXI J1535-571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 20 September to 14 October 2017. The source was quasi-simultaneously observed over the frequency range 0.84-19 GHz by UTMOST (the upgraded Molonglo Observatory Synthesis Telescope), the Australian Square Kilometre Array Pathfinder, the Australia Telescope Compact Array (ATCA), and the Australian Long Baseline Array (LBA). Using the LBA observations from 23 September 2017, we measured the source size to be $34\pm1$ mas. During the brightest radio flare on 21 September 2017, the source was detected down to 119 MHz by the MWA, and the radio spectrum indicates a turnover between 250 and 500 MHz, which is most likely due to synchrotron self-absorption (SSA). By fitting the radio spectrum with a SSA model and using the LBA size measurement, we determined various physical parameters of the jet knot (identified in ATCA data), including the jet opening angle (= $4.5\pm1.2^{\circ}$) and the magnetic field strength (= $104^{+80}_{-78}$ mG). Our fitted magnetic field strength agrees reasonably well with that inferred from the standard equipartition approach, suggesting the jet knot to be close to equipartition. Our study highlights the capabilities of the Australian suite of radio telescopes to jointly probe radio jets in black hole X-ray binaries (BH-XRBs) via simultaneous observations over a broad frequency range, and with differing angular resolutions. This suite allows us to determine the physical properties of XRB jets. Finally, our study emphasizes the potential contributions that can be made by the low-frequency part of the Square Kilometre Array (SKA-Low) in the study of BH-XRBs. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: This article has been accepted for publication in PASA. This article has 16 pages, 7 figures and 1 table

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.

Comments: 17 pages, 18 figures, accepted to MNRAS

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.

Comments: Accepted for publication in Publications of the Astronomical Society of Australia (PASA). 10 pages, 9 figures, 1 table

Journal ref: Publ. Astron. Soc. Aust. 38 (2021) e028

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.

Comments: 20 pages, 11 figures, accepted for publication in PASA

The All-Sky SignAl Short-Spacing INterferometer (ASSASSIN) I: Global sky measurements with the Engineering Development Array-2

Authors: B. McKinley, C. M. Trott, M. Sokolowski, R. B. Wayth, A. Sutinjo, N. Patra, J. Nambissan T., D. C. X. Ung

Abstract: Aiming to fill a crucial gap in our observational knowledge of the early Universe, experiments around the world continue to attempt to verify the claimed detection of the redshifted 21-cm signal from Cosmic Dawn by the EDGES experiment. This sky-averaged or 'global' signal from neutral hydrogen should be detectable at low radio frequencies (50-200 MHz), but is difficult to measure due to bright fo… ▽ More Aiming to fill a crucial gap in our observational knowledge of the early Universe, experiments around the world continue to attempt to verify the claimed detection of the redshifted 21-cm signal from Cosmic Dawn by the EDGES experiment. This sky-averaged or 'global' signal from neutral hydrogen should be detectable at low radio frequencies (50-200 MHz), but is difficult to measure due to bright foreground emission and difficulties in reaching the required levels of instrumental-calibration precision. In this paper we outline our progress toward using a novel new method to measure the global redshifted 21-cm signal. Motivated by the need to use alternative methods with very different systematic errors to EDGES for an independent result, we employ an array of closely-spaced antennas to measure the global sky signal interferometrically, rather than using the conventional approach with a single antenna. We use simulations to demonstrate our newly-developed methods and show that, for an idealised instrument, a 21-cm signal could theoretically be extracted from the visibilities of an array of closely-spaced dipoles. We verify that our signal-extraction methods work on real data using observations made with a Square Kilometre-Array-like prototype; the Engineering Development Array-2. Finally, we use the lessons learned in both our simulations and observations to lay out a clear plan for future work, which will ultimately lead to a new global redshifted 21-cm instrument: the All-Sky SignAl Short-Spacing INterferometer (ASSASSIN). △ Less

Submitted 13 September, 2020; originally announced September 2020.

Comments: 18 pages, 13 figures

Exploring reionisation and high-z galaxy observables with recent multi-redshift MWA upper limits on the 21-cm signal

Authors: Bradley Greig, Cathryn M. Trott, Nichole Barry, Simon J. Mutch, Bart Pindor, Rachel L. Webster, J. Stuart B. Wyithe

Abstract: We use the latest multi-redshift ($z=6.5-8.7$) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. These upper limits are achieved using 298 h of carefully excised data over four observing seasons. To explore these upper limits in the context of reionisation astrophysics, we use 21CMMC. We then connect the… ▽ More We use the latest multi-redshift ($z=6.5-8.7$) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. These upper limits are achieved using 298 h of carefully excised data over four observing seasons. To explore these upper limits in the context of reionisation astrophysics, we use 21CMMC. We then connect the disfavoured regions of parameter space to existing observational constraints on reionisation such as high-$z$ galaxy ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of models disfavoured by the MWA limits are already inconsistent with existing observational constraints. These inconsistent models arise from two classes of models: (i) `cold' reionisation and (ii) pure matter density fluctuations (i.e. no reionisation). However, a small subsample of models are consistent implying the existing MWA limits provide unique information in disfavouring models of reionisation, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding $1σ$ lower limits of $ε_{{\rm X},0.5-2~{\rm keV}}\gtrsim10^{34.5}$ erg s$^{-1}$ Mpc$^{-3}$. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, 2.4 K at $z=6.5, 6.8, 7.1, 7.8, 8.2, 8.7$. With this we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on $ε_{{\rm X},0.5-2~{\rm keV}}$ and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Comments: 14 pages, 7 figures, 1 table. Submitted to MNRAS, comments welcome. arXiv admin note: text overlap with arXiv:2006.03203

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.

Comments: 17 pages, 11 Figures, accepted for publication in Publications of the Astronomical Society of Australia (18/05/2020). "For the SHAMFI code, see: https://github.com/JLBLine/SHAMFI" . "For the SHAMFI documentation, see: https://shamfi.readthedocs.io/" . "For the WODEN code and documentation see: https://github.com/JLBLine/WODEN"

The bispectrum and 21cm foregrounds during the Epoch of Reionization

Authors: Catherine A Watkinson, Cathryn M Trott, Ian Hothi

Abstract: Numerous studies have established the theoretical potential of the 21cm bispectrum to boost our understanding of the Epoch of Reionization (EoR), and therefore early generation of stars and galaxies. In this paper we look at the impact of foregrounds and instrumental effects on the 21cm bispectrum and our ability to measure it. Unlike the power spectrum for which (in the absence of instrumental ef… ▽ More Numerous studies have established the theoretical potential of the 21cm bispectrum to boost our understanding of the Epoch of Reionization (EoR), and therefore early generation of stars and galaxies. In this paper we look at the impact of foregrounds and instrumental effects on the 21cm bispectrum and our ability to measure it. Unlike the power spectrum for which (in the absence of instrumental effects) there is a window clear of smooth-spectrum foregrounds in which the 21cm power spectrum may be detectable, there is no such "EoR window" for the bispectrum. Instead, on smaller scales EoR structure modulates that of the foregrounds (FG) to alter the EoR+FG bispectrum from that of the foregrounds in a complex manner. On larger scales the EoR structures are completely swamped by those of the foregrounds, and the EoR+FG bispectrum is entirely dominated by that of the foregrounds. It is therefore unlikely that the bispectrum will be useful in cases where we are observationally restricted to using foreground avoidance techniques. We also find that there is potential for instrumental effects to seriously corrupt the bispectrum, possibly even rendering the bispectrum useless for parameter estimation. On larger scales, foreground removal using GMCA is found to recover the EoR bispectrum to a reasonable level of accuracy (even better than the power spectrum for certain configurations). Further studies are necessary to understand the error and/or bias associated with foreground removal before the 21cm bispectrum can be practically applied in analysis of future data. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 14 pages, 13 figures + 2 appendix figures. Submitted to MNRAS (14th February 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.

Comments: 19 pages, 29 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Calibration and 21-cm Power Spectrum Estimation in the Presence of Antenna Beam Variations

Authors: Ronniy C. Joseph, C. M. Trott, R. B. Wayth, A. Nasirudin

Abstract: Detecting a signal from the Epoch of Reionisation (EoR) requires an exquisite understanding of galactic and extra-galactic foregrounds, low frequency radio instruments, instrumental calibration, and data analysis pipelines. In this work we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibra… ▽ More Detecting a signal from the Epoch of Reionisation (EoR) requires an exquisite understanding of galactic and extra-galactic foregrounds, low frequency radio instruments, instrumental calibration, and data analysis pipelines. In this work we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibration errors have the potential to inhibit EoR detections by introducing additional spectral features that mimic the structure of EoR signals. We present a straightforward way to estimate the impact of a wide variety of modelling residuals in EoR PS estimation. We apply this framework to the specific case of broken dipoles in Murchison Widefield Array (MWA) to understand its effect and estimate its impact on PS estimation. Combining an estimate of the percentage of MWA tiles that have at least one broken dipole (15%-40%) with an analytic description of beam errors induced by such dipoles, we compute the residuals of the foregrounds after calibration and source subtraction. We find that that incorrect beam modelling introduces bias in the 2D-PS on the order of $\sim 10^3\, \mathrm{mK}^2 \,h^{-3}\, \mathrm{Mpc}^{3}$. Although this is three orders of magnitude lower than current lowest limits, it is two orders of magnitude higher than the expected signal. Determining the accuracy of both current beam models and direction dependent calibration pipelines is therefore crucial in our search for an EoR signal. △ Less

Submitted 25 February, 2020; v1 submitted 29 November, 2019; originally announced November 2019.

Comments: 14 pages, 8 figures. Accepted for publication in MNRAS December 2019. Erratum submitted on 20 February 2020

First Season MWA Phase II EoR Power Spectrum Results at Redshift 7

Authors: W. Li, J. C. Pober, N. Barry, B. J. Hazelton, M. F. Morales, C. M. Trott, A. Lanman, M. Wilensky, I. Sullivan, A. P. Beardsley, T. Booler, J. D. Bowman, R. Byrne, B. Crosse, D. Emrich, T. M. O. Franzen, K. Hasegawa, L. Horsley, M. Johnston-Hollitt, D. C. Jacobs, C. H. Jordan, R. C. Joseph, T. Kaneuji, D. L. Kaplan, D. Kenney , et al. (22 additional authors not shown)

Abstract: The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of… ▽ More The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $Δ^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$. △ Less

Submitted 20 December, 2019; v1 submitted 22 November, 2019; originally announced November 2019.

Comments: 26 pages, 15 figures, 3 tables. Updated acknowledgements. Published in ApJ

Journal ref: ApJ (2019) 887 141

Science with the Murchison Widefield Array: Phase I Results and Phase II Opportunities

Authors: A. P. Beardsley, M. Johnston-Hollitt, C. M. Trott, J. C. Pober, J. Morgan, D. Oberoi, D. L. Kaplan, C. R. Lynch, G. E. Anderson, P. I. McCauley, S. Croft, C. W. James, O. I. Wong, C. D. Tremblay, R. P. Norris, I. H. Cairns, C. J. Lonsdale, P. J. Hancock, B. M. Gaensler, N. D. R. Bhat, W. Li, N. Hurley-Walker, J. R. Callingham, N. Seymour, S. Yoshiura , et al. (34 additional authors not shown)

Abstract: The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR),… ▽ More The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 38 pages, 12 figures, accepted for publication in PASA

A VOEvent based automatic trigger system for the Murchison Widefield Array

Authors: P. J. Hancock, G. E. Anderson, A. Williams, M. Sokolowski, S. E. Tremblay, A. Rowlinson, B. Crosse, B. W. Meyers, C. R. Lynch, A. Zic, A. P. Beardsley, D. Emrich, T. M. O. Franzen, L. Horsley, M. Johnston-Hollitt, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, K. Steele, S. J. Tingay, C. M. Trott, M. Walker, R. B. Wayth, C. Wu

Abstract: The Murchison Widefield Array (MWA) is an electronically steered low frequency ($<300$\,MHz) radio interferometer, with a `slew' time less than 8seconds. Low frequency ($\sim 100$\,MHz) radio telescopes are ideally suited for rapid-response follow-up of transients due to their large field of view, the inverted spectrum of coherent emission, and the fact that the dispersion delay between a 1GHz and… ▽ More The Murchison Widefield Array (MWA) is an electronically steered low frequency ($<300$\,MHz) radio interferometer, with a `slew' time less than 8seconds. Low frequency ($\sim 100$\,MHz) radio telescopes are ideally suited for rapid-response follow-up of transients due to their large field of view, the inverted spectrum of coherent emission, and the fact that the dispersion delay between a 1GHz and 100MHz pulse is on the order of $1-10$\,min for dispersion measures of $100-2000$\,pc/cm$^3$. The MWA has previously been used to provide fast follow up for transient events including gamma-ray bursts, fast radio bursts, and gravitational waves, using systems that respond to gamma-ray coordinates network (GCN) packet-based notifications. We describe a system for automatically triggering MWA observations of such events, based on VOEvent triggers, which is more flexible, capable, and accurate than previous systems. The system can respond to external multi-messenger triggers, which makes it well-suited to searching for prompt coherent radio emission from gamma-ray bursts, the study of fast radio bursts and gravitational waves, single pulse studies of pulsars, and rapid follow-up of high-energy superflares from flare stars. The new triggering system has the capability to trigger observations in both the regular correlator mode (limited to $\geq 0.5$\,s integrations) or using the Voltage Capture System (VCS, $0.1$\,ms integration) of the MWA, and represents a new mode of operation for the MWA. The upgraded standard correlator triggering capability has been in use since MWA observing semester 2018B (July-Dec 2018), and the VCS and buffered mode triggers will become available for observing in a future semester. △ Less

Submitted 6 October, 2019; originally announced October 2019.

Comments: 14 pages, 5 figures, accepted for publication in PASA

The status of 21cm interferometric experiments

Authors: Cathryn M. Trott, Jonathan C. Pober

Abstract: Interferometric experiments of the reionization era offer the advantages of measuring power in spatial modes with increased sensitivity afforded by multiple independent sky measurements. Here we review early work to measure this signal, current experiments, and future opportunities, highlighting the lessons learned along the way that have shaped the research field and experimental design. In parti… ▽ More Interferometric experiments of the reionization era offer the advantages of measuring power in spatial modes with increased sensitivity afforded by multiple independent sky measurements. Here we review early work to measure this signal, current experiments, and future opportunities, highlighting the lessons learned along the way that have shaped the research field and experimental design. In particular, this chapter discusses the history, progress, challenges and forecasts for detection and exploration of the spatial structure of the 21~cm brightness temperature signal in the Epoch of Reionisation using interferometric experiments. We discuss GMRT, PAPER, LOFAR, MWA, and the future HERA and SKA. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: To appear as a book chapter in "The Cosmic 21-cm Revolution: Charting the first billion years of our Universe", Ed Andrei Mesinger (Bristol: IOP Publishing Ltd) AAS-IOP ebooks http://www.iopscience.org/books/aas

Improving the Epoch of Reionization Power Spectrum Results from Murchison Widefield Array Season 1 Observations

Authors: N. Barry, M. Wilensky, C. M. Trott, B. Pindor, A. P. Beardsley, B. J. Hazelton, I. S. Sullivan, M. F. Morales, J. C. Pober, J. Line, B. Greig, R. Byrne, A. Lanman, W. Li, C. H. Jordan, R. C. Joseph, B. McKinley, M. Rahimi, S. Yoshiura, J. D. Bowman, B. M. Gaensler, J. N. Hewitt, D. C. Jacobs, D. A. Mitchell, N. Udaya Shankar , et al. (5 additional authors not shown)

Abstract: 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 techn… ▽ 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 $Δ^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. △ Less

Submitted 8 October, 2019; v1 submitted 2 September, 2019; originally announced September 2019.

Comments: 21 pages, 10 figures, published in ApJ, updated title and text

A high time resolution study of the millisecond pulsar J2241-5236 at frequencies below 300 MHz

Authors: D. Kaur, N. D. R. Bhat, S. E. Tremblay, R. M. Shannon, S. J. McSweeney, S. M. Ord, A. P. Beardsley, B. Crosse, D. Emrich, T. M. O. Franzen, L. Horsley, M. Johnston-Hollitt, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, K. Steele, S. J. Tingay, C. M. Trott, M. Walker, R. B. Wayth1, A. Williams, C. Wu

Abstract: One of the major challenges for pulsar timing array (PTA) experiments is the mitigation of the effects of the turbulent interstellar medium (ISM) from timing data. These can potentially lead to measurable delays and/or distortions in the pulse profiles and scale strongly with the inverse of the radio frequency. Low-frequency observations are therefore highly appealing for characterizing them. Howe… ▽ More One of the major challenges for pulsar timing array (PTA) experiments is the mitigation of the effects of the turbulent interstellar medium (ISM) from timing data. These can potentially lead to measurable delays and/or distortions in the pulse profiles and scale strongly with the inverse of the radio frequency. Low-frequency observations are therefore highly appealing for characterizing them. However, in order to achieve the necessary time resolution to resolve profile features of short-period millisecond pulsars, phase-coherent de-dispersion is essential, especially at frequencies below $300$ MHz. We present the lowest-frequency ($80$-$220$ MHz), coherently de-dispersed detections of one of the most promising pulsars for current and future PTAs, PSR J2241$-$5236, using our new beam-former software for the MWA's voltage capture system (VCS), which reconstructs the time series at a much higher time resolution of $\sim 1 μ$s by re-synthesizing the recorded voltage data at $10$-kHz/$100$-$μ$s native resolutions. Our data reveal a dual-precursor type feature in the pulse profile that is either faint or absent in high-frequency observations from Parkes. The resultant high-fidelity detections have enabled dispersion measure (DM) determinations with very high precision, of the order of $(2$-$6)\times10^{-6}$ $\rm pc\,cm^{-3}$, owing to the microsecond level timing achievable for this pulsar at the MWA's low frequencies. This underscores the usefulness of low-frequency observations for probing the ISM toward PTA pulsars and informing optimal observing strategies for PTA experiments. △ Less

Submitted 21 July, 2019; originally announced July 2019.

Comments: 09 pages, 4 figures, 2 tables

Gridded and direct Epoch of Reionisation bispectrum estimates using the Murchison Widefield Array

Authors: Cathryn M. Trott, Catherine A. Watkinson, Christopher H. Jordan, Shintaro Yoshiura, Suman Majumdar, N. Barry, R. Byrne, B. J. Hazelton, K. Hasegawa, R. Joseph, T. Kaneuji, K. Kubota, W. Li, J. Line, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales, S. Murray, B. Pindor, J. C. Pober, M. Rahimi, J. Riding, K. Takahashi, S. J. Tingay , et al. (20 additional authors not shown)

Abstract: We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct… ▽ More We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours. △ Less

Submitted 17 May, 2019; originally announced May 2019.

Comments: 19 pages, 10 figures, accepted for publication in PASA

Robust statistics toward detection of the 21 cm signal from the Epoch of Reionisation

Authors: Cathryn M. Trott, Shih Ching Fu, Steven Murray, Christopher Jordan, Jack Line, N. Barry, R. Byrne, B. J. Hazelton, K. Hasegawa, R. Joseph, T. Kaneuji, K. Kubota, W. Li, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales, B. Pindor, J. C. Pober, M. Rahimi, K. Takahashi, S. J. Tingay, R. B. Wayth, R. L. Webster, M. Wilensky , et al. (4 additional authors not shown)

Abstract: We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the t… ▽ More We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Comments: 19 pages, 18 figures, accepted for publication in MNRAS

Hunting for radio emission from the intermittent pulsar J1107-5907 at low frequencies

Authors: B. W. Meyers, S. E. Tremblay, N. D. R. Bhat, C. Flynn, V. Gupta, R. M. Shannon, S. G. Murray, C. Sobey, S. M. Ord, S. Osłowski, B. Crosse, A. Williams, F. Jankowski, W. Farah, V. Venkatraman Krishnan, T. Bateman, M. Bailes, A. Beardsley, D. Emrich, T. M. O. Franzen, B. M. Gaensler, L. Horsley, M. Johnston-Hollitt, D. L. Kaplan, D. Kenney , et al. (8 additional authors not shown)

Abstract: The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molongl… ▽ More The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molonglo Observatory Synthesis Telescope (UTMOST) at 835 MHz, as part of an on-going observing campaign. During a 30-minute simultaneous observation, we detected the pulsar in its bright emission state for approximately 15 minutes, where 86 and 283 pulses were detected above a signal-to-noise threshold of 6 with the MWA and UTMOST, respectively. Of the detected pulses, 51 had counterparts at both frequencies and exhibited steep spectral indices for both the bright main pulse component and the precursor component. We find that the bright state pulse energy distribution is best parameterised by a log-normal distribution at both frequencies, contrary to previous results which suggested a power law distribution. Further low-frequency observations are required in order to explore in detail aspects such as pulse-to-pulse variability, intensity modulations and to better constrain the signal propagation effects due to the interstellar medium and intermittency characteristics at these frequencies. The spectral index, extended profile emission covering a large fraction of pulse longitude, and the broadband intermittency of PSR J1107-5907 suggests that future low-frequency pulsar searches, for instance those planned with SKA-Low, will be in an excellent position to find and investigate new pulsars of this type. △ Less

Submitted 19 December, 2018; v1 submitted 2 November, 2018; originally announced November 2018.

Comments: 19 pages, 7 figures, 4 tables, accepted for publication in ApJ; minor updates to keep in-line with published version

The Bias and Uncertainty of Redundant and Sky-Based Calibration under Realistic Sky and Telescope Conditions

Authors: Ronniy C. Joseph, Cathryn M. Trott, Randall B. Wayth

Abstract: The advent of a new generation of low frequency interferometers has opened a direct window into the Epoch of Reionisation (EoR). However, key to a detection of the faint 21-cm signal, and reaching the sensitivity limits of these arrays, is a detailed understanding of the instruments and their calibration. In this work we use simulations to investigate the bias and uncertainty of redundancy based c… ▽ More The advent of a new generation of low frequency interferometers has opened a direct window into the Epoch of Reionisation (EoR). However, key to a detection of the faint 21-cm signal, and reaching the sensitivity limits of these arrays, is a detailed understanding of the instruments and their calibration. In this work we use simulations to investigate the bias and uncertainty of redundancy based calibration. Specifically, we study the influence of the flux distribution of the radio sky and the impact of antenna position offsets on the complex calibration solutions. We find that the position offsets introduce a bias into the phase component of the calibration solutions. This phase bias increases with the distance between bright radio sources and the pointing center, and with the flux density of these sources. This is potentially problematic for redundant calibration on MWA observations of EoR fields 1 and 2. EoR field 0, however, lacks such sources. We also compared the simulations with theoretical estimates for the bias and uncertainty in sky model based calibration on incomplete sky models for the redundant antenna tiles in the MWA. Our results indicate that redundant calibration outperforms sky based calibration due to the high positional precision of the MWA antenna tiles. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: 15 pages, 10 figures, Accepted for publication in the Astronomical Journal

The effect of baseline layouts on the EoR foreground wedge: a semi-analytical approach

Authors: Steven G. Murray, Cathryn M. Trott

Abstract: The 2D power spectrum is a cornerstone of the modern toolkit for analysis of the low-frequency radio interferometric observations of the 21 cm signal arising from the early Universe. Its familiar form disentangles a great deal of systematic information concerning both the sky and telescope, and displays as a foreground-dominated `brick' and `wedge' on large line-of-sight scales, and a complementar… ▽ More The 2D power spectrum is a cornerstone of the modern toolkit for analysis of the low-frequency radio interferometric observations of the 21 cm signal arising from the early Universe. Its familiar form disentangles a great deal of systematic information concerning both the sky and telescope, and displays as a foreground-dominated `brick' and `wedge' on large line-of-sight scales, and a complementary `window' on smaller scales. This paper builds on many previous works in the literature which seek to elucidate the varied instrumental and foreground factors which contribute to these familiar structures in the 2D power spectrum. In particular, we consider the effects of $uv$-sampling on the emergence of the wedge. Our results verify the expectation that arbitrarily dense instrument layouts in principal restore the missing information that leads to mode-mixing, and can therefore mitigate the wedge feature. We derive rule-of-thumb estimates for the required baseline density for complete wedge mitigation, showing that these will be unachievable in practice. We also discuss the optimal shape of the layout, showing that logarithmic regularity in the radial separation of baselines is favourable. While complete suppression of foreground leakage into the wedge is practically unachievable, we find that designing layouts which promote radial density and regularity is able to reduce the amplitude of foreground power by 1-3 orders of magnitude. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: 21 pages + 4 Appendix Pages, 16 figures, Accepted to ApJ

No low-frequency emission from extremely bright Fast Radio Bursts

Authors: M. Sokolowski, N. D. R. Bhat, J. P. Macquart, R. M. Shannon, K. W. Bannister, R. D. Ekers, D. R. Scott, A. P. Beardsley, B. Crosse, D. Emrich, T. M. O. Franzen, B. M. Gaensler, L. Horsley, M. Johnston-Hollitt, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, G. Sleap, K. Steele, S. J. Tingay, C. M. Trott, M. Walker, R. B. Wayth, A. Williams , et al. (1 additional authors not shown)

Abstract: We present the results of a coordinated campaign conducted with the Murchison Widefield Array (MWA) to shadow Fast Radio Bursts (FRBs) detected by the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz, which resulted in simultaneous MWA observations of seven ASKAP FRBs. We de-dispersed the $24$ $\times$ $1.28$ MHz MWA images across the $170-200$ MHz band taken at 0.5 second time reso… ▽ More We present the results of a coordinated campaign conducted with the Murchison Widefield Array (MWA) to shadow Fast Radio Bursts (FRBs) detected by the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz, which resulted in simultaneous MWA observations of seven ASKAP FRBs. We de-dispersed the $24$ $\times$ $1.28$ MHz MWA images across the $170-200$ MHz band taken at 0.5 second time resolution at the known dispersion measures (DMs) and arrival times of the bursts and searched both within the ASKAP error regions (typically $\sim$ $10$ arcmin $\times$ $10$ arcmin), and beyond ($4$ deg $\times$ $4$ deg). We identified no candidates exceeding a $5σ$ threshold at these DMs in the dynamic spectra. These limits are inconsistent with the mean fluence scaling of $α=-1.8 \pm 0.3$ (${\cal F}_ν\propto ν^α$, where $ν$ is the observing frequency) that is reported for ASKAP events, most notably for the three high fluence (${\cal F}_{1.4\,{\rm GHz}} \gtrsim 100$ Jy ms) FRBs 171020, 180110 and 180324. Our limits show that pulse broadening alone cannot explain our non-detections, and that there must be a spectral turnover at frequencies above 200 MHz. We discuss and constrain parameters of three remaining plausible spectral break mechanisms: free-free absorption, intrinsic spectral turn-over of the radiative processes, and magnification of signals at ASKAP frequencies by caustics or scintillation. If free-free absorption were the cause of the spectral turnover, we constrain the thickness of the absorbing medium in terms of the electron temperature, $T$, to $< 0.03$ $(T/10^4 K)^{-1.35}$ pc for FRB 171020. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Comments: 9 pages, 3 figures , 2 tables (accepted for publication in ApJL)

Journal ref: The Astrophysical Journal Letters, Vol. 867 (2018), Number 1

Assessment of ionospheric activity tolerances for Epoch of Reionisation science with the Murchison Widefield Array

Authors: Cathryn M. Trott, C. H. Jordan, S. G. Murray, B. Pindor, D. A. Mitchell, R. B. Wayth, J. Line, B. McKinley, A. Beardsley, J. Bowman, F. Briggs, B. J. Hazelton, J. Hewitt, D. Jacobs, M. F. Morales, J. C. Pober, S. Sethi, U. Shankar, R. Subrahmanyan, M. Tegmark, S. J. Tingay, R. L. Webster, J. S. B. Wyithe

Abstract: Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky sign… ▽ More Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky signal. We present a general framework for understanding the signatures of ionospheric behaviour in the two-dimensional (2D) neutral hydrogen power spectrum measured by a low-frequency radio interferometer. Two primary classes of ionospheric behaviour are considered, corresponding to dominant modes observed in Murchison Widefield Array (MWA) EoR data; namely, anisotropic structured wave behaviour, and isotropic turbulence. Analytic predictions for power spectrum bias due to this contamination are computed, and compared with simulations. We then apply the ionospheric metric described in Jordan et al. (2017) to study the impact of ionospheric structure on MWA data, by dividing MWA EoR datasets into classes with good and poor ionospheric conditions, using sets of matched 30-minute observations from 2014 September. The results are compared with the analytic and simulated predictions, demonstrating the observed bias in the power spectrum when the ionosphere is active (displays coherent structures or isotropic turbulence). The analysis demonstrates that unless ionospheric activity can be quantified and corrected, active data should not be included in EoR analysis in order to avoid systematic biases in cosmological power spectra. When data are corrected with a model formed from the calibration information, bias reduces below the expected 21~cm signal level. Data are considered `quiet' when the median measured source position offsets are less than 10-15~arcseconds. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: 19 pages, 11 figures, Accepted for publication in the Astrophysical Journal

The Phase II Murchison Widefield Array: Design Overview

Authors: Randall B. Wayth, Steven J. Tingay, Cathryn M. Trott, David Emrich, Melanie Johnston-Hollitt, Ben McKinley, B. M. Gaensler, A. P. Beardsley, T. Booler, B. Crosse, T. M. O. Franzen, L. Horsley, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, G. Sleap, K. Steele, M. Walker, A. Williams, C. Wu, Iver. H. Cairns, M. D. Filipovic, S. Johnston, T. Murphy , et al. (4 additional authors not shown)

Abstract: We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These n… ▽ More We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These new tiles enhance the surface brightness sensitivity of the MWA and will improve the ability of the MWA to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ~3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u,v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of MWA continuum images. The upgrade retains all of the features that have underpinned the MWA's success (large field-of-view, snapshot image quality, pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: 11 pages, 8 figs. Accepted for publication in PASA

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.

Comments: 13 pages, 7 figures, accepted in MNRAS

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.

Comments: 13 pages, 7 figures Accepted to PASA

Journal ref: Publ. Astron. Soc. Aust. 35 (2018) e045

Observations of Low-Frequency Radio Emission from Millisecond Pulsars and Multipath Propagation in the Interstellar Medium

Authors: N. D. R. Bhat, S. E. Tremblay, F. Kirsten, B. W. Meyers, M. Sokolowski, W. van Straten, S. J. McSweeney, S. M. Ord, R. M. Shannon, A. Beardsley, B. Crosse, D. Emrich, T. M. O. Franzen, L. Horsley, M. Johnston-Hollitt, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, K. Steele, S. J. Tingay, C. M. Trott, M. Walker, R. B. Wayth, A. Williams , et al. (1 additional authors not shown)

Abstract: Studying the gravitational-wave sky with pulsar timing arrays (PTAs) is a key science goal for the Square Kilometre Array (SKA) and its pathfinder telescopes. With current PTAs reaching sub-microsecond timing precision, making accurate measurements of interstellar propagation effects and mitigating them effectively has become increasingly important to realise PTA goals. As these effects are much s… ▽ More Studying the gravitational-wave sky with pulsar timing arrays (PTAs) is a key science goal for the Square Kilometre Array (SKA) and its pathfinder telescopes. With current PTAs reaching sub-microsecond timing precision, making accurate measurements of interstellar propagation effects and mitigating them effectively has become increasingly important to realise PTA goals. As these effects are much stronger at longer wavelengths, low-frequency observations are most appealing for characterizing the interstellar medium (ISM) along the sight lines toward PTA pulsars. The Murchison Widefield Array (MWA) and the Engineering Development Array (EDA), which utilizes MWA technologies, present promising opportunities for undertaking such studies, particularly for PTA pulsars located in the southern sky. Such pulsars are also the prime targets for PTA efforts planned with the South African MeerKAT, and eventually with the SKA. In this paper we report on observations of two bright southern millisecond pulsars PSRs J0437-4715 and J2145-0750 made with these facilities; MWA observations sampling multiple frequencies across the 80-250 MHz frequency range, while the EDA providing direct-sampled baseband data to yield a large instantaneous usable bandwidth of $\sim$200 MHz. Using these observations, we investigate various aspects relating to pulsar emission and ISM properties, such as spectral evolution of the mean pulse shape, scintillation as a function of frequency, chromaticity in interstellar dispersion, and flux density spectra at low frequencies. Systematic and regular monitoring observations will help ascertain the role of low-frequency measurements in PTA experiments, while simultaneously providing a detailed characterization of the ISM toward the pulsars, which will be useful in devising optimal observing strategies for future PTA experiments. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 23 pages, 10 Figures, 4 Tables, Accepted to ApJS

Survey parameters for detecting 21cm - Ly$α$ emitter cross correlations with the Square Kilometre Array

Authors: Anne Hutter, Cathryn M. Trott, Pratika Dayal

Abstract: Detections of the cross correlation signal between the 21cm signal during reionization and high-redshift Lyman Alpha emitters (LAEs) are subject to observational uncertainties which mainly include systematics associated with radio interferometers and LAE selection. These uncertainties can be reduced by increasing the survey volume and/or the survey luminosity limit, i.e. the faintest detectable Ly… ▽ More Detections of the cross correlation signal between the 21cm signal during reionization and high-redshift Lyman Alpha emitters (LAEs) are subject to observational uncertainties which mainly include systematics associated with radio interferometers and LAE selection. These uncertainties can be reduced by increasing the survey volume and/or the survey luminosity limit, i.e. the faintest detectable Lyman Alpha (Ly$α$) luminosity. We use our model of high-redshift LAEs and the underlying reionization state to compute the uncertainties of the 21cm-LAE cross correlation function at $z\simeq6.6$ for observations with SKA1-Low and LAE surveys with $Δz=0.1$ for three different values of the average IGM ionization state ($\langleχ_\mathrm{HI}\rangle$=0.1, 0.25, 0.5). At $z\simeq6.6$, we find SILVERRUSH type surveys, with a field of view of 21 deg$^2$ and survey luminosity limits of $L_α\geq7.9\times10^{42}$erg~s$^{-1}$, to be optimal to distinguish between an inter-galactic medium (IGM) that is 50%, 25% and 10% neutral, while surveys with smaller fields of view and lower survey luminosity limits, such as the 5 and 10 deg$^2$ surveys with WFIRST, can only discriminate between a 50% and 10% neutral IGM. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS

Limits on radio emission from meteors using the MWA

Authors: Xiang Zhang, Paul Hancock, Hadrien A. R. Devillepoix, Randall B. Wayth, A. Beardsley, B. Crosse, D. Emrich, T. M. O. Franzen, B. M. Gaensler, L. Horsley, 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, Jianghui Ji, Yuehua Ma

Abstract: Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_ν\propto ν^α$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteo… ▽ More Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_ν\propto ν^α$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteor emission and other transient events using the Murchison Widefield Array (MWA) at 72-103 MHz. In our 322-hour survey, down to a $5σ$ detection threshold of 3.5 Jy/beam, no transient candidates were identified as intrinsic emission from meteors. We derived an upper limit of -3.7 (95% confidence limit) on the spectral index in our frequency range. We also report detections of other transient events, like reflected FM broadcast signals from small satellites, conclusively demonstrating the ability of the MWA to detect and track space debris on scales as small as 0.1 m in low Earth orbits. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: Accepted for Publication in MNRAS on 11 April 2018, 11 pages, 11 figures, 5 tables

Study of systematics effects on the Cross Power Spectrum of 21 cm Line and Cosmic Microwave Background using Murchison Widefield Array Data

Authors: S. Yoshiura, K. Ichiki, B. Pindor, K. Takahashi, H. Tashiro, C. M. Trott

Abstract: Observation of the 21cm line signal from neutral hydrogen during the Epoch of Reionization is challenging due to extremely bright Galactic and extragalactic foregrounds and complicated instrumental calibration. A reasonable approach for mitigating these problems is the cross correlation with other observables. In this work, we present the first results of the cross power spectrum (CPS) between rad… ▽ More Observation of the 21cm line signal from neutral hydrogen during the Epoch of Reionization is challenging due to extremely bright Galactic and extragalactic foregrounds and complicated instrumental calibration. A reasonable approach for mitigating these problems is the cross correlation with other observables. In this work, we present the first results of the cross power spectrum (CPS) between radio images observed by the Murchison Widefield Array and the cosmic microwave background (CMB), measured by the Planck experiment. We study the systematics due to the ionospheric activity, the dependence of CPS on group of pointings, and frequency. The resulting CPS is consistent with zero because the error is dominated by the foregrounds in the 21cm observation. Additionally, the variance of the signal indicates the presence of unexpected systematics error at small scales. Furthermore, we reduce the error by one order of magnitude with application of a foreground removal using a polynomial fitting method. Based on the results, we find that the detection of the 21cm-CMB CPS with the MWA Phase I requires more than 99.95% of the foreground signal removed, 2000 hours of deep observation and 50% of the sky fraction coverage. △ Less

Submitted 27 November, 2018; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 15 pages, 16 figures, accepted to MNRAS