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.

Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($\zeta$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $\zeta$ to be biased by up to $5\sigma$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1\sigma$.

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

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

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.

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 $\Delta^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$.

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.

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

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

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

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

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

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.

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.

Precise instrument calibration is critical to the success of 21 cm Cosmology experiments. Unmitigated errors in calibration contaminate the Epoch of Reionization (EoR) signal, precluding a detection. Barry et al. 2016 characterizes one class of inherent errors that emerge from calibrating to an incomplete sky model, however it has been unclear if errors in the sky model affect the calibration of redundant arrays. In this paper, we show that redundant calibration is vulnerable to errors from sky model incompleteness even in the limit of perfect antenna positioning and identical beams. These errors are at a level that can overwhelm the EoR signal and prevent a detection. Finally, we suggest error mitigation strategies with implications for the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA).

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

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

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

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

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

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.

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.

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.

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.

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

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.

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

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

We examine data from the Murchison Widefield Array (MWA) in the frequency range 72 -- 102 MHz for a field-of-view that serendipitously contained the interstellar object 'Oumuamua on 2017 November 28. Observations took place with time resolution of 0.5 s and frequency resolution of 10 kHz. %This observation was undertaken for another purpose but due to the MWA's extremely large field-of-view, 'Oumuamua was serendipitously observed simultaneously. Based on the interesting but highly unlikely suggestion that 'Oumuamua is an interstellar spacecraft, due to some unusual orbital and morphological characteristics, we examine our data for signals that might indicate the presence of intelligent life associated with 'Oumuamua. We searched our radio data for: 1) impulsive narrow-band signals; 2) persistent narrow-band signals; and 3) impulsive broadband signals. We found no such signals with non-terrestrial origins and make estimates of the upper limits on Equivalent Isotropic Radiated Power (EIRP) for these three cases of approximately 7 kW, 840 W, and 100 kW, respectively. These transmitter powers are well within the capabilities of human technologies, and are therefore plausible for alien civilizations. While the chances of positive detection in any given Search for Extraterrestrial Intelligence (SETI) experiment are vanishingly small, the characteristics of new generation telescopes such as the MWA (and in the future, the Square Kilometre Array) make certain classes of SETI experiment easy, or even a trivial by-product of astrophysical observations. This means that the future costs of SETI experiments are very low, allowing large target lists to partially balance the low probability of a positive detection.

The Square Kilometre Array (SKA) Epoch of Reionisation and Cosmic Dawn (EoR/CD) experiments aim to explore the growth of structure and production of ionising radiation in the first billion years of the Universe. Here I describe the experiments planned for the future low-frequency components of the Observatory, and work underway to define, design and execute these programs.

We review an improved statistical model of extra-galactic point-source foregrounds first introduced in [Murray 2017], in the context of the Epoch of Reionization. This model extends the instrumentally-convolved foreground covariance used in inverse-covariance foreground mitigation schemes, by considering the cosmological clustering of the sources. In this short work, we show that over scales of $k \sim (0.6, 40.) h {\rm Mpc}^{-1}$, ignoring source clustering is a valid approximation. This is in contrast to [Murray 2017], who found a possibility of false detection if the clustering was ignored. The dominant cause for this change is the introduction of a Galactic synchrotron component which shadows the clustering of sources.

In the era of the SKA precursors, telescopes are producing deeper, larger images of the sky on increasingly small time-scales. The greater size and volume of images place an increased demand on the software that we use to create catalogues, and so our source finding algorithms need to evolve accordingly. In this paper we discuss some of the logistical and technical challenges that result from the increased size and volume of images that are to be analysed, and demonstrate how the Aegean source finding package has evolved to address these challenges. In particular we address the issues of source finding on spatially correlated data, and on images in which the background, noise, and point spread function, vary across the sky. We also introduce the concept of forced or priorized fitting.

We test the hypothesis that limitations in the sky model used to calibrate an interferometric radio telescope, where the model contains extended radio sources, will generate bias in the Epoch of Reionisation (EoR) power spectrum. The information contained in a calibration model about the spatial and spectral structure of an extended source is incomplete because a radio telescope cannot sample all Fourier components. Application of an incomplete sky model to calibration of EoR data will imprint residual error in the data, which propagates forward to the EoR power spectrum. This limited information is studied in the context of current and future planned instruments and surveys at EoR frequencies, such as the Murchison Widefield Array (MWA), Giant Metrewave Radio Telescope (GMRT) and the Square Kilometre Array (SKA1-Low). For the MWA EoR experiment, we find that both the additional short baseline $uv$-coverage of the compact EoR array, and the additional long baselines provided by TGSS and planned MWA expansions, are required to obtain sufficient information on all relevant scales. For SKA1-Low, arrays with maximum baselines of 49~km and 65~km yield comparable performance at 50~MHz and 150~MHz, while 39~km, 14~km and 4~km arrays yield degraded performance.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment.

We detail new techniques for analysing ionospheric activity, using Epoch of Reionisation (EoR) datasets obtained with the Murchison Widefield Array (MWA), calibrated by the `Real-Time System' (RTS). Using the high spatial- and temporal-resolution information of the ionosphere provided by the RTS calibration solutions over 19 nights of observing, we find four distinct types of ionospheric activity, and have developed a metric to provide an `at a glance' value for data quality under differing ionospheric conditions. For each ionospheric type, we analyse variations of this metric as we reduce the number of pierce points, revealing that a modest number of pierce points is required to identify the intensity of ionospheric activity; it is possible to calibrate in real-time, providing continuous information of the phase screen. We also analyse temporal correlations, determine diffractive scales, examine the relative fractions of time occupied by various types of ionospheric activity, and detail a method to reconstruct the total electron content responsible for the ionospheric data we observe. These techniques have been developed to be instrument agnostic, useful for application on LOFAR and SKA-Low.

We describe the design and performance of the Engineering Development Array (EDA), which is a low frequency radio telescope comprising 256 dual-polarisation dipole antennas working as a phased-array. The EDA was conceived of, developed, and deployed in just 18 months via re-use of Square Kilometre Array (SKA) precursor technology and expertise, specifically from the Murchison Widefield Array (MWA) radio telescope. Using drift scans and a model for the sky brightness temperature at low frequencies, we have derived the EDA's receiver temperature as a function of frequency. The EDA is shown to be sky-noise limited over most of the frequency range measured between 60 and 240 MHz. By using the EDA in interferometric mode with the MWA, we used calibrated visibilities to measure the absolute sensitivity of the array. The measured array sensitivity matches very well with a model based on the array layout and measured receiver temperature. The results demonstrate the practicality and feasibility of using MWA-style precursor technology for SKA-scale stations. The modular architecture of the EDA allows upgrades to the array to be rolled out in a staged approach. Future improvements to the EDA include replacing the second stage beamformer with a fully digital system, and to transition to using RF-over-fibre for the signal output from first stage beamformers.

The current generation of experiments aiming to detect the neutral hydrogen signal from the Epoch of Reionisation (EoR) is likely to be limited by systematic effects associated with removing foreground sources from target fields. In this paper we develop a model for the compact foreground sources in one of the target fields of the MWA's EoR key science experiment: the `EoR1' field. The model is based on both the MWA's GLEAM survey and GMRT 150 MHz data from the TGSS survey, the latter providing higher angular resolution and better astrometric accuracy for compact sources than is available from the MWA alone. The model contains 5049 sources, some of which have complicated morphology in MWA data, Fornax A being the most complex. The higher resolution data show that 13% of sources that appear point-like to the MWA have complicated morphology such as double and quad structure, with a typical separation of 33~arcsec. We derive an analytic expression for the error introduced into the EoR two-dimensional power spectrum due to peeling close double sources as single point sources and show that for the measured source properties, the error in the power spectrum is confined to high $k_\bot$ modes that do not affect the overall result for the large-scale cosmological signal of interest. The brightest ten mis-modelled sources in the field contribute 90% of the power bias in the data, suggesting that it is most critical to improve the models of the brightest sources. With this hybrid model we reprocess data from the EoR1 field and show a maximum of 8% improved calibration accuracy and a factor of two reduction in residual power in $k$-space from peeling these sources. Implications for future EoR experiments including the SKA are discussed in relation to the improvements obtained.

We present a sophisticated statistical point-source foreground model for low-frequency radio Epoch of Reionization (EoR) experiments using the 21 cm neutral hydrogen emission line. Motivated by our understanding of the low-frequency radio sky, we enhance the realism of two model components compared with existing models: the source count distributions as a function of flux density and spatial position (source clustering), extending current formalisms for the foreground covariance of 2D power spectral modes in 21 cm EoR experiments. The former we generalise to an arbitrarily broken power-law, and the latter to an arbitrary isotropically-correlated field. This paper presents expressions for the modified covariance under these extensions, and shows that for a more realistic source spatial distribution, extra covariance arises in the EoR window which was previously unaccounted for. Failure to include this contribution can yield bias in the final power spectrum and under-estimate uncertainties, potentially leading to a false detection of signal. The extent of this effect is uncertain, owing to ignorance of physical model parameters, but we show that it is dependent on the relative abundance of faint sources, to the effect that our extension will become more important for future deep surveys. Finally, we show that under some parameter choices, ignoring source clustering can lead to false detections on large scales, due to both the induced bias and an artificial reduction in the estimated measurement uncertainty.

We test the bandpass smoothness performance of two prototype Square Kilometre Array (SKA) SKA1-Low log-periodic dipole antennas, the SKALA2 and SKALA3 (`SKA Log-periodic Antenna'), and the current dipole from the Murchison Widefield Array (MWA) precursor telescope. Throughout this paper, we refer to the output complex-valued voltage response of an antenna when connected to a low noise amplifier (LNA), as the dipole bandpass. In Paper I (de Lera Acedo et al. 2017), the bandpass spectral response of the log-periodic antenna being developed for the SKA1-Low was estimated using numerical electromagnetic simulations and analyzed using low-order polynomial fittings and it was compared with the HERA antenna against the delay spectrum metric. In this work, realistic simulations of the SKA1-Low instrument, including frequency-dependent primary beams and array configuration, are used with a weighted least-squares polynomial estimator to assess the ability of prototype antennas to perform the SKA Epoch of Reionisation (EoR) statistical experiments. This work complements the ideal estimator tolerances computed for the proposed EoR science experiments in Trott & Wayth (2016), with the realised performance of an optimal and standard estimation (calibration) procedure. With a sufficient sky calibration model at higher frequencies, all antennas have bandpasses that are sufficiently smooth to meet the tolerances described in Trott & Wayth (2016) to perform the EoR statistical experiments, and these are primarily limited by an adequate sky calibration model, and the thermal noise level in the calibration data. At frequencies of the Cosmic Dawn (CD), which is of principal interest to SKA as one of the first next-generation telescopes capable of accessing higher redshifts, the MWA dipole and SKALA3 antenna have adequate performance, while the SKALA2 design will impede the ability to explore this era.

This paper is the first in a series of papers describing the impact of antenna instrumental artefacts on the 21-cm cosmology experiments to be carried out by the low frequency instrument (SKA1-LOW) of the Square Kilometre Array telescope (SKA), i.e., the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The smoothness of the passband response of the current log-periodic antenna being developed for the SKA1-LOW is analyzed using numerical electromagnetic simulations. The amplitude variations over the frequency range are characterized using low-order polynomials defined locally, in order to study the impact of the passband smoothness in the instrument calibration and CD/EoR Science. A solution is offered to correct a fast ripple found at 60~MHz during a test campaign at the SKA site at the Murchison Radio-astronomy Observatory, Western Australia in September 2015 with a minor impact on the telescope's performance and design. A comparison with the Hydrogen Epoch of Reionization Array antenna is also shown demonstrating the potential use of the SKA1-LOW antenna for the Delay Spectrum technique to detect the EoR.

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

We study the signatures of reionization and ionizing properties of the early galaxies in the cross-correlations between the 21cm emission from the spin-flip transition of neutral hydrogen (HI) and the underlying galaxy population, in particular a sub-population of galaxies visible as Lyman Alpha Emitters (LAEs). With both observables simultaneously derived from a $z\simeq6.6$ hydrodynamical simulation (GADGET-2) snapshot post-processed with a radiative transfer code (pCRASH) and a dust model, we perform a parameter study and aim to constrain both the average intergalactic medium (IGM) ionization state ($1-\langle \chi_{HI} \rangle$) and the reionization topology (outside-in versus inside-out). We find that in our model LAEs occupy the densest and most-ionized regions resulting in a very strong anti-correlation between the LAEs and the 21cm emission. A 1000h SKA-LOW1 - Subaru Hyper Suprime Cam experiment can provide constraints on $\langle \chi_{HI} \rangle$, allowing us to distinguish between IGM ionization levels of 50%, 25%, 10% and fully ionized at scales $r<10$ comoving Mpc (assuming foreground avoidance for SKA). Our results support the inside-out reionization scenario where the densest knots (under-dense voids) are ionized first (last) for $\langle \chi_{HI} \rangle >= 0.1$. Further, 1000h SKA-LOW1 observations should be able to confirm the inside-out scenario by detecting a lower 21cm brightness temperature (by about 2-10 mK) in the densest regions ($> 2$ arcminute scales) hosting LAEs compared to lower-density regions devoid of them.

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.

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.

We present and evaluate several strategies to search for prompt, low-frequency radio emission associated with gravitational wave transients using the Murchison Widefield Array (MWA). As we are able to repoint the MWA on timescales of tens of seconds, we can search for the dispersed radio signal that has been predicted to originate along with or shortly after a neutron star-neutron star merger. We find that given the large, 600 deg^2 instantaneous field-of-view of the MWA we can cover a significant fraction of the predicted gravitational wave error region, although due to the complicated geometry of the latter we only cover >50% of the error region for approximately 5% of events, and roughly 15% of events will be located <10 deg from the MWA pointing center such that they will be covered in the radio images. For optimal conditions our limiting flux density for a 10-s long transient would be 0.1 Jy, increasing to about 1 Jy for a wider range of events. This corresponds to luminosity limits of 1e38-1e39 erg/s based on expectations for the distances of the gravitational wave transients, which should be sufficient to detect or significantly constrain a range of models for prompt emission.

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.

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.

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.

The V-FASTR experiment on the Very Long Baseline Array was designed to detect dispersed pulses of milliseconds duration, such as fast radio bursts (FRBs). We use all V-FASTR data through February 2015 to report V-FASTR's upper limits on the rates of FRBs, and compare these with re-derived rates from Parkes FRB detection experiments. V-FASTR's operation at lambda=20 cm allows direct comparison with the 20 cm Parkes rate, and we derive a power-law limit of \gamma<-0.4 (95% confidence limit) on the index of FRB source counts, N(>S)\propto S^\gamma. Using the previously measured FRB rate and the unprecedented amount of survey time spent searching for FRBs at a large range of wavelengths (0.3 cm > \lambda > 90 cm), we also place frequency-dependent limits on the spectral distribution of FRBs. The most constraining frequencies place two-point spectral index limits of \alpha_{20cm}^{4cm} < 5.8 and \alpha_{90cm}^{20cm} > -7.6, where fluence F \propto f^\alpha if we assume true the burst rate reported by Champion et al. (2016) of R(F~0.6 Jy ms) = 7 x 10^3 sky^{-1} day^{-1} (for bursts of ~3 ms duration). This upper limit on \alpha suggests that if FRBs are extragalactic but non-cosmological, that on average they are not experiencing excessive free-free absorption due to a medium with high optical depth (assuming temperature ~8,000 K), which excessively invert their low-frequency spectrum. This in turn implies that the dispersion of FRBs arises in either or both of the intergalactic medium or the host galaxy, rather than from the source itself.

The Cosmic Dawn and Epoch of Reionisation, during which collapsed structures produce the first ionising photons and proceed to reionise the intergalactic medium, span a large range in redshift (z~30-6) and time (t_{age} ~ 0.1-1.0~Gyr). Exploration of these epochs using the redshifted 21~cm emission line from neutral hydrogen is currently limited to statistical detection and estimation metrics (e.g., the power spectrum) due to the weakness of the signal. Brightness temperature fluctuations in the line-of-sight (LOS) dimension are probed by observing the emission line at different frequencies, and their structure is used as a primary discriminant between the cosmological signal and contaminating foreground extragalactic and Galactic continuum emission. Evolution of the signal over the observing bandwidth leads to the `line cone effect' whereby the HI structures at the start and end of the observing band are not statistically consistent, yielding a biased estimate of the signal power, and potential reduction in signal detectability. We implement a wavelet transform to wide bandwidth radio interferometry experiments to probe the local statistical properties of the signal. We show that use of the wavelet transform yields estimates with improved estimation performance, compared with the standard Fourier Transform over a fixed bandwidth. With the suite of current and future large bandwidth reionisation experiments, such as with the 300~MHz instantaneous bandwidth of the Square Kilometre Array, a transform that retains local information will be important.

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.

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.

Electron density irregularities in the ionosphere are known to be magnetically anisotropic, preferentially elongated along the lines of force. While many studies of their morphology have been undertaken by topside sounding and whistler measurements, it is only recently that detailed regional-scale reconstructions have become possible, enabled by the advent of widefield radio telescopes. Here we present a new approach for visualising and studying field-aligned irregularities (FAIs), which involves transforming interferometric measurements of TEC gradients onto a magnetic shell tangent plane. This removes the perspective distortion associated with the oblique viewing angle of the irregularities from the ground, facilitating the decomposition of dynamics along and across magnetic field lines. We apply this transformation to the dataset of Loi et al. [2015a], obtained on 15 October 2013 by the Murchison Widefield Array (MWA) radio telescope and displaying prominent FAIs. We study these FAIs in the new reference frame, quantifying field-aligned and field-transverse behaviour, examining time and altitude dependencies, and extending the analysis to FAIs on sub-array scales. We show that the inclination of the plane can be derived solely from the data, and verify that the best-fit value is consistent with the known magnetic inclination. The ability of the model to concentrate the fluctuations along a single spatial direction may find practical application to future calibration strategies for widefield interferometry, by providing a compact representation of FAI-induced distortions.

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.

Spectral features introduced by instrumental chromaticity of radio interferometers have the potential to negatively impact the ability to perform Epoch of Reionisation (EoR) and Cosmic Dawn (CD) science using the redshifted neutral hydrogen emission line from the early Universe. We describe instrument calibration choices that influence the spectral characteristics of the science data, and assess their impact on EoR statistical and tomographic experiments. Principally, we consider the intrinsic spectral response of the receiving antennas, embedded within a complete frequency-dependent primary beam response, and frequency-dependent instrument sampling. We assess different options for bandpass calibration. The analysis is applied to the proposed SKA1-Low EoR/CD experiments. We provide tolerances on the smoothness of the SKA station primary beam bandpass, to meet the scientific goals of statistical and tomographic (imaging) EoR programs. Two calibration strategies are tested: (1) fitting of each fine channel independently, and (2) fitting of an nth-order polynomial for each ~1~MHz coarse channel with (n+1)th-order residuals (n=2,3,4). Strategy (1) leads to uncorrelated power in the 2D power spectrum proportional to the thermal noise power, thereby reducing the overall array sensitivity. Strategy (2) leads to correlated residuals from the fitting, and residual signal power with (n+1)th-order curvature. For the residual power to be less than the thermal noise, the fractional amplitude of a fourth-order term in the bandpass across a single coarse channel must be <2.5% (50~MHz), <0.5% (150~MHz), <0.8% (200~MHz). The tomographic experiment places stringent constraints on phase residuals in the bandpass. We find that the root-mean-square variability over all stations of the change in phase across any fine channel (4.578~kHz) should not exceed 0.2 degrees.

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.

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.

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.

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$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

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$\sigma$ 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 $\alpha<-1.2$ ($S\propto\nu^{\alpha}$). 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.

Interplanetary Scintillation (IPS) induces intensity fluctuations in small angular size astronomical radio sources via the distortive effects of spatially and temporally varying electron density associated with outflows from the Sun. These radio sources are a potential foreground contaminant signal for redshifted HI emission from the Epoch of Reionisation (EoR) because they yield time-dependent flux density variations in bright extragalactic point sources. Contamination from foreground continuum sources complicates efforts to discriminate the cosmological signal from other sources in the sky. In IPS, at large angles from the Sun applicable to EoR observations, weak scattering induces spatially and temporally correlated fluctuations in the measured flux density of sources in the field, potentially affecting the detectability of the EoR signal by inducing non-static variations in the signal strength. In this work, we explore the impact of interplanetary weak scintillation on EoR power spectrum measurements, accounting for the instrumental spatial and temporal sampling. We use published power spectra of electron density fluctuations and parameters of EoR experiments to derive the IPS power spectrum in the wavenumber phase space of EoR power spectrum measurements. The contrast of IPS power to expected cosmological power is used as a metric to assess the impact of IPS. We show that IPS has a different spectral structure to power from foregrounds alone, but the additional leakage into the EoR observation parameter space is negligible under typical IPS conditions, unless data are used from deep within the foreground contamination region.

We present observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to >1e6 km. The serendipitous night-time nature of these detections illustrates the new regime that the MWA has opened for IPS studies with sensitive night-time, wide-field, low-frequency observations. This regime complements traditional dedicated strategies for observing IPS and can be utilized in real-time to facilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact (arcsec) structures in low-frequency radio sources despite the 2 arcmin resolution of the array.

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celestial sources are proportional to spatial gradients in the electron column density transverse to the line of sight. These can be used to probe plasma structures and waves in the ionosphere. The regional (10-100 km) scales probed by the MWA, determined by the size of its FoV and the spatial density of radio sources (typically thousands in a single FoV), complement the global (100-1000 km) scales of GPS studies and local (0.01-1 km) scales of radar scattering measurements. Our data exhibit a range of complex structures and waves. Some fluctuations have the characteristics of travelling ionospheric disturbances (TIDs), while others take the form of narrow, slowly-drifting bands aligned along the Earth's magnetic field.

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40 kHz frequency resolution and 0.5 s time resolution; the second year using 10 kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing and image weighting scheme. At 154 MHz the image resolution is approximately 2.5 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA's primary beam.

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.

Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures.

We present a three epoch survey for transient and variables in the extended Chandra Deep Field South at 5.5 GHz with the Australia Telescope Compact Array. A region covering $\sim$0.3 deg$^{2}$ was observed on timescales of 2.5 months and 2.5 years and typical sensitivities 12.1 $-$ 17.1 $\mu$Jy beam$^{-1}$ (1$\sigma$) were achieved. This survey represents the deepest search for transient and variable radio sources at 5.5 GHz. In total 124 sources were detected above the 5.5$\sigma$ level. One highly variable radio source was found with $\Delta S > 50%$ implying a surface density of $\sim$3 deg$^{-2}$. A further three radio sources were found with lower levels of variability equating to a surface density of $\sim$13 deg$^{-2}$ above a detection threshold of 82.3 $\mu$Jy. All of the variable sources have inverted radio spectra (between 1.4 and 5.5 GHz) and are associated with active galactic nuclei. We conclude that these variables are young gigahertz peaked-spectrum sources with active and self-absorbed radio jets. We explore the variability completeness of this sample and conclude that the fairly low levels of variability would only be detectable in 3$-$25% of all sources within the field. No radio transients were detected in this survey and we place an upper limit on the surface density of transient events $ < 7.5$ deg$^{-2}$ above a detection threshold of 68.8~$\mu$Jy.

Impulsive radio bursts that are detectable across cosmological distances constitute extremely powerful probes of the ionized Inter-Galactic Medium (IGM), intergalactic magnetic fields, and the properties of space-time itself. Their dispersion measures (DMs) will enable us to detect the "missing" baryons in the low-redshift Universe and make the first measurements of the mean galaxy halo profile, a key parameter in models of galaxy formation and feedback. Impulsive bursts can be used as cosmic rulers at redshifts exceeding 2, and constrain the dark energy equation-of-state parameter, $w(z)$ at redshifts beyond those readily accessible by Type Ia SNe. Both of these goals are realisable with a sample of $\sim 10^4$ fast radio bursts (FRBs) whose positions are localized to within one arcsecond, sufficient to obtain host galaxy redshifts via optical follow-up. It is also hypothesised that gravitational wave events may emit coherent emission at frequencies probed by SKA1-LOW, and the localization of such events at cosmological distances would enable their use as cosmological standard sirens. To perform this science, such bursts must be localized to their specific host galaxies so that their redshifts may be obtained and compared against their dispersion measures, rotation measures, and scattering properties. The SKA can achieve this with a design that has a wide field-of-view, a substantial fraction of its collecting area in a compact configuration (80\% within a 3\,km radius), and a capacity to attach high-time-resolution instrumentation to its signal path.

The Murchison Widefield Array (MWA) is a new low-frequency interferometric radio telescope built in Western Australia at one of the locations of the future Square Kilometre Array (SKA). We describe the automated radio-frequency interference (RFI) detection strategy implemented for the MWA, which is based on the AOFlagger platform, and present 72-231-MHz RFI statistics from 10 observing nights. RFI detection removes 1.1% of the data. RFI from digital TV (DTV) is observed 3% of the time due to occasional ionospheric or atmospheric propagation. After RFI detection and excision, almost all data can be calibrated and imaged without further RFI mitigation efforts, including observations within the FM and DTV bands. The results are compared to a previously published Low-Frequency Array (LOFAR) RFI survey. The remote location of the MWA results in a substantially cleaner RFI environment compared to LOFAR's radio environment, but adequate detection of RFI is still required before data can be analysed. We include specific recommendations designed to make the SKA more robust to RFI, including: the availability of sufficient computing power for RFI detection; accounting for RFI in the receiver design; a smooth band-pass response; and the capability of RFI detection at high time and frequency resolution (second and kHz-scale respectively).

Using the new wideband capabilities of the Australia Telescope Compact Array (ATCA), we obtain spectra for PKS 1718-649, a well-known gigahertz-peaked spectrum radio source. The observations, between approximately 1 and 10 GHz over three epochs spanning approximately 21 months, reveal variability both above the spectral peak at ~3 GHz and below the peak. The combination of the low and high frequency variability cannot be easily explained using a single absorption mechanism, such as free-free absorption or synchrotron self-absorption. We find that the PKS 1718-649 spectrum and its variability are best explained by variations in the free-free optical depth on our line-of-sight to the radio source at low frequencies (below the spectral peak) and the adiabatic expansion of the radio source itself at high frequencies (above the spectral peak). The optical depth variations are found to be plausible when X-ray continuum absorption variability seen in samples of Active Galactic Nuclei is considered. We find that the cause of the peaked spectrum in PKS 1718-649 is most likely due to free-free absorption. In agreement with previous studies, we find that the spectrum at each epoch of observation is best fit by a free-free absorption model characterised by a power-law distribution of free-free absorbing clouds. This agreement is extended to frequencies below the 1 GHz lower limit of the ATCA by considering new observations with Parkes at 725 MHz and 199 MHz observations with the newly operational Murchison Widefield Array. These lower frequency observations argue against families of absorption models (both free-free and synchrotron self-absorption) that are based on simple homogenous structures.

For there to be a successful measurement of the 21 cm Epoch of Reionization (EoR) power spectrum, it is crucial that strong foreground contaminants be robustly suppressed. These foregrounds come from a variety of sources (such as Galactic synchrotron emission and extragalactic point sources), but almost all share the property of being spectrally smooth, and when viewed through the chromatic response of an interferometer, occupy a signature "wedge" region in cylindrical $k_\perp k_\parallel$ Fourier space. The complement of the foreground wedge is termed the "EoR window", and is expected to be mostly foreground-free, allowing clean measurements of the power spectrum. This paper is a sequel to a previous paper that established a rigorous mathematical framework for describing the foreground wedge and the EoR window. Here, we use our framework to explore statistical methods by which the EoR window can be enlarged, thereby increasing the sensitivity of a power spectrum measurement. We adapt the FKP approximation (commonly used in galaxy surveys) for 21 cm cosmology, and also compare the optimal quadratic estimator to simpler estimators that ignore covariances between different Fourier modes. The optimal quadratic estimator is found to suppress foregrounds by an extra factor of $\sim 10^5$ in power at the peripheries of the EoR window, boosting the detection of the cosmological signal from $12\sigma$ to $50\sigma$ at the midpoint of reionization in our fiducial models. If numerical issues can be finessed, decorrelation techniques allow the EoR window to be further enlarged, enabling measurements to be made deep within the foreground wedge. These techniques do not assume that foreground are Gaussian-distributed, and we additionally prove that a final round of foreground subtraction can be performed after decorrelation in a way that is guaranteed to have no cosmological signal loss.

The 21 cm line provides a powerful probe of astrophysics and cosmology at high redshifts, but unlocking the potential of this probe requires the robust mitigation of foreground contaminants that are typically several orders of magnitude brighter than the cosmological signal. Recent simulations and observations have shown that the smooth spectral structure of foregrounds combines with instrument chromaticity to contaminate a "wedge"-shaped region in cylindrical Fourier space. While previous efforts have explored the suppression of foregrounds within this wedge, as well as the avoidance of this highly contaminated region, all such efforts have neglected a rigorous examination of the error statistics associated with the wedge. Using a quadratic estimator formalism applied to the interferometric measurement equation, we provide a framework for such a rigorous analysis (incorporating a fully covariant treatment of errors). Additionally, we find that there are strong error correlations at high spatial wavenumbers that have so far been neglected in sensitivity derivations. These error correlations substantially degrade the sensitivity of arrays relying on contributions from long baselines, compared to what one would estimate assuming uncorrelated errors.

The detection of six Fast Radio Bursts (FRBs) has recently been reported. FRBs are short duration ($\sim$ 1 ms), highly dispersed radio pulses from astronomical sources. The physical interpretation for the FRBs remains unclear but is thought to involve highly compact objects at cosmological distance. It has been suggested that a fraction of FRBs could be physically associated with gamma-ray bursts (GRBs). Recent radio observations of GRBs have reported the detection of two highly dispersed short duration radio pulses using a 12 m radio telescope at 1.4 GHz. Motivated by this result, we have performed a systematic and sensitive search for FRBs associated with GRBs. We have observed five GRBs at 2.3 GHz using a 26 m radio telescope located at the Mount Pleasant Radio Observatory, Hobart. The radio telescope was automated to rapidly respond to Gamma-ray Coordination Network notifications from the Swift satellite and slew to the GRB position within $\sim$ 140 s. The data were searched for pulses up to 5000 pc $\rm cm^{-3}$ in dispersion measure and pulse widths ranging from 640 $\rm \mu$s to 25.60 ms. We did not detect any events $\rm \geq 6 \sigma$. An in-depth statistical analysis of our data shows that events detected above $\rm 5 \sigma$ are consistent with thermal noise fluctuations only. A joint analysis of our data with previous experiments shows that previously claimed detections of FRBs from GRBs are unlikely to be astrophysical. Our results are in line with the lack of consistency noted between the recently presented FRB event rates and GRB event rates.

(Abridged) Noise considerations for experiments that aim to statistically estimate the 21 cm signal from high redshift neutral hydrogen during the Epoch of Reionisation (EoR) using interferometric data are typically computed assuming a tracked observation, where the telescope pointing centre and instrument phase centre are the same over the observation. Current low frequency interferometers use aperture arrays of fixed dipoles, which are steered electronically on the sky, and have different properties to mechanically-steered single apertures, such as reduced sensitivity away from zenith, and discrete pointing positions on the sky. Use of two additional observing modes is encouraged: (1) zenith drift, where the pointing centre remains fixed at the zenith, and the phase centre tracks the sky, and (2) drift+shift, where the telescope uses discrete pointing centres, and the sky drifts during each pointing. The coherence of measurements in these modes dictates the optimal reduction in thermal noise by combination of coherent modes, and the reduction in cosmic variance by combination of incoherent modes (views of different patches of the sky). The balance between these noise components provides one measure for the utility of these three modes for measuring a statistical EoR signature. We provide a general framework for estimating the uncertainty in the power spectrum for a given observing mode, telescope beam shape, and antenna distribution. We apply this framework to the Murchison Widefield Array (MWA) using an analysis of the 2D and 1D power spectra for 900 hours of observing. We demonstrate that zenith drift scans can yield marginally lower uncertainty in the signal power compared with tracked scans for the MWA EoR experiment.

We present a search for transient and variable radio sources at 154 MHz with the Murchison Widefield Array 32-tile prototype. Fifty-one images were obtained that cover a field of view of 1430 deg^2 centred on Hydra A. The observations were obtained over three days in 2010 March and three days in 2011 April and May. The mean cadence of the observations was 26 minutes and there was additional temporal information on day and year timescales. We explore the variability of a sample of 105 low frequency radio sources within the field. Four bright (S > 6 Jy) candidate variable radio sources were identified that displayed low levels of short timescale variability (26 minutes). We conclude that this variability is likely caused by simplifications in the calibration strategy or ionospheric effects. On the timescale of one year we find two sources that show significant variability. We attribute this variability to either refractive scintillation or intrinsic variability. No radio transients were identified and we place an upper limit on the surface density of sources rho < 7.5 x 10^-5 deg^-2 with flux densities > 5.5 Jy, and characteristic timescales of both 26 minutes and one year.

Fast Radio Bursts (FRBs) are short timescale (<<1 s) astrophysical radio signals, presumed to be a signature of cataclysmic events of extragalactic origin. The discovery of six high-redshift events at ~1400 MHz from the Parkes radio telescope suggests that FRBs may occur at a high rate across the sky. The Murchison Widefield Array (MWA) operates at low radio frequencies (80-300 MHz) and is expected to detect FRBs due to its large collecting area (~2500 m^2) and wide field-of-view (FOV, ~1000 square degrees at nu=200 MHz). We compute the expected number of FRB detections for the MWA assuming a source population consistent with the reported detections. Our formalism properly accounts for the frequency-dependence of the antenna primary beam, the MWA system temperature, and unknown spectral index of the source population, for three modes of FRB detection: coherent; incoherent; and fast imaging. We find that the MWA's sensitivity and large FOV combine to provide the expectation of multiple detectable events per week in all modes, potentially making it an excellent high time resolution science instrument. Deviations of the expected number of detections from actual results will provide a strong constraint on the assumptions made for the underlying source population and intervening plasma distribution.

The VLBA was used at 1.6 GHz to observe a target field 50 arcminutes in diameter including the core of M31. Novel VLBI correlation techniques were used to observe 200 sources simultaneously, of which 16 were detected. We classify all 16 as background AGN based on their X-ray properties and arcsecond- and mas-scale morphology. The detected sources were then analyzed for evidence of scatter-broadening due to the ionized ISM of M31. The detection of a compact background source only 0.25 kpc projected distance from M31 places a constraint on the extent of any extreme scattering region associated with center of M31. However, the two sources closest to the core show evidence of scatter broadening consistent with that which would be observed if a compact source were observed through the inner disk of our Galaxy at the inclination of M31. We interpret this as a detection of the ionized ISM of M31 along two lines of sight. With the increases in bandwidth and sensitivity envisaged for future long-baseline interferometers, this should prove to be a remarkably powerful technique for understanding the ionized ISM in external galaxies.

We define a framework for determining constraints on the detection rate of fast transient events from a population of underlying sources, with a view to incorporating beam shape, frequency effects, scattering effects, and detection efficiency into the metric. We then demonstrate a method for combining independent datasets into a single event rate constraint diagram, using a probabilistic approach to the limits on parameter space. We apply this new framework to present the latest results from the V-FASTR experiment, a commensal fast transients search using the Very Long Baseline Array (VLBA). In the 20 cm band, V-FASTR now has the ability to probe the regions of parameter space of importance for the observed Lorimer and Keane fast radio transient candidates, by combining the information from observations with differing bandwidths, and properly accounting for the source dispersion measure, VLBA antenna beam shape, experiment time sampling, and stochastic nature of events. We then apply the framework to combine the results of the V-FASTR and ATA Fly's Eye experiments, demonstrating their complementarity. Expectations for fast transients experiments for the SKA Phase I dish array are then computed, and the impact of large differential bandwidths is discussed.

We investigate the S/N of a new incoherent dedispersion algorithm optimized for FPGA-based architectures intended for deployment on ASKAP and other SKA precursors for fast transients surveys. Unlike conventional CPU- and GPU-optimized incoherent dedispersion algorithms, this algorithm has the freedom to maximize the S/N by way of programmable dispersion profiles that enable the inclusion of different numbers of time samples per spectral channel. This allows, for example, more samples to be summed at lower frequencies where intra-channel dispersion smearing is larger, or it could even be used to optimize the dedispersion sum for steep spectrum sources. Our analysis takes into account the intrinsic pulse width, scatter broadening, spectral index and dispersion measure of the signal, and the system's frequency range, spectral and temporal resolution, and number of trial dedispersions. We show that the system achieves better than 80% of the optimal S/N where the temporal resolution and the intra-channel smearing time are smaller than a quarter of the average width of the pulse across the system's frequency band (after including scatter smearing). Coarse temporal resolutions suffer a Delta_t^(-1/2) decay in S/N, and coarse spectral resolutions cause a Delta_nu^(-1/2) decay in S/N, where Delta_t and Delta_nu are the temporal and spectral resolutions of the system, respectively. We show how the system's S/N compares with that of matched filter and boxcar filter detectors. We further present a new algorithm for selecting trial dispersion measures for a survey that maintains a given minimum S/N performance across a range of dispersion measures.

We use information theory to derive fundamental limits on the capacity to calibrate next-generation radio interferometers, and measure parameters of point sources for instrument calibration, point source subtraction, and data deconvolution. We demonstrate the implications of these fundamental limits, with particular reference to estimation of the 21cm Epoch of Reionization power spectrum with next-generation low-frequency instruments (e.g., the Murchison Widefield Array -- MWA, Precision Array for Probing the Epoch of Reionization -- PAPER), where short time scale instrumental calibration is required due to the impact of the ionosphere on the signal wavefront. Finally, we explore the optimal point source precision available by using a combination of current and prior information. Estimation schemes that incorporate prior information may be advantageous when the measurement precision is comparable to the characteristic refraction scale of the ionosphere.

Precise subtraction of foreground sources is crucial for detecting and estimating 21cm HI signals from the Epoch of Reionization (EoR). We quantify how imperfect point source subtraction due to limitations of the measurement dataset yields structured residual signal in the dataset. We use the Cramer-Rao lower bound, as a metric for quantifying the precision with which a parameter may be measured, to estimate the residual signal in a visibility dataset due to imperfect point source subtraction. We then propagate these residuals into two metrics of interest for 21cm EoR experiments - the angular and two-dimensional power spectrum - using a combination of full analytic covariant derivation, analytic variant derivation, and covariant Monte Carlo simulations. This methodology differs from previous work in two ways: (1) it uses information theory to set the point source position error, rather than assuming a global root-mean-square error, and (2) it describes a method for propagating the errors analytically, thereby obtaining the full correlation structure of the power spectra. The methods are applied to two upcoming low-frequency instruments: the Murchison Widefield Array and the Precision Array for Probing the Epoch of Reionization. In addition to the actual antenna configurations, we apply the methods to minimally-redundant and maximally-redundant configurations. We find that for peeling sources above 1 Jy, the amplitude of the residual signal, and its variance, will be smaller than the contribution from thermal noise for the observing parameters proposed for upcoming EoR experiments, and that optimal subtraction of bright point sources will not be a limiting factor for EoR parameter estimation. We then use the formalism to provide an ab initio analytic derivation motivating the 'wedge' feature in the two-dimensional power spectrum, complementing previous discussion in the literature.

The first Search for Extra-Terrestrial Intelligence (SETI) conducted with Very Long Baseline Interferometry (VLBI) is presented. By consideration of the basic principles of interferometry, we show that VLBI is efficient at discriminating between SETI signals and human generated radio frequency interference (RFI). The target for this study was the star Gliese 581, thought to have two planets within its habitable zone. On 2007 June 19, Gliese 581 was observed for 8 hours at 1230-1544 with the Australian Long Baseline Array. The dataset was searched for signals appearing on all interferometer baselines above five times the noise limit. A total of 222 potential SETI signals were detected and by using automated data analysis techniques, were ruled out as originating from the Gliese 581 system. From our results we place an upper limit of 7 MW/Hz on the power output of any isotropic emitter located in the Gliese 581 system, within this frequency range. This study shows that VLBI is ideal for targeted SETI, including follow-up observations. The techniques presented are equally applicable to next-generation interferometers, such as the long baselines of the Square Kilometre Array (SKA).

Transient radio signals of astrophysical origin present an avenue for studying the dynamic universe. With the next generation of radio interferometers being planned and built, there is great potential for detecting and studying large samples of radio transients. Currently-used image-based techniques for detecting radio sources have not been demonstrated to be optimal, and there is a need for development of more sophisticated algorithms, and methodology for comparing different detection techniques. A visibility-space detector benefits from our good understanding of visibility-space noise properties, and does not suffer from the image artifacts and need for deconvolution in image-space detectors. In this paper, we propose a method for designing optimal source detectors using visibility data, building on statistical decision theory. The approach is substantially different to conventional radio astronomy source detection. Optimal detection requires an accurate model for the data, and we present a realistic model for the likelihood function of radio interferometric data, including the effects of calibration, signal confusion and atmospheric phase fluctuations. As part of this process, we derive fundamental limits on the calibration of an interferometric array, including the case where many relatively weak "in-beam" calibrators are used. These limits are then applied, along with a model for atmospheric phase fluctuations, to determine the limits on measuring source position, flux density and spectral index, in the general case. We then present an optimal visibility-space detector using realistic models for an interferometer.

We construct a mass model for the spiral lens galaxy 2237+0305, at redshift z_l=0.04, based on gravitational-lensing constraints, HI rotation, and new stellar-kinematic information, based on data taken with the ESI spectrograph on the 10m Keck-II Telescope. High resolution rotation curves and velocity dispersion profiles along two perpendicular directions, close to the major and minor axes of the lens galaxy, were obtained by fitting the Mgb-Fe absorption line region. The stellar rotation curve rises slowly and flattens at r~1.5" (~1.1 kpc). The velocity dispersion profile is approximately flat. A combination of photometric, kinematic and lensing information is used to construct a mass model for the four major mass components of the system -- the dark matter halo, disc, bulge, and bar. The best-fitting solution has a dark matter halo with a logarithmic inner density slope of gamma=0.9+/-0.3 for rho_DM propto r^-gamma, a bulge with M/L_B=6.6+/-0.3 Upsilon_odot, and a disc with M/L_B =1.2+/-0.3 Upsilon_odot, in agreement with measurements of late-type spirals. The bulge dominates support in the inner regions where the multiple images are located and is therefore tightly constrained by the observations. The disc is sub-maximal and contributes 45+/-11 per cent of the rotational support of the galaxy at 2.2r_d. The halo mass is (2.0+/-0.6) x 10^12 M_odot, and the stellar to virial mass ratio is 7.0+/-2.3 per cent, consistent with typical galaxies of the same mass.

The statistical mechanics of N cold dark matter (CDM) particles interacting via a softened gravitational potential is reviewed in the microcanonical ensemble and mean-field limit. A phase diagram for the system is computed as a function of the total energy E and gravitational softening length epsilon. For softened systems, two stable phases exist: a collapsed phase, whose radial density profile rho(r) is a central Dirac cusp, and an extended phase, for which rho(r) has a central core and rho(r) ~ r^{-2.2} at large r. It is shown that many N-body simulations of CDM haloes in the literature inadvertently sample the collapsed phase only, even though this phase is unstable when there is zero softening. Consequently, there is no immediate reason to expect agreement between simulated and observed profiles unless the gravitational potential is appreciably softened in nature.

We have studied the mass distribution in the lensing galaxy 2237+0305 using constraints from both gravitational lensing and photometric and spectroscopic observations. We find that with sufficient dynamical information we can constrain the mass-to-light ratio of the luminous components and determine the dynamical contribution of the disk. In addition, future observations should allow us to place constraints on the shape of the dark matter halo and its inner slope.

We determine the mass distribution of a spiral galaxy, 2237+0305 using both gravitational lensing and dynamical constraints. We find that lensing can break the disc-halo degeneracy. 2237+0305 has a sub-maximal disc, contributing 57+/-3 per cent of the rotational support at the disc maximum. The disc mass-to-light ratio is 1.1+/-0.2 in the I-band and the bulge, 2.9+/-0.5. The dark matter halo, modelled as a softened isothermal sphere, has a large core radius (13.4+/-0.4 kpc, equivalent to 1.4r_d) to high accuracy for the best-fitting solution. The image positions are reasonably well fitted, but require further rotation information to obtain a unique solution.