A new high time resolution observing mode for the Murchison Widefield Array (MWA) is described, enabling full polarimetric observations with up to 30.72 MHz of bandwidth and a time resolution of ~0.8 $\mu$s. This mode makes use of a polyphase synthesis filter to "undo" the polyphase analysis filter stage of the standard MWA's Voltage Capture System (VCS) observing mode. Sources of potential error in the reconstruction of the high time resolution data are identified and quantified, with the S/N loss induced by the back-to-back system not exceeding -0.65 dB for typical noise-dominated samples. The system is further verified by observing three pulsars with known structure on microsecond timescales.

We present an analysis of Murchison Widefield Array radio telescope data from $\omega$ Cen, possibly a stripped dwarf spheroidal galaxy core captured by our Galaxy. Recent interpretations of Fermi-LAT $\gamma$-ray data by Brown {\it et al.} (2019) and Reynoso-Cordova {\it et al.} (2019) suggest that $\omega$ Cen may contain significant Dark Matter. We utilise their best-fit Dark Matter annihilation models, and an estimate of the magnetic field strength in $\omega$ Cen, to calculate the expected radio synchrotron signal from annihilation, and show that one can usefully rule out significant parts of the magnetic field - diffusion coefficient plane using the current observational limits. Improvement by a factor of 10-100 on these limits could constrain the models even more tightly.

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 carousel model of pulsar emission attributes the phenomenon of subpulse drifting to a set of discrete sparks located very near the stellar surface rotating around the magnetic axis. Here, we investigate the subpulse drifting behaviour of PSR B0031-07 in the context of the carousel model. We show that B0031-07's three drift modes (A, B, and C) can be understood in terms of a single carousel rotation rate if the number of sparks is allowed to change by an integral number, and where the different drift rates are due to (first-order) aliasing effects. This also results in harmonically-related values for P 3 (the time it takes a subpulse to reappear at the same pulse phase), which we confirm for B0031-07. A representative solution has [n_A, n_B, n_C] = [15, 14, 13] sparks and a carousel rotation period of P_4 = 16.4 P_1. We also investigate the frequency dependence of B0031-07's subpulse behaviour. We extend the carousel model to include the dual effects of aberration and retardation, including the time it takes the information about the surface spark configuration to travel from the surface up to the emission point. Assuming these effects dominate at B0031-07's emission heights, we derive conservative emission height differences of $\lesssim 2000$ km for mode A and $\lesssim 1000$ km for modes B and C as seen between 185 MHz and 610 MHz. This new method of measuring emission heights is independent of others that involve average profile components or the polarisation position angle curve, and thus provides a potentially strong test of the carousel model.

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.

In this paper we demonstrate improved techniques to extend coherent processing intervals for passive radar processing, with the Murchison Widefield Array. Specifically, we apply a two stage linear range and Doppler migration compensation by utilising Keystone Formatting and a recent dechirping method. These methods are used to further demonstrate the potential for the surveillance of space with the Murchison Widefield Array using passive radar, by detecting objects orders of magnitude smaller than previous work. This paper also demonstrates how the linear Doppler migration methods can be extended to higher order compensation to further increase potential processing intervals.

The Murchison Widefield Array is a low-frequency Square Kilometre Array precursor located at the Murchison Radio-astronomy Observatory in Western Australia. Primarily designed as an imaging telescope, but with a flexible signal path, the capabilities of this telescope have recently been extended to include off-line incoherent and tied-array beam formation using recorded antenna voltages. This has provided the capability for high-time and frequency resolution observations, including a pulsar science program. This paper describes the algorithms and pipeline that we have developed to form the tied array beam products from the summation of calibrated signals of the antenna elements, and presents example polarimetric profiles for PSRs J0437-4715 and J1900-2600 at 185 MHz.

Polarimetric studies of pulsars at low radio frequencies provide important observational insights into the pulsar emission mechanism and beam models, and probe the properties of the magneto-ionic interstellar medium (ISM). Aperture arrays are the main form of next-generation low-frequency telescopes, including the Murchison Widefield Array (MWA). These require a distinctly different approach to data processing (e.g. calibration and beamforming) compared to traditional dish antennas. As the second paper of this series, we present a verification of the MWA's pulsar polarimetry capability, using two bright southern pulsars, PSRs J0742-2822 and J1752-2806. Our observations simultaneously cover multiple frequencies (76-313 MHz) and were taken at multiple zenith angles during a single night for each pulsar. We show that the MWA can be reliably calibrated for zenith angles < 45 degree and frequencies < 270 MHz. We present the polarimetric profiles for PSRs J0742-2822 and J1752-2806 at frequencies lower than 300 MHz for the first time, along with an analysis of the linear polarisation degree and pulse profile evolution with frequency. For PSR J0742-2822, the measured degree of linear polarisation shows a rapid decrease at low frequencies, in contrast with the generally expected trend, which can be attributed to depolarisation effects from small-scale, turbulent, magneto-ionic ISM components. This effect has not been widely explored for pulsars in general, and will be further investigated in future work.

The phenomenon of subpulse drifting offers unique insights into the emission geometry of pulsars, and is commonly interpreted in terms of a rotating carousel of "spark" events near the stellar surface. We develop a detailed geometric model for the emission columns above a carousel of sparks that is entirely calculated in the observer's inertial frame, and which is consistent with the well-understood rotational effects of aberration and retardation. We explore the observational consequences of the model, including (1) the appearance of the reconstructed beam pattern via the cartographic transform and (2) the morphology of drift bands and how they might evolve as a function of frequency. The model, which is implemented in the software package PSRGEOM, is applicable to a wide range of viewing geometries, and we illustrate its implications using PSRs B0809+74 and B2034+19 as examples. Some specific predictions are made with respect to the difference between subpulse evolution and microstructure evolution, which provides a way to further test our model.

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.

The Murchison Widefield Array (MWA), and its recently-developed Voltage Capture System (VCS), facilitates extending the low-frequency range of pulsar observations at high-time and -frequency resolution in the Southern Hemisphere, providing further information about pulsars and the ISM. We present the results of an initial time-resolved census of known pulsars using the MWA. To significantly reduce the processing load, we incoherently sum the detected powers from the 128 MWA tiles, which yields ~10% of the attainable sensitivity of the coherent sum. This preserves the large field-of-view (~450 deg2 at 185 MHz), allowing multiple pulsars to be observed simultaneously. We developed a WIde-field Pulsar Pipeline (WIPP) that processes the data from each observation and automatically folds every known pulsar located within the beam. We have detected 50 pulsars to date, 6 of which are millisecond pulsars. This is consistent with our expectation, given the telescope sensitivity and the sky coverage of the processed data (~17,000 deg2). For ten pulsars, we present the lowest-frequency detections published. For a subset of the pulsars, we present multi-frequency pulse profiles by combining our data with published profiles from other telescopes. Since the MWA is a low-frequency precursor to the Square Kilometre Array (SKA), we use our census results to forecast that a survey using Phase 1 of SKA-Low (SKA1-Low) can potentially detect around 9400 pulsars.

The phenomenon of subpulse drifting may hold the key to understanding the pulsar emission mechanism. Here, we report on new observations of PSR J0034-0721 (B0031-07), carried out with the Murchison Widefield Array at 185 MHz. We observe three distinct drift modes whose "vertical" drift band separations ($P_3$) and relative abundances are consistent with previous studies at similar and higher frequencies. The driftbands, however, are observed to change their slopes over the course of individual drift modes, which can be interpreted as a continuously changing drift rate. The implied acceleration of the intrinsic carousel rotation cannot easily be explained by plasma models based on ExB drift. Furthermore, we find that methods of classifying the drift modes by means of $P_3$ measurements can sometimes produce erroneous identifications in the presence of a changing drift rate. The "horizontal" separation between driftbands ($P_2$) is found to be larger at later rotation phases within the pulse window, which is inconsistent with the established effects of retardation, aberration, and the motion of the visible point. Longer observations spanning at least ~10,000 pulses are required to determine how the carousel rotation parameters change from one drift sequence to the next.

Low-frequency observations of pulsars provide a powerful means for probing the microstructure in the turbulent interstellar medium (ISM). Here we report on high-resolution dynamic spectral analysis of our observations of the timing-array millisecond pulsar J0437-4715 with the Murchison Widefield Array (MWA), enabled by our recently commissioned tied-array beam processing pipeline for voltage data recorded from the high time resolution mode of the MWA. A secondary spectral analysis reveals faint parabolic arcs, akin to those seen in high-frequency observations of pulsars with the Green Bank and Arecibo telescopes. Data from Parkes observations at a higher frequency of 732 MHz reveal a similar parabolic feature, with a curvature that scales approximately as the square of the observing wavelength ($\lambda^2$) to the MWA's frequency of 192 MHz. Our analysis suggests that scattering toward PSR J0437-4715 predominantly arises from a compact region about 115 pc from the Earth, which matches well with the expected location of the edge of the Local Bubble that envelopes the local Solar neighborhood. As well as demonstrating new and improved pulsar science capabilities of the MWA, our analysis underscores the potential of low-frequency pulsar observations for gaining valuable insights into the local ISM and for characterizing the ISM toward timing-array pulsars.

We report on observations of giant pulses from the Crab pulsar performed simultaneously with the Parkes radio telescope and the incoherent combination of the Murchison Widefield Array (MWA) antenna tiles. The observations were performed over a duration of approximately one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center frequency of 193 MHz with 15 MHz bandwidth at the MWA. Our analysis has led to the detection of 55 giant pulses at the MWA and 2075 at Parkes above a threshold of 3.5$\sigma$ and 6.5$\sigma$ respectively. We detected 51$\%$ of the MWA giant pulses at the Parkes radio telescope, with spectral indices in the range of $-3.6>\alpha> -4.9$ ($S_{\rm \nu} \propto \nu^\alpha$). We present a Monte Carlo analysis supporting the conjecture that the giant pulse emission in the Crab is intrinsically broadband, the less than $100\%$ correlation being due to the relative sensitivities of the two instruments and the width of the spectral index distribution. Our observations are consistent with the hypothesis that the spectral index of giant pulses is drawn from normal distribution of standard deviation 0.6, but with a mean that displays an evolution with frequency from -3.00 at 1382 MHz, to -2.85 at 192 MHz.