We present observations and models of the kinematics and the distribution of the neutral hydrogen (HI) in the isolated dwarf irregular galaxy, Wolf-Lundmark-Melotte (WLM). We observed WLM with the Green Bank Telescope (GBT) and as part of the MeerKAT Early Science Programme, where 16 dishes were available. The HI disc of WLM extends out to a major axis diameter of 30 arcmin (8.5 kpc), and a minor axis diameter of 20 arcmin (5.6 kpc) as measured by the GBT. We use the MeerKAT data to model WLM using the TiRiFiC software suite, allowing us to fit different tilted-ring models and select the one that best matches the observation. Our final best-fitting model is a flat disc with a vertical thickness, a constant inclination and dispersion, and a radially-varying surface brightness with harmonic distortions. To simulate bar-like motions, we include second-order harmonic distortions in velocity in the tangential and the vertical directions. We present a model with only circular motions included and a model with non-circular motions. The latter describes the data better. Overall, the models reproduce the global distribution and the kinematics of the gas, except for some faint emission at the 2-sigma level. We model the mass distribution of WLM with a pseudo-isothermal (ISO) and a Navarro-Frenk-White (NFW) dark matter halo models. The NFW and the ISO models fit the derived rotation curves within the formal errors, but with the ISO model giving better reduced chi-square values. The mass distribution in WLM is dominated by dark matter at all radii.

MeerKATHI is the current development name for a radio-interferometric data reduction pipeline, assembled by an international collaboration. We create a publicly available end-to-end continuum- and line imaging pipeline for MeerKAT and other radio telescopes. We implement advanced techniques that are suitable for producing high-dynamic-range continuum images and spectroscopic data cubes. Using containerization, our pipeline is platform-independent. Furthermore, we are applying a standardized approach for using a number of different of advanced software suites, partly developed within our group. We aim to use distributed computing approaches throughout our pipeline to enable the user to reduce larger data sets like those provided by radio telescopes such as MeerKAT. The pipeline also delivers a set of imaging quality metrics that give the user the opportunity to efficiently assess the data quality.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

Accurate instrumental bandpass corrections are essential for the reliable interpretation of spectral lines from targeted and survey-mode observations with radio interferometers. Bandpass correction is typically performed by comparing measurements of a strong calibrator source to an assumed model, typically an isolated point source. The wide field-of-view and high sensitivity of modern interferometers means that additional sources are often detected in observations of calibrators. This can introduce errors into bandpass corrections and subsequently the target data if not properly accounted for. Focusing on the standard calibrator PKS B1934-638, we perform simulations to asses this effect by constructing a wide-field sky model. The cases of ASKAP (0.7-1.9 GHz), MeerKAT (UHF: 0.58-1.05 GHz; L-band: 0.87-1.67 GHz) and Band 2 (0.95-1.76 GHz) of SKA-MID are examined. The use of a central point source model during bandpass calibration is found to impart amplitude errors into spectra measured by the precursor instruments at the ~0.2-0.5% level dropping to ~0.01% in the case of SKA-MID. This manifests itself as ripples in the source spectrum, the behaviour of which is coupled to the distribution of the array baselines, the solution interval, the primary beam size, the hour-angle of the calibration scan, as well as the weights used when imaging the target. Calibration pipelines should routinely employ complete field models for standard calibrators to remove this potentially destructive contaminant from the data, a recommendation we validate by comparing our simulation results to a MeerKAT scan of PKS B1934-638, calibrated with and without our expanded sky model.

We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.

We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

The Galactic Centre contains a supermassive black hole with a mass of 4 million suns within an environment that differs markedly from that of the Galactic disk. While the black hole is essentially quiescent in the broader context of active galactic nuclei, X-ray observations have provided evidence for energetic outbursts from its surroundings. Also, while the levels of star formation in the Galactic Centre have been approximately constant over the last few hundred Myr, there is evidence of elevated short-duration bursts, strongly influenced by interaction of the black hole with the enhanced gas density present within the ring-like Central Molecular Zone at Galactic longitude |l| < 0.7 degrees and latitude |b| < 0.2 degrees. The inner 200 pc region is characterized by large amounts of warm molecular gas, a high cosmic ray ionization rate, unusual gas chemistry, enhanced synchrotron emission, and a multitude of radio-emitting magnetised filaments, the origin of which has not been established. Here we report radio imaging that reveals bipolar bubbles spanning 1 degree x 3 degrees (140 parsecs x 430 parsecs), extending above and below the Galactic plane and apparently associated with the Galactic Centre. The structure is edge-brightened and bounded, with symmetry implying creation by an energetic event in the Galactic Centre. We estimate the age of the bubbles to be a few million years, with a total energy of 7 x 10^52 ergs. We postulate that the progenitor event was a major contributor to the increased cosmic-ray density in the Galactic Centre, and is in turn the principal source of the relativistic particles required to power the synchrotron emission of the radio filaments within and in the vicinity of the bubble cavities.

We present MeerKAT observations of neutral hydrogen gas (HI) in the nearby merger remnant NGC 1316 (Fornax A), the brightest member of a galaxy group which is falling into the Fornax cluster. We find HI on a variety of scales, from the galaxy centre to its large-scale environment. For the first time we detect HI at large radii (70 - 150 kpc in projection), mostly distributed on two long tails associated with the galaxy. Gas in the tails dominates the HI mass of NGC 1316: 7e+8 Msun -- 14 times more than in previous observations. The total HI mass is comparable to the amount of neutral gas found inside the stellar body, mostly in molecular form. The HI tails are associated with faint optical tidal features thought to be the remnant of a galaxy merger occurred a few billion years ago. They demonstrate that the merger was gas-rich. During the merger, tidal forces pulled some gas and stars out to large radii, where we now detect them in the form of optical tails and, thanks to our new data, HI tails; while torques caused the remaining gas to flow towards the centre of the remnant, where it was converted into molecular gas and fuelled the starburst revealed by the galaxy's stellar populations. Several of the observed properties of NGC 1316 can be reproduced by a ~ 10:1 merger between a dominant, gas-poor early-type galaxy and a smaller, gas-rich spiral occurred 1 - 3 Gyr ago, likely followed by subsequent accretion of satellite galaxies.

We present VLA HI observations of JO206, a prototypical ram-pressure stripped galaxy in the GASP sample. This massive galaxy (M$_{\ast} =$ 8.5 $\times$ 10$^{10}$ M$_{\odot}$) is located at a redshift of $z =$ 0.0513, near the centre of the low-mass galaxy cluster, IIZw108 ($\sigma \sim575$ km/s). JO206 is characterised by a long tail ($\geq$90 kpc) of ionised gas stripped away by ram-pressure. We find a similarly long HI tail in the same direction as the ionised gas tail and measure a total HI mass of $3.2 \times 10^{9}$ M$_{\odot}$. This is about half the expected HI mass given the stellar mass and surface density of JO206. A total of $1.8 \times 10^{9}$ M$_{\odot}$ (60%) of the detected HI is in the gas stripped tail. An analysis of the star formation rate shows that the galaxy is forming more stars compared to galaxies with the same stellar and HI mass. On average we find a HI gas depletion time of $\sim$0.5 Gyr which is about four times shorter than that of "normal" spiral galaxies. We performed a spatially resolved analysis of the relation between star formation rate density and gas density in the disc and tail of the galaxy at the resolution of our HI data. The star formation efficiency of the disc is about 10 times higher than that of the tail at fixed HI surface densities. Both the inner and outer parts of JO206 show an enhanced star formation compared to regions of similar HI surface density in field galaxies. The enhanced star formation is due to ram-pressure stripping during the galaxy's first infall into the cluster.

After a decade of design and construction, South Africa's SKA-MID precursor MeerKAT has begun its science operations. To make full use of the widefield capability of the array, it is imperative that we have an accurate model of the primary beam of its antennas. We have used an available L-band full-polarization astro-holographic observation and electromagnetic simulation to create sparse representations of the beam using principal components and Zernike polynomials. The spectral behaviour of the spatial coefficients has been modelled using discrete cosine transform. We have provided the Zernike-based model over a diameter of 10 degrees in an associated software tool that can be useful for direction dependent calibration and imaging. The model is more accurate for the diagonal elements of the beam Jones matrix and at lower frequencies. As we get more accurate beam measurements and simulations in the future, especially for the cross-polarization patterns, our pipeline can be used to create more accurate sparse representations of MeerKAT beam.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.

The new generation of radio interferometers is characterized by high sensitivity, wide fields of view and large fractional bandwidth. To synthesize the deepest images enabled by the high dynamic range of these instruments requires us to take into account the direction-dependent Jones matrices, while estimating the spectral properties of the sky in the imaging and deconvolution algorithms. In this paper we discuss and implement a wide-band wide-field spectral deconvolution framework (DDFacet) based on image plane faceting, that takes into account generic direction-dependent effects. Specifically, we present a wide-field co-planar faceting scheme, and discuss the various effects that need to be taken into account to solve for the deconvolution problem (image plane normalization, position-dependent PSF, etc). We discuss two wide-band spectral deconvolution algorithms based on hybrid matching pursuit and sub-space optimisation respectively. A few interesting technical features incorporated in our imager are discussed, including baseline dependent averaging, which has the effect of improving computing efficiency. The version of DDFacet presented here can account for any externally defined Jones matrices and/or beam patterns.

Almost two decades of observations of radio emission in galaxy clusters have proven the existence of relativistic particles and magnetic fields that generate extended synchrotron emission in the form of radio halos. In the current scenario, radio halos are generated through re--acceleration of relativistic electrons by turbulence generated by cluster mergers. Although this theoretical framework has received increasingly supporting observational evidence over the last ten years, observations of statistically complete samples are needed in order to fundamentally test model predictions. In this paper we briefly review our 7--element Karoo Radio Telescope observations of a sample of nearby clusters aimed to test the predictions of the turbulent re--acceleration model in small systems ($M_{500} > 4 \times 10^{14}$ M$_{\odot}$). We conclude by presenting two galaxy cluster surveys to be carried out with MeerKAT in order to provide crucial test of models of radio halo formation in nearby ($z < 0.1$) and high redshift ($z > 0.4$) systems respectively.