Search

A first statistical detection of the 21-cm Epoch of Reionization (EoR) is on the horizon, as cosmological volumes of the Universe become accessible via the adoption of low-frequency interferometers. We explore the impact which non-identical instrumental beam responses can have on the calibrated power spectrum and a future EoR detection. All-sky satellite measurements of Murchison Widefield Array (MWA) beams have revealed significant sidelobe deviations from cutting-edge electromagnetic simulations at the ~10% zenith power level. By generating physically motivated deformed beam models, we emulate real measurements of the MWA which inherently encode the imprints of varied beams. We explore two calibration strategies: using a single beam model across the array, or using a full set of deformed beams. Our simulations demonstrate beam-induced leakage of foreground power into theoretically uncontaminated modes, at levels which exceed the expected cosmological signal by factors of over ~1000 between the modes k=0.1-1 $hMpc^{-1}$. We also show that this foreground leakage can be mitigated by including measured models of varied beams into calibration frameworks, reducing the foreground leakage to a sub-dominant effect and potentially unveiling the EoR. Finally, we outline the future steps necessary to make this approach applicable to real measurements by radio interferometers.

PKS 2004-447 is a narrow-line Seyfert 1 (NLS1) harbouring a relativistic jet with gamma-ray emission. On 2019-10-25, the Fermi-Large Area Telescope captured a $\gamma$-ray flare from this source, offering a chance to study the broad-line region (BLR) and jet during such violent events. This can provide insights to the BLR structure and jet interactions, which are important for active galactic nuclei and host galaxy coevolution. We report X-Shooter observations of enhancements in the broad line components of Balmer, Paschen and He I lines seen only during the post-flare and vanishing 1.5 years after. These features are biased redward up to $\sim$250 km s$^{-1}$ and are narrower than the pre-existing broad line profiles. This indicates a connection between the relativistic jet and the BLR of a young AGN, and how $\gamma$-ray production can lead to localised addition of broad emission lines., Comment: 19 pages, 14 figures, 3 tables, including appendix. Accepted for publication on A&A Letters

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\sigma$ upper limit at z=6.5 of $\Delta^2 \leq (73.78 ~\mathrm{mK)^2}$ at $k=0.13~\mathrm{h~ Mpc^{-1}}$, improving on previous EoR1 measurement results., Comment: 18 pages, 17 figures, Accepted for publication in MNRAS

The locations of Ly-$\alpha$ 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 $\Delta{\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., Comment: 9 pages, 6 figures, accepted for publication in MNRAS

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., Comment: 17 pages, 18 figures, accepted to MNRAS

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., Comment: 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

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., Comment: 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"

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$\sigma$ 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., Comment: 19 pages, 29 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

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

Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $\Delta^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude., Comment: 21 pages, 10 figures, published in ApJ, updated title and text

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., Comment: 19 pages, 10 figures, accepted for publication in PASA

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., Comment: 19 pages, 18 figures, accepted for publication in MNRAS

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., Comment: 19 pages, 11 figures

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., Comment: 19 pages, 11 figures, Accepted for publication in the Astrophysical Journal

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

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., Comment: 13 pages, 7 figures, accepted in MNRAS

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., Comment: 13 pages, 7 figures Accepted to PASA

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., Comment: 20 pages, 11 figures. Accepted to ApJ

For the first time both H$\alpha$ and far-ultraviolet (FUV) observations from an HI-selected sample are used to determine the dust-corrected star formation rate density (SFRD: $\dot{\rho}$) in the local Universe. Applying the two star formation rate indicators on 294 local galaxies we determine log($\dot{\rho}$$ _{H\alpha}) = -1.68~^{+0.13}_{-0.05}$ [M$_{\odot} $ yr$^{-1} $ Mpc$^{-3}]$ and log($\dot{\rho}_{FUV}$) $ = -1.71~^{+0.12}_{-0.13}$ [M$_\odot $ yr$^{-1} $ Mpc$^{-3}]$. These values are derived from scaling H$\alpha$ 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$\alpha$/FUV flux ratios than the remaining galaxy populations, consistent with the earlier results of Meurer. The near-identical H$\alpha$- and FUV-derived SFRD values arise with the low H$\alpha$/FUV flux ratios of some galaxies being offset by enhanced H$\alpha$ 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., Comment: 17 pages, 7 figures

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., Comment: Accepted for publication in MNRAS, 10 pages, 9 figures

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-$\alpha$ 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}}$., Comment: 13pages, 13figures, 3tables, accepted to MNRAS

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)., Comment: Accepted for publication in Astrophysical Journal on 06 February 2017

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., Comment: 16 pages, 14 figures, accepted to AJ

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., Comment: Accepted for publication in Publications of the Astronomical Society of Australia (PASA) - 34 pages - 20 figures - PUMA github repository: https://github.com/JLBLine/PUMA

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., Comment: 17 pages, 10 figures, accepted for publication in The Astrophysical Journal

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 $\Delta^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., Comment: 21 pages, 16 figures, accepted to ApJ