Showing 1–50 of 162 results for author: Webster, R L

The accretion of a solar mass per day by a 17-billion solar mass black hole

Authors: Christian Wolf, Samuel Lai, Christopher A. Onken, Neelesh Amrutha, Fuyan Bian, Wei Jeat Hon, Patrick Tisserand, Rachel L. Webster

Abstract: Around a million quasars have been catalogued in the Universe by probing deeper and using new methods for discovery. However, the hardest ones to find seem to be the rarest and brightest specimen. In this work, we study the properties of the most luminous of all quasars found so far. It has been overlooked until recently, which demonstrates that modern all-sky surveys have much to reveal. The blac… ▽ More Around a million quasars have been catalogued in the Universe by probing deeper and using new methods for discovery. However, the hardest ones to find seem to be the rarest and brightest specimen. In this work, we study the properties of the most luminous of all quasars found so far. It has been overlooked until recently, which demonstrates that modern all-sky surveys have much to reveal. The black hole in this quasar accretes around one solar mass per day onto an existing mass of $\sim$17 billion solar masses. In this process its accretion disc alone releases a radiative energy of $2\times 10^{41}$ Watts. If the quasar is not strongly gravitationally lensed, then its broad line region (BLR) is expected to have the largest physical and angular diameter occurring in the Universe, and will allow the Very Large Telescope Interferometer to image its rotation and measure its black hole mass directly. This will be an important test for BLR size-luminosity relations, whose extrapolation has underpinned common black-hole mass estimates at high redshift. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: This version of the article has been accepted for publication after peer review and is subject to Springer Nature's AM terms of use, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at http://dx.doi.org/10.1038/s41550-024-02195-x

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 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

Dual Polarization Measurements of MWA Beampatterns at 137 MHz

Authors: A. Chokshi, J. L. B. Line, N. Barry, D. Ung, D. Kenney, A. McPhail, A. Williams, R. L. Webster

Abstract: The wide adoption of low-frequency radio interferometers as a tool for deeper and higher resolution astronomical observations has revolutionised radio astronomy. Despite their construction from static, relatively simple dipoles, the sheer number of distinct elements introduces new, complicated instrumental effects. Their necessary remote locations exacerbate failure rates, while electronic interac… ▽ More The wide adoption of low-frequency radio interferometers as a tool for deeper and higher resolution astronomical observations has revolutionised radio astronomy. Despite their construction from static, relatively simple dipoles, the sheer number of distinct elements introduces new, complicated instrumental effects. Their necessary remote locations exacerbate failure rates, while electronic interactions between the many adjacent receiving elements can lead to non-trivial instrumental effects. The Murchison Widefield Array (MWA) employs phased array antenna elements (tiles), which improve collecting area at the expense of complex beam shapes. Advanced electromagnetic simulations have produced the Fully Embedded Element (FEE) simulated beam model which has been highly successful in describing the ideal beam response of MWA antennas. This work focuses on the relatively unexplored aspect of various in-situ, environmental perturbations to beam models and represents the first large-scale, in-situ, all-sky measurement of MWA beam shapes at multiple polarizations and pointings. Our satellite based beam measurement approach enables all-sky beam response measurements with a dynamic range of \sim 50 dB, at 137 MHz. △ Less

Submitted 31 January, 2021; originally announced February 2021.

Comments: 16 pages, 15 figures, accepted for publication by the Monthly Notices of the Royal Astronomical Society (15/01/2021). For EMBERS code see: https://github.com/amanchokshi/EMBERS . For EMBERS documentation see: https://embers.readthedocs.io

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

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"

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

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

Determining Quasar Orientation

Authors: Suk Yee Yong, Rachel L. Webster, Anthea L. King, Nicholas F. Bate, Kathleen Labrie, Matthew J. O'Dowd

Abstract: Since the discovery of active galactic nuclei (AGN) and their subclasses, a unification scheme of AGN has been long sought. Orientation-based unified models predict that some of the diversity within AGN subclasses can be explained by the different viewing angles of the observer. Several orientation categorisations have been suggested, but a widely applicable measure has yet to be found. Using the… ▽ More Since the discovery of active galactic nuclei (AGN) and their subclasses, a unification scheme of AGN has been long sought. Orientation-based unified models predict that some of the diversity within AGN subclasses can be explained by the different viewing angles of the observer. Several orientation categorisations have been suggested, but a widely applicable measure has yet to be found. Using the properties of the ultraviolet and optical broad emission lines of quasars, in particular the velocity offsets and line widths of high-ionisation CIV and low-ionisation MgII lines, a correlation has been measured. It is postulated that this correlation is due to the viewing angle of the observer. Comparison with other orientation tracers shows consistency with this interpretation. Using a simulation of a wide angle disk-wind model for the broad emission line region, we successfully replicate the observed correlation with inclination. Future more detailed modelling will not only enable improved accuracy in the determination of the viewing angle to individual AGN, but will also substantially increase our understanding of the emitting regions of AGN. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS. 15 pages, 17 figures, 2 table

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

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

The Murchison Widefield Array Transients Survey (MWATS). A search for low frequency variability in a bright Southern hemisphere sample

Authors: M. E. Bell, Tara Murphy, P. J. Hancock, J. R. Callingham, S. Johnston, D. L. Kaplan, R. W. Hunstead, E. M. Sadler, S. Croft, S. V. White, N. Hurley-Walker, R. Chhetri, J. S. Morgan, P. G. Edwards, A. Rowlinson, A. R. Offringa, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt , et al. (16 additional authors not shown)

Abstract: We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>… ▽ More We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 19 pages, 11 figures

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

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

Comparing Redundant and Sky Model Based Interferometric Calibration: A First Look with Phase II of the MWA

Authors: W. Li, J. C. Pober, B. J. Hazelton, N. Barry, M. F. Morales, I. Sullivan, A. R. Parsons, Z. S. Ali, J. S. Dillon, A. P. Beardsley, J. D. Bowman, F. Briggs, R. Byrne, P. Carroll, B. Crosse, D. Emrich, A. Ewall-Wice, L. Feng, T. M. O. Franzen, J. N. Hewitt, L. Horsley, D. C. Jacobs, M. Johnston-Hollitt, C. Jordan, R. C. Joseph , et al. (31 additional authors not shown)

Abstract: Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly u… ▽ More Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: 20 pages, 11 figures. Accepted to ApJ

HST imaging of four gravitationally lensed quasars

Authors: Nicholas F. Bate, Georgios Vernardos, Matthew J. O'Dowd, Daniel M. Neri-Larios, Rachel L. Webster, David J. E. Floyd, Robert L. Barone-Nugent, Kathleen Labrie, Anthea L. King, Suk Yee Yong

Abstract: We present new HST WFPC3 imaging of four gravitationally lensed quasars: MG 0414+0534; RXJ 0911+0551; B 1422+231; WFI J2026-4536. In three of these systems we detect wavelength-dependent microlensing, which we use to place constraints on the sizes and temperature profiles of the accretion discs in each quasar. Accretion disc radius is assumed to vary with wavelength according to the power-law rela… ▽ More We present new HST WFPC3 imaging of four gravitationally lensed quasars: MG 0414+0534; RXJ 0911+0551; B 1422+231; WFI J2026-4536. In three of these systems we detect wavelength-dependent microlensing, which we use to place constraints on the sizes and temperature profiles of the accretion discs in each quasar. Accretion disc radius is assumed to vary with wavelength according to the power-law relationship $r\propto λ^p$, equivalent to a radial temperature profile of $T\propto r^{-1/p}$. The goal of this work is to search for deviations from standard thin disc theory, which predicts that radius goes as wavelength to the power $p=4/3$. We find a wide range of power-law indices, from $p=1.4^{+0.5}_{-0.4}$ in B 1422+231 to $p=2.3^{+0.5}_{-0.4}$ in WFI J2026-4536. The measured value of $p$ appears to correlate with the strength of the wavelength-dependent microlensing. We explore this issue with mock simulations using a fixed accretion disc with $p=1.5$, and find that cases where wavelength-dependent microlensing is small tend to under-estimate the value of $p$. This casts doubt on previous ensemble single-epoch measurements which have favoured low values using samples of lensed quasars that display only moderate chromatic effects. Using only our systems with strong chromatic microlensing we prefer $p>4/3$, corresponding to shallower temperature profiles than expected from standard thin disc theory. △ Less

Submitted 10 July, 2018; originally announced July 2018.

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

Using the Properties of Broad Absorption Line Quasars to Illuminate Quasar Structure

Authors: Suk Yee Yong, Anthea L. King, Rachel L. Webster, Nicholas F. Bate, Matthew J. O'Dowd, Kathleen Labrie

Abstract: A key to understanding quasar unification paradigms is the emission properties of broad absorption line quasars (BALQs). The fact that only a small fraction of quasar spectra exhibit deep absorption troughs blueward of the broad permitted emission lines provides a crucial clue to the structure of quasar emitting regions. To learn whether it is possible to discriminate between the BALQ and non-BALQ… ▽ More A key to understanding quasar unification paradigms is the emission properties of broad absorption line quasars (BALQs). The fact that only a small fraction of quasar spectra exhibit deep absorption troughs blueward of the broad permitted emission lines provides a crucial clue to the structure of quasar emitting regions. To learn whether it is possible to discriminate between the BALQ and non-BALQ populations given the observed spectral properties of a quasar, we employ two approaches: one based on statistical methods and the other supervised machine learning classification, applied to quasar samples from the Sloan Digital Sky Survey. The features explored include continuum and emission line properties, in particular the absolute magnitude, redshift, spectral index, line width, asymmetry, strength, and relative velocity offsets of high-ionisation CIV $λ1549$ and low-ionisation MgII $λ2798$ lines. We consider a complete population of quasars, and assume that the statistical distributions of properties represent all angles where the quasar is viewed without obscuration. The distributions of the BALQ and non-BALQ sample properties show few significant differences. None of the observed continuum and emission line features are capable of differentiating between the two samples. Most published narrow disk-wind models are inconsistent with these observations, and an alternative disk-wind model is proposed. The key feature of the proposed model is a disk-wind filling a wide opening angle with multiple radial streams of dense clumps. △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: Accepted for publication in MNRAS. 19 pages, 3 figures, 6 table

Near-identical star formation rate densities from H$α$ and FUV at redshift zero

Authors: Fiona M. Audcent-Ross, Gerhardt R. Meurer, O. I. Wong, Z. Zheng, D. Hanish, M. A. Zwaan, J. Bland-Hawthorn, A. Elagali, M. Meyer, M. E. Putman, E. V. Ryan-Webber, S. M. Sweet, D. A. Thilker, M. Seibert, R. Allen, M. A. Dopita, M. T. Doyle-Pegg, M. Drinkwater, H. C. Ferguson, K. C. Freeman, T. M. Heckman, R. C. Kennicutt Jr, V. A. Kilborn, J. H. Kim, P. M. Knezek , et al. (5 additional authors not shown)

Abstract: For the first time both H$α$ and far-ultraviolet (FUV) observations from an HI-selected sample are used to determine the dust-corrected star formation rate density (SFRD: $\dotρ$) in the local Universe. Applying the two star formation rate indicators on 294 local galaxies we determine log($\dotρ$$ _{Hα}) = -1.68~^{+0.13}_{-0.05}$ [M$_{\odot} $ yr$^{-1} $ Mpc$^{-3}]$ and log($\dotρ_{FUV}$)… ▽ More For the first time both H$α$ and far-ultraviolet (FUV) observations from an HI-selected sample are used to determine the dust-corrected star formation rate density (SFRD: $\dotρ$) in the local Universe. Applying the two star formation rate indicators on 294 local galaxies we determine log($\dotρ$$ _{Hα}) = -1.68~^{+0.13}_{-0.05}$ [M$_{\odot} $ yr$^{-1} $ Mpc$^{-3}]$ and log($\dotρ_{FUV}$) $ = -1.71~^{+0.12}_{-0.13}$ [M$_\odot $ yr$^{-1} $ Mpc$^{-3}]$. These values are derived from scaling H$α$ and FUV observations to the HI mass function. Galaxies were selected to uniformly sample the full HI mass (M$_{HI}$) range of the HI Parkes All-Sky Survey (M$_{HI} \sim10^{7}$ to $\sim10^{10.7}$ M$_{\odot}$). The approach leads to relatively larger sampling of dwarf galaxies compared to optically-selected surveys. The low HI mass, low luminosity and low surface brightness galaxy populations have, on average, lower H$α$/FUV flux ratios than the remaining galaxy populations, consistent with the earlier results of Meurer. The near-identical H$α$- and FUV-derived SFRD values arise with the low H$α$/FUV flux ratios of some galaxies being offset by enhanced H$α$ from the brightest and high mass galaxy populations. Our findings confirm the necessity to fully sample the HI mass range for a complete census of local star formation to include lower stellar mass galaxies which dominate the local Universe. △ Less

Submitted 15 June, 2018; originally announced June 2018.

Comments: 17 pages, 7 figures

The Intrinsic Far-UV Spectrum of B1422+231

Authors: M. O'Dowd, N. F. Bate, R. L. Webster, K. Labrie, A. L. King, S-. Y. Yong

Abstract: We present new spectroscopy of the z=3.62 gravitationally lensed quasar B1422+117 from the Gemini North GMOS integral field spectrograph. We observe significant differential magnifications between the broad emission lines and the continuum, as well as across the velocity structure of the Lyman-alpha line. We take advantage of this differential microlensing to algebraically decompose the quasar spe… ▽ More We present new spectroscopy of the z=3.62 gravitationally lensed quasar B1422+117 from the Gemini North GMOS integral field spectrograph. We observe significant differential magnifications between the broad emission lines and the continuum, as well as across the velocity structure of the Lyman-alpha line. We take advantage of this differential microlensing to algebraically decompose the quasar spectrum into the absorbed broad emission line and absorbed continuum components. We use the latter to derive the intrinsic Ly-alpha forest absorption spectrum. The proximity effect is clearly detected, with a proximity zone edge of 8.6-17.3 Mpc from the quasar, implying (perhaps intermittent) activity over at least 28 Myrs. The Ly-alpha line profile exhibits a blue excess that is inconsistent with a symmetric fit to the unabsorbed red side. This has important implications for the use of this fitting technique in estimating the absorbed blue Ly-alpha wings of Gunn-Peterson trough quasars. △ Less

Submitted 4 October, 2017; v1 submitted 24 August, 2017; originally announced August 2017.

Comments: Accepted for publication in MNRAS, 10 pages, 9 figures

The Kinematics of Quasar Broad Emission Line Regions Using a Disk-wind Model

Authors: Suk Yee Yong, Rachel L. Webster, Anthea L. King, Nicholas F. Bate, Matthew J. O'Dowd, Kathleen Labrie

Abstract: The structure and kinematics of the broad line region (BLR) in quasars are still not well established. One popular BLR model is the disk-wind model that offers a geometric unification of a quasar based on the angle of viewing. We construct a simple kinematical disk-wind model with a narrow outflowing wind angle. The model is combined with radiative transfer in the Sobolev, or high velocity, limit.… ▽ More The structure and kinematics of the broad line region (BLR) in quasars are still not well established. One popular BLR model is the disk-wind model that offers a geometric unification of a quasar based on the angle of viewing. We construct a simple kinematical disk-wind model with a narrow outflowing wind angle. The model is combined with radiative transfer in the Sobolev, or high velocity, limit. We examine how angle of viewing affects the observed characteristics of the emission line, especially the line widths and velocity offsets. The line profiles exhibit distinct properties depending on the orientation, wind opening angle, and region of the wind where the emission arises. At low inclination angle (close to face-on), we find the shape of the emission line is asymmetric with narrow width and significantly blueshifted. As the inclination angle increases (close to edge-on), the line profile becomes more symmetric, broader, and less blueshifted. Additionally, lines that arise close to the base of the disk wind, near the accretion disk, tend to be broad and symmetric. The relative increase in blueshift of the emission line with increasing wind vertical distance is larger for polar winds compared with equatorial winds. By considering the optical thickness of the wind, single-peaked line profiles are recovered for the intermediate and equatorial outflowing wind. The model is able to reproduce a faster response in either the red and blue sides of the line profile found in reverberation mapping studies. A quicker response in the red side is achieved in the model with a polar wind and intermediate wind opening angle at low viewing angle. The blue side response is faster for an equatorial wind seen at high inclination. △ Less

Submitted 23 August, 2017; originally announced August 2017.

Comments: Accepted for publication in PASA. 15 pages, 6 figures, 1 table

Detectability of 21cm-signal during the Epoch of Reionization with 21cm-Lyman-$α$ emitter cross-correlation. I

Authors: Kenji Kubota, Shintaro Yoshiura, Keitaro Takahashi, Kenji Hasegawa, Hidenobu Yajima, Masami Ouchi, B. Pindor, R. L. Webster

Abstract: Detection of the redshifted 21cm-line signal from neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) is complicated by intense foregrounds such as galactic synchrotron and extragalactic radio galaxies. The 21cm-Lyman-$α$ emitter(LAE) cross-correlation is one of the tools available to reduce the foreground effects because the foreground emission from such radi… ▽ More Detection of the redshifted 21cm-line signal from neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) is complicated by intense foregrounds such as galactic synchrotron and extragalactic radio galaxies. The 21cm-Lyman-$α$ emitter(LAE) cross-correlation is one of the tools available to reduce the foreground effects because the foreground emission from such radio sources is statistically independent of LAE distribution. LAE surveys during the EoR at redshifts $z=6.6$ and $7.3$ are ongoing by the Subaru Hyper Suprime-Cam (HSC). Additionally, Prime Focus Spectrograph (PFS) will provide precise redshift information of the LAEs discovered by the HSC survey. In this paper, we investigate the detectability of the 21cm signal with the 21cm-LAE cross-correlation by using our improved reionization simulations. We also focus on the error budget and evaluate it quantitatively in order to consider a strategy to improve the signal-to-noise ratio. In addition, we explore an expansion of the LAE survey to suggest optimal survey parameters and show a potential to measure a characteristic size of ionized bubbles via the turnover scale of the cross-power spectrum. As a result, we find that the Murchison Widefield Array (MWA) has ability to detect the cross-power spectrum signal on large scales by combining LAE Deep field survey of HSC. We also show that the sensitivity is improved dramatically at small scales by adding redshift information from the PFS measurements. The Square Kilometre Array (SKA) has a potential to measure the turnover scale with an accuracy of $6\times10^{-3}~{\rm Mpc^{-1}}$. △ Less

Submitted 21 June, 2018; v1 submitted 21 August, 2017; originally announced August 2017.

Comments: 13pages, 13figures, 3tables, accepted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 479, Issue 2, 11 September 2018, Pages 2754

Spectral energy distribution and radio halo of NGC 253 at low radio frequencies

Authors: A. D. Kapinska, L. Staveley-Smith, R. Crocker, G. R. Meurer, S. Bhandari, N. Hurley-Walker, A. R. Offringa, D. J. Hanish, N. Seymour, R. D. Ekers, M. E. Bell, J. R. Callingham, K. S. Dwarakanath, B. -Q. For, B. M. Gaensler, P. J. Hancock, L. Hindson, M. Johnston-Hollitt, E. Lenc, B. McKinley, J. Morgan, P. Procopio, R. B. Wayth, C. Wu, Q. Zheng , et al. (45 additional authors not shown)

Abstract: We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the sta… ▽ More We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis). △ Less

Submitted 19 February, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: Accepted for publication in Astrophysical Journal on 06 February 2017

A Matched Filter Technique For Slow Radio Transient Detection And First Demonstration With The Murchison Widefield Array

Authors: L. Feng, R. Vaulin, J. N. Hewitt, R. Remillard, D. L. Kaplan, Tara Murphy, N. Kudryavtseva, P. Hancock, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, D. Oberoi, S. M. Ord, T. Prabu , et al. (8 additional authors not shown)

Abstract: Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble… ▽ More Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between ~20--200 mJy, providing the best limits to date for hour- and month-long transients. △ Less

Submitted 12 January, 2017; originally announced January 2017.

Comments: 16 pages, 14 figures, accepted to AJ

Wavelet-Based Characterization of Small-Scale Solar Emission Features at Low Radio Frequencies

Authors: Akshay Suresh, Rohit Sharma, Divya Oberoi, Srijan B. Das, Victor Pankratius, Brian Timar, Colin J. Lonsdale, Judd D. Bowman, Frank Briggs, Roger J. Cappallo, Brian E. Corey, Avinash A. Deshpande, David Emrich, Robert Goeke, Lincoln J. Greenhill, Bryna J. Hazelton, Melanie Johnston-Hollitt, David L. Kaplan, Justin C. Kasper, Eric Kratzenberg, Mervyn J. Lynch, S. Russell McWhirter, Daniel A. Mitchell, Miguel F. Morales, Edward Morgan , et al. (14 additional authors not shown)

Abstract: Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detecti… ▽ More Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts. △ Less

Submitted 2 June, 2017; v1 submitted 3 December, 2016; originally announced December 2016.

Comments: 15 pages, 14 figures, Accepted for publication in ApJ

PUMA: The Positional Update and Matching Algorithm

Authors: J. L. B. Line, R. L. Webster, B. Pindor, D. A. Mitchell, C. M. Trott

Abstract: We present new software to cross-match low-frequency radio catalogues: the Positional Update and Matching Algorithm (PUMA). PUMA combines a positional Bayesian probabilistic approach with spectral matching criteria, allowing for confusing sources in the matching process. We go on to create a radio sky model using PUMA based on the Murchison Widefield Array Commissioning Survey, and are able to aut… ▽ More We present new software to cross-match low-frequency radio catalogues: the Positional Update and Matching Algorithm (PUMA). PUMA combines a positional Bayesian probabilistic approach with spectral matching criteria, allowing for confusing sources in the matching process. We go on to create a radio sky model using PUMA based on the Murchison Widefield Array Commissioning Survey, and are able to automatically cross-match 98.5% of sources. Using the characteristics of this sky model, we create simple simulated mock catalogues on which to test PUMA, and find that PUMA can reliably find the correct spectral indices of sources, along with being able to recover ionospheric offsets. Finally, we use this sky model to calibrate and remove foreground sources from simulated interferometric data, generated using OSKAR (the Oxford University visibility generator). We demonstrate that there is a substantial improvement in foreground source removal when using higher frequency and higher resolution source positions, even when correcting positions by an average of 0.3 given a synthesized beam-width of 2.3. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: Accepted for publication in Publications of the Astronomical Society of Australia (PASA) - 34 pages - 20 figures - PUMA github repository: https://github.com/JLBLine/PUMA

GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue

Authors: Natasha Hurley-Walker, Joseph R. Callingham, Paul J. Hancock, Thomas M. O. Franzen, Luke Hindson, Anna D. Kapinska, John Morgan, Andre R. Offringa, Randall B. Wayth, Chen Wu, Q. Zheng, Tara Murphy, Martin E. Bell, K. S. Dwarakanath, Bi-Qing For, Bryan M. Gaensler, Melanie Johnston-Hollitt, Emil Lenc, Pietro Procopio, Lister Staveley-Smith, Ron Ekers, Judd D. Bowman, Frank Briggs, R. J. Cappallo, Avinash A. Deshpande , et al. (18 additional authors not shown)

Abstract: Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic… ▽ More Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307,455 radio sources with 20 separate flux density measurements across 72--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5σ$, which itself is typically 50mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low frequency flux density scale of the southern sky to better than 10%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: 30 pages, 18 figures, 6 tables, published in Monthly Notices of the Royal Astronomical Society

Delay Spectrum with Phase-Tracking Arrays: Extracting the HI power spectrum from the Epoch of Reionization

Authors: Sourabh Paul, Shiv K. Sethi, Miguel F. Morales, K. S. Dwarkanath, N. Udaya Shankar, Ravi Subrahmanyan, N. Barry, A. P. Beardsley, Judd D. Bowman, F. Briggs, P. Carroll, A. de Oliveira-Costa, Joshua S. Dillon, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. M. Gaensler, B. J. Hazelton, J. N. Hewitt, N. Hurley-Walker, D. J. Jacobs, Han-Seek Kim, P. Kittiwisit, E. Lenc, J. Line , et al. (29 additional authors not shown)

Abstract: The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our met… ▽ More The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data. △ Less

Submitted 22 October, 2016; originally announced October 2016.

Comments: 17 pages, 10 figures, accepted for publication in The Astrophysical Journal

Journal ref: ApJ, 2016, 833, 2, 213

First Season MWA EoR Power Spectrum Results at Redshift 7

Authors: A. P. Beardsley, B. J. Hazelton, I. S. Sullivan, P. Carroll, N. Barry, M. Rahimi, B. Pindor, C. M. Trott, J. Line, Daniel C. Jacobs, M. F. Morales, J. C. Pober, G. Bernardi, Judd D. Bowman, M. P. Busch, F. Briggs, R. J. Cappallo, B. E. Corey, A. de Oliveira-Costa, Joshua S. Dillon, D. Emrich, A. Ewall-Wice, L. Feng, B. M. Gaensler, R. Goeke , et al. (41 additional authors not shown)

Abstract: The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and qu… ▽ More The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of $Δ^2 \leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: 21 pages, 16 figures, accepted to ApJ

Low frequency observations of linearly polarized structures in the interstellar medium near the south Galactic pole

Authors: Emil Lenc, B. M. Gaensler, X. H. Sun, E. M. Sadler, A. G. Willis, N. Barry, A. P. Beardsley, M. E. Bell, G. Bernardi, J. D. Bowman, F. Briggs, J. R. Callingham, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, A. A. Deshpande, J. S. Dillon, K. S. Dwarkanath, D. Emrich, A. Ewall-Wice, L. Feng, B. -Q. For, R. Goeke, L. J. Greenhill , et al. (54 additional authors not shown)

Abstract: We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but… ▽ More We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ~54' resolution with a sensitivity of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz. △ Less

Submitted 3 August, 2016; v1 submitted 19 July, 2016; originally announced July 2016.

Comments: 32 pages, 20 figures, accepted to ApJ

A High Reliability Survey of Discrete Epoch of Reionization Foreground Sources in the MWA EoR0 Field

Authors: P. A. Carroll, J. Line, M. F. Morales, N. Barry, A. P. Beardsley, B. J. Hazelton, D. C. Jacobs, J. C. Pober, I. S. Sullivan, R. L. Webster, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. de Oliveira-Costa, J. S. Dillon, D. Emrich, A. Ewall-Wice, L. Feng, B. M. Gaensler, R. Goeke, L. J. Greenhill, J. N. Hewitt, N. Hurley-Walker , et al. (39 additional authors not shown)

Abstract: Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision… ▽ More Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction. △ Less

Submitted 13 July, 2016; originally announced July 2016.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal on June 30, 2016

The radio spectral energy distribution of infrared-faint radio sources

Authors: A. Herzog, R. P. Norris, E. Middelberg, N. Seymour, L. R. Spitler, B. H. C. Emonts, T. M. O. Franzen, R. Hunstead, H. T. Intema, J. Marvil, Q. A. Parker, S. K. Sirothia, N. Hurley-Walker, M. Bell, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, J. R. Callingham, A. A. Deshpande, K. S. Dwarakanath, B. -Q. For, L. J. Greenhill, P. Hancock, B. J. Hazelton , et al. (35 additional authors not shown)

Abstract: Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) s… ▽ More Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year. △ Less

Submitted 10 July, 2016; originally announced July 2016.

Comments: 25 pages, 7 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 593, A130 (2016)

Time-domain and spectral properties of pulsars at 154 MHz

Authors: M. E. Bell, Tara Murphy, S. Johnston, D. L. Kaplan, S. Croft, P. Hancock, J. R. Callingham, A. Zic, D. Dobie, J. K. Swiggum, A. Rowlinson, N. Hurley-Walker, A. R. Offringa, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell , et al. (13 additional authors not shown)

Abstract: We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < δ < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and P… ▽ More We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < δ < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715. △ Less

Submitted 30 May, 2016; originally announced May 2016.

The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology

Authors: Daniel C. Jacobs, B. J. Hazelton, C. M. Trott, Joshua S. Dillon, B. Pindor, I. S. Sullivan, J. C. Pober, N. Barry, A. P. Beardsley, G. Bernardi, Judd D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, D. Emrich, A. Ewall-Wice, L. Feng, B. M. Gaensler, R. Goeke, L. J. Greenhill, J. N. Hewitt, N. Hurley-Walker, M. Johnston-Hollitt , et al. (39 additional authors not shown)

Abstract: We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial nig… ▽ More We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Comments: accepted to ApJ

A Large Scale, Low Frequency Murchison Widefield Array Survey of Galactic HII regions between 260< l <\340

Authors: L. Hindson, M. Johnston-Hollitt, N. Hurley-Walker, J. R. Callingham, H. Su, J. Morgan, M. Bell, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, K. S. Dwarakanath, B. -Q For, B. M. Gaensler, L. J. Greenhill, P. Hancock, B. J. Hazelton, A. D. Kapinska, D. L. Kaplan, E. Lenc, C. J. Lonsdale, B. Mckinley, S. R. McWhirter, D. A. Mitchell , et al. (18 additional authors not shown)

Abstract: We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spect… ▽ More We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spectral indices of these HII regions. By identifying the point at which HII regions transition from the optically thin to thick regime we derive the physical properties including the electron density, ionised gas mass and ionising photon flux, towards 61 HII regions. This catalogue of HII regions represents the most extensive and uniform low frequency survey of HII regions in the Galaxy to date. △ Less

Submitted 21 June, 2016; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: 20 pages, 9 figures, 3 tables

First Limits on the 21 cm Power Spectrum during the Epoch of X-ray heating

Authors: A. Ewall-Wice, Joshua S. Dillon, J. N. Hewitt, A. Loeb, A. Mesinger, A. R. Neben, A. R. Offringa, M. Tegmark, N. Barry, A. P. Beardsley, G. Bernardi, Judd D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, D. Emrich, L. Feng, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, N. Hurley-Walker, M. Johnston-Hollitt , et al. (40 additional authors not shown)

Abstract: We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic e… ▽ More We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window. △ Less

Submitted 16 May, 2016; v1 submitted 29 April, 2016; originally announced May 2016.

Comments: 31 pages, 23 figures, corrected spelling errors in metadata

The 154 MHz radio sky observed by the Murchison Widefield Array: noise, confusion and first source count analyses

Authors: T. M. O. Franzen, C. A. Jackson, A. R. Offringa, R. D. Ekers, R. B. Wayth, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, J. Morgan, D. Oberoi, S. M. Ord, T. Prabu, N. Seymour , et al. (8 additional authors not shown)

Abstract: We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the c… ▽ More We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the centre of the image is 4-5 mJy/beam. The MWA counts are in excellent agreement with counts from other instruments and are the most precise ever derived in the flux density range 30-200 mJy due to the sky area covered. Using the deepest available source count data, we find that the MWA image is affected by sidelobe confusion noise at the ~3.5 mJy/beam level, due to incompletely-peeled and out-of-image sources, and classical confusion becomes apparent at ~1.7 mJy/beam. This work highlights that (i) further improvements in ionospheric calibration and deconvolution imaging techniques would be required to probe to the classical confusion limit and (ii) the shape of low-frequency source counts, including any flattening towards lower flux densities, must be determined from deeper ~150 MHz surveys as it cannot be directly inferred from higher frequency data. △ Less

Submitted 13 April, 2016; originally announced April 2016.

Comments: 13 pages, 13 figures, accepted for publication in MNRAS

High-energy sources at low radio frequency: the Murchison Widefield Array view of Fermi blazars

Authors: M. Giroletti, F. Massaro, R. D'Abrusco, R. Lico, D. Burlon, N. Hurley-Walker, M. Johnston-Hollitt, J. Morgan, V. Pavlidou, M. Bell, G. Bernardi, R. Bhat, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, A. Ewall-Rice, D. Emrich, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, L. Hindson, D. L. Kaplan , et al. (32 additional authors not shown)

Abstract: Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blaz… ▽ More Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle α_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population. △ Less

Submitted 29 February, 2016; originally announced February 2016.

Comments: accepted for publication in A&A

Limits on Fast Radio Bursts and other transient sources at 182 MHz using the Murchison Widefield Array

Authors: A. Rowlinson, M. E. Bell, T. Murphy, C. M. Trott, N. Hurley-Walker, S. Johnston, S. J. Tingay, D. L. Kaplan, D. Carbone, P. J. Hancock, L. Feng, A. R. Offringa, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales , et al. (11 additional authors not shown)

Abstract: We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for… ▽ More We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: MNRAS Accepted, 17 pages, 10 figures

Beamforming Errors in Murchison Widefield Array Phased Array Antennas and their effects on Epoch of Reionization Science

Authors: A. R. Neben, J. N. Hewitt, R. F. Bradley, J. S. Dillon, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, D. Oberoi, S. M. Ord, T. Prabu, N. Udaya Shankar, K. S. Srivani , et al. (6 additional authors not shown)

Abstract: Accurate antenna beam models are critical for radio observations aiming to isolate the redshifted 21cm spectral line emission from the Dark Ages and the Epoch of Reionization and unlock the scientific potential of 21cm cosmology. Past work has focused on characterizing mean antenna beam models using either satellite signals or astronomical sources as calibrators, but antenna-to-antenna variation d… ▽ More Accurate antenna beam models are critical for radio observations aiming to isolate the redshifted 21cm spectral line emission from the Dark Ages and the Epoch of Reionization and unlock the scientific potential of 21cm cosmology. Past work has focused on characterizing mean antenna beam models using either satellite signals or astronomical sources as calibrators, but antenna-to-antenna variation due to imperfect instrumentation has remained unexplored. We characterize this variation for the Murchison Widefield Array (MWA) through laboratory measurements and simulations, finding typical deviations of order +/- 10-20% near the edges of the main lobe and in the sidelobes. We consider the ramifications of these results for image- and power spectrum-based science. In particular, we simulate visibilities measured by a 100m baseline and find that using an otherwise perfect foreground model, unmodeled beamforming errors severely limit foreground subtraction accuracy within the region of Fourier space contaminated by foreground emission (the "wedge"). This region likely contains much of the cosmological signal, and accessing it will require measurement of per-antenna beam patterns. However, unmodeled beamforming errors do not contaminate the Fourier space region expected to be free of foreground contamination (the "EOR window"), showing that foreground avoidance remains a viable strategy. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: Accepted to ApJ

Black Hole Mass Estimation: How good is the virial estimate?

Authors: Suk Yee Yong, Rachel L. Webster, Anthea L. King

Abstract: Black hole mass is a key factor in determining how a black hole interacts with its environment. However, the determination of black hole masses at high redshifts depends on secondary mass estimators, which are based on empirical relationships and broad approximations. A dynamical disk wind broad line region (BLR) model of active galactic nuclei (AGN) is built in order to test the impact on the bla… ▽ More Black hole mass is a key factor in determining how a black hole interacts with its environment. However, the determination of black hole masses at high redshifts depends on secondary mass estimators, which are based on empirical relationships and broad approximations. A dynamical disk wind broad line region (BLR) model of active galactic nuclei (AGN) is built in order to test the impact on the black hole mass calculation due to different BLR geometries and the inclination of the AGN. Monte Carlo simulations of two disk wind models are constructed to recover the virial scale factor, $f$, at various inclination angles. The resulting $f$ values strongly correlate with inclination angle, with large $f$ values associated with small inclination angles (close to face-on) and small $f$ values with large inclination angles (close to edge-on). The $f$ factors are consistent with previously determined $f$ values, found from empirical relationships. Setting $f$ as a constant may introduce a bias into virial black hole mass estimates for a large sample of AGN. However, the extent of the bias depends on the line width characterisation (e.g. full width at half maximum (FWHM) or line dispersion). Masses estimated using $f_{\text{FWHM}}$ tend to biased towards larger masses, but this can be corrected by calibrating for the width or shape of the emission line. △ Less

Submitted 15 February, 2016; originally announced February 2016.

Comments: Accepted for publication in PASA. 9 pages, 5 figures, 2 table

Parametrising Epoch of Reionization foregrounds: A deep survey of low-frequency point-source spectra with the MWA

Authors: A. R. Offringa, C. M. Trott, N. Hurley-Walker, M. Johnston-Hollitt, B. McKinley, N. Barry, A. P. Beardsley, J. D. Bowman, F. Briggs, P. Carroll, J. S. Dillon, A. Ewall-Wice, L. Feng, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, J. N. Hewitt, D. C. Jacobs, H. -S. Kim, P. Kittiwisit, E. Lenc, J. Line, A. Loeb, D. A. Mitchell, M. F. Morales , et al. (16 additional authors not shown)

Abstract: Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the MWA EoR fields using 45 h of observation. A novel… ▽ More Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the MWA EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by PSF sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source. △ Less

Submitted 6 February, 2016; originally announced February 2016.

Comments: Accepted for publication in MNRAS. 15 pages, 16 figures. Catalogue of sources externally available as CSV file

The Importance of Wide-field Foreground Removal for 21 cm Cosmology: A Demonstration With Early MWA Epoch of Reionization Observations

Authors: J. C. Pober, B. J. Hazelton, A. P. Beardsley, N. A. Barry, Z. E. Martinot, I. S. Sullivan, M. F. Morales, M. E. Bell, G. Bernardi, N. D. R. Bhat, J. D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, A. A. Deshpande, Joshua. S. Dillon, D. Emrich, A. M. Ewall-Wice, L. Feng, R. Goeke, L. J. Greenhill, J. N. Hewitt, L. Hindson , et al. (44 additional authors not shown)

Abstract: In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements us… ▽ More In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum. △ Less

Submitted 22 January, 2016; originally announced January 2016.

Comments: 13 pages, 9 figures, matches version accepted to ApJ

CHIPS: The Cosmological HI Power Spectrum Estimator

Authors: Cathryn M. Trott, Bart Pindor, Pietro Procopio, Randall B. Wayth, Daniel A. Mitchell, Benjamin McKinley, Steven J. Tingay, N. Barry, A. P. Beardsley, G. Bernardi, Judd D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, A. de Oliveira-Costa, Joshua S. Dillon, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. J. Hazelton, J. N. Hewitt, N. Hurley-Walker, M. Johnston-Hollitt, Daniel C. Jacobs, D. L. Kaplan , et al. (27 additional authors not shown)

Abstract: Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dime… ▽ More Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$σ$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $Δ_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates. △ Less

Submitted 8 January, 2016; originally announced January 2016.

Comments: 34 pages, accepted for publication in ApJ

A search for Fast Radio Bursts at low frequencies with Murchison Widefield Array high time resolution imaging

Authors: S. J. Tingay, C. M. Trott, R. B. Wayth, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, L. Feng, B. M. Gaensler, L. J. Greenhill, P. J. Hancock, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, T. Murphy, D. Oberoi, T. Prabu, N. Udaya Shankar, K. S. Srivani , et al. (4 additional authors not shown)

Abstract: We present the results of a pilot study search for Fast Radio Bursts (FRBs) using the Murchison Widefield Array (MWA) at low frequencies (139 - 170 MHz). We utilised MWA data obtained in a routine imaging mode from observations where the primary target was a field being studied for Epoch of Reionisation detection. We formed images with 2 second time resolution and 1.28~MHz frequency resolution for… ▽ More We present the results of a pilot study search for Fast Radio Bursts (FRBs) using the Murchison Widefield Array (MWA) at low frequencies (139 - 170 MHz). We utilised MWA data obtained in a routine imaging mode from observations where the primary target was a field being studied for Epoch of Reionisation detection. We formed images with 2 second time resolution and 1.28~MHz frequency resolution for 10.5 hours of observations, over 400 square degrees of the sky. We de-dispersed the dynamic spectrum in each of 372,100 resolution elements of 2$\times$2 arcmin$^{2}$, between dispersion measures of 170 and 675~pc~cm$^{-3}$. Based on the event rate calculations in Trott, Tingay & Wayth (2013), which assumes a standard candle luminosity of $8\times10^{37}$ Js$^{-1}$, we predict that with this choice of observational parameters, the MWA should detect ($\sim10$,$\sim2$,$\sim0$) FRBs with spectral indices corresponding to ($-$2, $-$1, 0), based on a 7$σ$ detection threshold. We find no FRB candidates above this threshold from our search, placing an event rate limit of $<700$ above 700 Jy.ms per day per sky and providing evidence against spectral indices $α<-1.2$ ($S\proptoν^α$). We compare our event rate and spectral index limits with others from the literature. We briefly discuss these limits in light of recent suggestions that supergiant pulses from young neutron stars could explain FRBs. We find that such supergiant pulses would have to have much flatter spectra between 150 and 1400 MHz than have been observed from Crab giant pulses to be consistent with the FRB spectral index limit we derive. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 14 pages, 5 figures

Direct Shear Mapping: Prospects for weak lensing studies of individual galaxy-galaxy lensing systems

Authors: Catherine O. de Burgh-Day, Edward N. Taylor, Rachel L. Webster, Andrew M. Hopkins

Abstract: We have investigated, using both a theoretical and an empirical approach, the frequency of low redshift galaxy-galaxy lensing systems in which the signature of weak lensing might be directly detectable. We find good agreement between these two approaches. In order to make a theoretical estimate of the weak lensing shear, $γ$, for each galaxy in a catalogue, we have made an estimate of the asymptot… ▽ More We have investigated, using both a theoretical and an empirical approach, the frequency of low redshift galaxy-galaxy lensing systems in which the signature of weak lensing might be directly detectable. We find good agreement between these two approaches. In order to make a theoretical estimate of the weak lensing shear, $γ$, for each galaxy in a catalogue, we have made an estimate of the asymptotic circular velocity from the stellar mass using three different approaches: from a simulation based relation, from an empirically-derived relation, and using the baryonic Tully-Fisher relation. Using data from the Galaxy and Mass Assembly redshift survey we estimate the frequency of detectable weak lensing at low redshift. We find that to a redshift of $z\sim 0.6$, the probability of a galaxy being weakly lensed by at least $γ= 0.02$ is $\sim 0.01$. A scatter in the $M_*-M_h$ relation results in a shift towards higher measured shears for a given population of galaxies. Given this, and the good probability of weak lensing at low redshifts, we have investigated the feasibility of measuring the scatter in the $M_*-M_h$ relation using shear statistics. This is a novel measurement, and is made possible because DSM is able to make individual \itshape direct~\upshape shear measurements, in contrast to traditional weak lensing techniques which can only make statistical measurements. We estimate that for a shear measurement error of $Δγ= 0.02$ (consistent with the sensitivity of DSM), a sample of $\sim$50,000 spatially and spectrally resolved galaxies would allow a measurement of the scatter in the $M_*-M_h$ relation to be made. While there are no currently existing IFU surveys of this size, there are upcoming surveys which will provide this data (e.g The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), surveys with Hector, and the Square Kilometre Array (SKA)). △ Less

Submitted 20 September, 2015; originally announced September 2015.

Comments: Accepted for publication in PASA. 18 pages, 12 figures, 1 table

Quantifying ionospheric effects on time-domain astrophysics with the Murchison Widefield Array

Authors: Shyeh Tjing Loi, Tara Murphy, Martin E. Bell, David L. Kaplan, Emil Lenc, Andre R. Offringa, Natasha Hurley-Walker, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, D. Emrich, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, J. C. Kasper, E. Kratzenberg, C. J. Lonsdale, M. J. Lynch, S. R. McWhirter, D. A. Mitchell , et al. (17 additional authors not shown)

Abstract: Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) rad… ▽ More Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We directly examine 51 hours of data observed over 10 nights under quiet geomagnetic conditions (global storm index Kp < 2), analysing the behaviour of short-timescale angular position and peak flux density variations of around ten thousand unresolved sources. We find that while much of the variation in angular position can be attributed to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching under typical conditions. By examining bulk trends in amplitude variability, we place upper limits on the modulation index associated with ionospheric scintillation of 1-3% for the various nights. For sources fainter than ~1 Jy, this variation is below the image noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is not a significant impediment to the goals of time-domain science with the MWA at 154 MHz. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS