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.

We present HI observations of the jellyfish galaxy, JO201. This massive galaxy (M$_{\ast} = 3.5 \times 10^{10}$ M$_\odot$) is falling along the line-of-sight towards the centre of a rich cluster (M$_{200} \sim 1.6 \times 10^{15}$ M$_\odot$, $\sigma_{cl} \sim 982$ km/s) at a high velocity $\geq$3363 km/s. Its H$\alpha$ emission shows a $\sim$40 kpc tail confined closely to its stellar disc and a $\sim$100 kpc tail extending further out. We find HI emission coinciding only with the shorter clumpy H$\alpha$ tail. In total, we measure an HI mass of M$_{\rm HI} = 1.65 \times 10^{9}$ M$_\odot$, which is about 60% lower than expected based on its stellar mass and stellar surface density. We compared JO201 to another jellyfish in the GASP sample, JO206 (of similar mass but residing in a 10$\times$ less massive cluster), and find that they are similarly HI-deficient. Of the total HI mass in JO201, about 30% lies outside the galaxy disc in projection. This HI fraction is probably a lower limit since most of the HI is redshifted relative to the stellar disc and could be outside the disc. The global star formation rate (SFR) analysis of JO201 suggests that its observed SFR would be expected if it had 10$\times$ its current HI mass. The disc is the main contributor of the high star formation efficiency at a given HI gas density for both galaxies, but their tails also show higher star formation efficiencies compared to the outer regions of field galaxies. Generally, we find that JO201 and JO206 are similar based on their HI content, stellar mass and star formation rate. This finding is unexpected considering their different environments. A toy model comparing the ram pressure of the ICM versus the restoring forces of these galaxies suggests that the ram pressure strength exerted on them could be comparable if we consider their 3D orbital velocities and radial distances relative to the clusters.

Modern radio telescopes are routinely reaching depths where normal starforming galaxies are the dominant observed population. Realising the potential of radio as a tracer of star formation and black hole activity over cosmic time involves achieving such depths over representative volumes, with radio forming part of a larger multiwavelength campaign. In pursuit of this we used the Karl G. Jansky Very Large Array (VLA) to image $\sim$5 deg$^{2}$ of the VIDEO/XMM-LSS extragalactic deep field at 1--2 GHz. We achieve a median depth of 16 $\mu$Jy beam$^{-1}$ with an angular resolution of 4.5\arcsec. Comparisons with existing radio observations of XMM-LSS showcase the improved survey speed of the upgraded VLA: we cover 2.5 times the area and increase the depth by $\sim$20\% in 40\% of the time. Direction-dependent calibration and wide-field imaging were required to suppress the error patterns from off-axis sources of even modest brightness. We derive a catalogue containing 5,762 sources from the final mosaic. Sub-band imaging provides in-band spectral indices for 3,458 (60\%) sources, with the average spectrum becoming flatter than the canonical synchrotron slope below 1 mJy. Positional and flux-density accuracy of the observations, and the differential source counts are in excellent agreement with those of existing measurements. A public release of the images and catalogue accompanies this article.

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 present MeerKAT 1000 MHz and 1400 MHz observations of a bright radio galaxy in the southern hemisphere, ESO~137-006. The galaxy lies at the centre of the massive and merging Norma galaxy cluster. The MeerKAT continuum images (rms ~0.02 mJy/beam at ~10" resolution) reveal new features that have never been seen in a radio galaxy before: collimated synchrotron threads of yet unknown origin, which link the extended and bent radio lobes of ESO~137-006. The most prominent of these threads stretches in projection for about 80 kpc and is about 1 kpc in width. The radio spectrum of the threads is steep, with a spectral index of up to $\alpha\simeq 2$ between 1000 MHz and 1400 MHz.

We present new observations of Fornax A taken at 1 GHz with the MeerKAT telescope and at 6 GHz with the Sardinia Radio Telescope (SRT). The sensitive (noise ~16 micro-Jy beam$^{-1}$), high resolution ( < 10'') MeerKAT images show that the lobes of Fornax A have a double-shell morphology, where dense filaments are embedded in a diffuse and extended cocoon. We study the spectral properties of these components by combining the MeerKAT and SRT observations with archival data between 84 MHz and 217 GHz. For the first time, we show that multiple episodes of nuclear activity must have formed the extended radio lobes. The modelling of the radio spectrum suggests that the last episode of injection of relativistic particles into the lobes started ~ 24 Myr ago and stopped approximately 12 Myr ago. More recently (~ 3 Myr ago), a less powerful and short ( < 1 Myr) phase of nuclear activity generated the central jets. Currently, the core may be in a new active phase. It appears that Fornax A is rapidly flickering. The dense environment in which Fornax A lives has lead to a complex recent merger history for this galaxy, including mergers spanning a range of gas contents and mass ratios, as shown by the analysis of the galaxy's stellar- and cold-gas phases. This complex recent history may be the cause of the rapid, recurrent nuclear activity of Fornax A.

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.

Modern interferometric imaging relies on advanced calibration that incorporates direction-dependent effects. Their increasing number of antennas (e.g. in LOFAR, VLA, MeerKAT/SKA) and sensitivity are often tempered with the accuracy of their calibration. Beam accuracy drives particularly the capability for high dynamic range imaging (HDR - contrast > 1:$10^6$). The Radio Interferometric Measurement Equation (RIME) proposes a refined calibration framework for wide field of views (i.e. beyond the primary lobe and first null) using beam models. We have used holography data taken on 12 antennas of the Very Large Array (VLA) with two different approaches: a `data-driven' representation derived from Principal Component Analysis (PCA) and a projection on the Zernike polynomials. We determined sparse representations of the beam to encode its spatial and spectral variations. For each approach, we compressed the spatial and spectral distribution of coefficients using low-rank approximations. The spectral behaviour was encoded with a Discrete Cosine Transform (DCT). We compared our modelling to that of the Cassbeam software which provides a parametric model of the antenna and its radiated field. We present comparisons of the beam reconstruction fidelity vs. `compressibility'. We found that the PCA method provides the most accurate model. In the case of VLA antennas, we discuss the frequency ripple over L-band which is associated with a standing wave between antenna reflectors. The results are a series of coefficients that can easily be used `on-the-fly' in calibration pipelines to generate accurate beams at low computing costs.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

Kliko is a Docker-based container specification for running one or multiple related compute jobs. The key concepts of Kliko are the encapsulation of data processing software into a container and the formalization of the input, output and task parameters. By formalizing the parameters, the software is represented as abstract building blocks with a uniform and consistent interface. The main advantage is enhanced scriptability and empowering pipeline composition. Formalization is realized by bundling a container with a Kliko file, which describes the IO and task parameters. This Kliko container can then be opened and run by a Kliko runner. The Kliko runner will parse the Kliko definition and gather the values for these parameters, for example by requesting user input or retrieving pre-defined values from disk. Parameters can be various primitive types, for example: float, int or the path to a file. This paper will also discuss the implementation of a support library named Kliko which can be used to create Kliko containers, parse Kliko definitions, chain Kliko containers in workflows using a workflow manager library such as Luigi. The Kliko library can be used inside the container to interact with the Kliko runner. Finally, to illustrate the applicability of the Kliko definition, this paper will discuss two reference implementations based on the Kliko library: RODRIGUES, a web-based Kliko container scheduler, and output visualizer specifically for astronomical data, and VerMeerKAT, a multi-container workflow data reduction pipeline which is being used as a prototype pipeline for the commissioning of the MeerKAT radio telescope.

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.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

The MeerKAT telescope will be one of the most sensitive radio arrays in the pre-SKA era. Here we discuss a low-frequency SZ-selected cluster survey with MeerKAT, the MeerKAT Extended Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey. The primary goal of this survey is to detect faint signatures of diffuse cluster emission, specifically radio halos and relics. SZ-selected cluster samples offer a homogeneous, mass-limited set of targets out to higher redshift than X-ray samples. MeerKAT is sensitive enough to detect diffuse radio emission at the faint levels expected in low-mass and high-redshift clusters, thereby enabling radio halo and relic formation theories to be tested with a larger statistical sample over a significantly expanded phase space. Complementary multiwavelength follow-up observations will provide a more complete picture of any clusters found to host diffuse emission, thereby enhancing the scientific return of the MERGHERS survey.

The high sensitivity of the new generation of radio telescopes such as the Square Kilometre Array (SKA) will allow cosmological weak lensing measurements at radio wavelengths that are competitive with optical surveys. We present an adaptation to radio data of lensfit, a method for galaxy shape measurement originally developed and used for optical weak lensing surveys. This likelihood method uses an analytical galaxy model and makes a Bayesian marginalisation of the likelihood over uninteresting parameters. It has the feature of working directly in the visibility domain, which is the natural approach to adopt with radio interferometer data, avoiding systematics introduced by the imaging process. As a proof of concept, we provide results for visibility simulations of individual galaxies with flux density S >= 10 muJy at the phase centre of the proposed SKA1-MID baseline configuration, adopting 12 frequency channels in the band 950 - 1190 MHz. Weak lensing shear measurements from a population of galaxies with realistic flux and scalelength distributions are obtained after natural gridding of the raw visibilities. Shear measurements are expected to be affected by "noise bias": we estimate the bias in the method as a function of signal-to-noise ratio (SNR). We obtain additive and multiplicative bias values that are comparable to SKA1 requirements for SNR > 18 and SNR > 30, respectively. The multiplicative bias for SNR > 10 is comparable to that found in ground-based optical surveys such as CFHTLenS, and we anticipate that similar shear measurement calibration strategies to those used for optical surveys may be used to good effect in the analysis of SKA radio interferometer data.

Observationally, weak lensing has been served so far by optical surveys due to the much larger number densities of background galaxies achieved, which is typically by two to three orders of magnitude compared to radio. However, the high sensitivity of the new generation of radio telescopes such as the Square Kilometre Array (SKA) will provide a density of detected galaxies that is comparable to that found at optical wavelengths, and with significant source shape measurements to make large area radio surveys competitive for weak lensing studies. This will lead weak lensing to become one of the primary science drivers in radio surveys too, with the advantage that they will access the largest scales in the Universe going beyond optical surveys, like LSST and Euclid, in terms of redshifts that are probed. RadioLensfit is an adaptation to radio data of "lensfit", a model-fitting approach for galaxy shear measurement, originally developed for optical weak lensing surveys. Its key advantage is working directly in the visibility domain, which is the natural approach to adopt with radio data, avoiding systematics due to the imaging process. We present results on galaxy shear measurements, including investigation of sensitivity to instrumental parameters such as the visibilities gridding size, based on simulations of individual galaxy visibilities performed by using SKA1-MID baseline configuration. We get an amplitude of the shear bias in the method comparable with SKA1 requirements for a population of galaxies with realistic flux and scalelength distributions estimated from the VLA SWIRE catalog.

Weak gravitational lensing measurements are traditionally made at optical wavelengths where many highly resolved galaxy images are readily available. However, the Square Kilometre Array (SKA) holds great promise for this type of measurement at radio wavelengths owing to its greatly increased sensitivity and resolution over typical radio surveys. The key to successful weak lensing experiments is in measuring the shapes of detected sources to high accuracy. In this document we describe a simulation pipeline designed to simulate radio images of the quality required for weak lensing, and will be typical of SKA observations. We provide as input, images with realistic galaxy shapes which are then simulated to produce images as they would have been observed with a given radio interferometer. We exploit this pipeline to investigate various stages of a weak lensing experiment in order to better understand the effects that may impact shape measurement. We first show how the proposed SKA1-Mid array configurations perform when we compare the (known) input and output ellipticities. We then investigate how making small changes to these array configurations impact on this input-outut ellipticity comparison. We also demonstrate how alternative configurations for SKA1-Mid that are smaller in extent, and with a faster survey speeds produce similar performance to those originally proposed. We then show how a notional SKA configuration performs in the same shape measurement challenge. Finally, we describe ongoing efforts to utilise our simulation pipeline to address questions relating to how applicable current (mostly originating from optical data) shape measurement techniques are to future radio surveys. As an alternative to such image plane techniques, we lastly discuss a shape measurement technique based on the shapelets formalism that reconstructs the source shapes directly from the visibility data.

Morphologically classifying radio sources in continuum images with the SKA has the potential to address some of the key questions in cosmology and galaxy evolution. In particular, we may use different classes of radio sources as independent tracers of the dark-matter density field, and thus overcome cosmic variance in measuring large-scale structure, while on the galaxy evolution side we could measure the mechanical feedback from FRII and FRI jets. This work makes use of a \texttt{MeqTrees}-based simulations framework to forecast the ability of the SKA to recover true source morphologies at high redshifts. A suite of high resolution images containing realistic continuum source distributions with different morphologies (FRI, FRII, starburst galaxies) is fed through an SKA Phase 1 simulator, then analysed to determine the sensitivity limits at which the morphologies can still be distinguished. We also explore how changing the antenna distribution affects these results.

Galaxy clusters are known to host a variety of extended radio sources: tailed radio galaxies whose shape is modelled by the interaction with the intra-cluster medium (ICM); radio bubbles filling cavities in the ICM distribution and rising buoyantly through the thermal gas; diffuse giant radio sources ("halos" and "relics") revealing the presence of relativistic electrons and magnetic fields in the intra-cluster volume. It is currently the subject of an active debate how the non-thermal components that we observe at radio wavelengths affect the physical properties of the ICM and depend on the dynamical state of galaxy clusters. In this work we start our SKA1 feasibility study of the "radio cluster zoo" through simulations of a typical radio-loud cluster, hosting several bright tailed radio galaxies and a diffuse radio halo. Realistic simulations of SKA1 observations are obtained through the MeqTrees software. A new deconvolution algorithm, based on sparse representations and optimised for the detection of faint diffuse astronomical sources, is tested and compared to the classical CLEAN method.