We present details of the Automated Radio Telescope Imaging Pipeline (ARTIP) and results of a sensitive blind search for HI and OH absorbers at $z<0.4$ and $z<0.7$, respectively. ARTIP is written in Python 3.6, extensively uses the Common Astronomy Software Application (CASA) tools and tasks, and is designed to enable the geographically-distributed MeerKAT Absorption Line Survey (MALS) team to collaboratively process large volumes of radio interferometric data. We apply it to the first MALS dataset obtained using the 64-dish MeerKAT radio telescope and 32K channel mode of the correlator. With merely 40 minutes on target, we present the most sensitive spectrum of PKS1830-211 ever obtained and characterize the known HI ($z=0.19$) and OH ($z=0.89$) absorbers. We further demonstrate ARTIP's capabilities to handle realistic observing scenarios by applying it to a sample of 72 bright radio sources observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) to blindly search for HI and OH absorbers. We estimate the numbers of HI and OH absorbers per unit redshift to be $n_{21}(z\sim0.18)<$0.14 and $n_{\rm OH}(z\sim0.40)<$0.12, respectively, and constrain the cold gas covering factor of galaxies at large impact parameters (50 kpc $<\rho<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $\Delta z\sim$13 for HI with column density$>5.4\times10^{19}$ cm$^{-2}$, the survey has probed only the outskirts of star-forming galaxies at $\rho>30$ kpc. MALS with the expected $\Delta z\sim10^{3-4}$ will overcome this limitation and provide stringent constraints on the cold gas fraction of galaxies in diverse environments over $0<z<1.5$.

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.

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 early observations of 20 galaxies in the MHONGOOSE survey galaxies using KAT-7, the GBT, and MeerKAT. We present the best calibrators for five of the MHONGOOSE galaxies observed with the KAT-7, and search for signs of gas accretion in the GBT data, down to $3\sigma$ column density levels of $\rm 2.2\times10^{18}\, cm^{-2}$ over a $\rm 20\, km\, s^{-1}$ linewidth, but identify none. Using the KAT-7 and MeerKAT data, we have derived rotation curves and mass models for NGC 3621 and NGC 7424 out to an unprecedented extent. As a precursor to the SKA, the MeerKAT telescope combines both a high spatial resolution and a large field of view, necessary to map the extended neutral hydrogen in local galaxies. The mass models of the two galaxies were constructed for both the Dark Matter (DM) models (the pseudo-isothermal model and the Navarro-Frenk-White model) and MOND. Overall, we find that the DM models provide a better fit than MOND to the galaxies' rotation curves. Furthermore, the pseudo-isothermal model is found to be the most consistent with the observations.

HI observations of the Local Group dwarf irregular galaxies Sextans A and B, obtained with the Karoo Array Telescope (KAT-7) are presented. The KAT-7 wide field of view and excellent surface brightness sensitivity allows us to verify the true HI extent of the galaxies. We derive HI extends of 30$^{\prime}$ and 20$^{\prime}$ and total HI fluxes of 181 $\pm$ 2.0 Jy.km.s$^{-1}$ and 105 $\pm$ 1.4 Jy.km.s$^{-1}$ for Sextans A and B respectively. This result shows clearly the overestimate of the HI extent and total flux of 54$^{\prime}$ and 264 Jy.km.s$^{-1}$ reported for Sextans A using the Effelsberg observations. Tilted ring models allow us to derive the rotation curves (RCs) of Sextans A and B out to 550$^{\prime \prime}$ ($\sim$ 3.5 kpc) and 650$^{\prime \prime}$ ($\sim$ 4 kpc) respectively. The RCs of the two galaxies are seen to decline in the outer parts. The dark matter distribution in Sextans A is better described by the pseudo-isothermal halo model when a M/L ratio of 0.2 is used. For Sextans B, the mass model fits are not as good but again an isothermal sphere with a M/L of 0.2 represents best the data. Using the MOdified Newtonian Dynamics (MOND), better fits are obtained when the constant a$_{0}$ is allowed to vary. The critical densities for gravitational instabilities are calculated using the Toomre-$Q$ and cloud-growth based on shear criterion. We find that in regions of star formation, the cloud growth criterion based on shear explains better the star formation in both Sextans A and B.

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 present high sensitivity HI observations of NGC 6822, obtained with the Karoo Array Telescope (KAT-7). We study the kinematics, the mass distribution, and the star formation thresholds. The KAT-7 short baselines and low system temperature make it sensitive to large-scale, low surface brightness emission. The observations detected $\sim$ 23$\%$ more flux than previous ATCA observations. We fit a tilted ring model to the HI velocity field to derive the rotation curve (RC). The KAT-7 observations allow the measurement of the rotation curve of NGC 6822 out to 5.8 kpc, $\sim$ 1 kpc further than existing measurements. NGC 6822 is seen to be dark matter dominated at all radii. The observationally motivated pseudo-isothermal dark matter (DM) halo model reproduces well the observed RC while the Navarro Frank-White DM model gives a poor fit to the data. We find the best fit mass to light ratio (M/L) of 0.12 $\pm$ 0.01 which is consistent with the literature. The Modified Newtonian Dynamics (MOND) gives a poor fit to our data. We derive the star formation threshold in NGC 6822 using the HI and H$\alpha$ data. The critical gas densities were calculated for gravitational instabilities using the Toomre-Q criterion and the cloud-growth criterion. We found that in regions of star formation, the cloud-growth criterion explains star formation better than the Toomre-Q criterion. This shows that the local shear rate could be a key player in cloud formation for irregular galaxies such as NGC 6822.

Extending deep observations of the neutral atomic hydrogen (HI) to the environment around galaxy groups can reveal a complex history of group interactions which is invisible to studies that focus on the stellar component. Hickson Compact Group 44 (HCG 44) is a nearby example and we have combined HI data from the Karoo Array Telescope, Westerbork Synthesis Radio Telescope, and Arecibo Legacy Fast ALFA survey, in order to achieve high column density sensitivity (N_HI < 2x10^18 cm^-2) to the neutral gas over a large field-of-view beyond the compact group itself. We find the giant HI tail north of HCG 44 contains 1.1x10^9 M_Sun of gas and extends 450 kpc from the compact group: twice as much mass and 33% further than previously detected. However, the additional gas is still unable to account for the known HI deficiency of HCG 44. The tail likely formed through a strong tidal interaction and HI clouds in the tail have survived for 1 Gyr or more after being stripped. This has important implications for understanding the survival of neutral clouds in the intragroup and circumgroup medium, and we discuss their survival in the context of simulations of cold gas in hot halos. HCG 44 is one of a growing number of galaxy groups found to have more extended HI in the intragroup and circumgroup medium than previously measured. Our results provide constraints for simulations on the properties of galaxy group halos, and reveal a glimpse of what will be seen by future powerful HI telescopes and surveys.

The construction of the KAT-7 array in the Karoo region of the Northern Cape in South Africa was intended primarily as an engineering prototype for technologies and techniques applicable to the MeerKAT telescope. This paper looks at the main engineering and scien- tific highlights from this effort, and discusses their applicability to both MeerKAT and other next-generation radio telescopes. In particular we found that the composite dish surface works well, but it becomes complicated to fabricate for a dish lacking circular symmetry; the Stir- ling cycle cryogenic system with ion pump to achieve vacuum works but demands much higher maintenance than an equivalent Gifford-McMahon cycle system; the ROACH (Recon- figurable Open Architecture Computing Hardware)-based correlator with SPEAD (Stream- ing Protocol for Exchanging Astronomical Data) protocol data transfer works very well and KATCP (Karoo Array Telescope Control Protocol) control protocol has proven very flexible and convenient. KAT-7 has also been used for scientific observations where it has a niche in mapping low surface-brightness continuum sources, some extended HI halos and OH masers in star-forming regions. It can also be used to monitor continuum source variability, observe pulsars, and make VLBI observations

The Antlia Cluster is a nearby, dynamically young structure, and its proximity provides a valuable opportunity for detailed study of galaxy and group accretion onto clusters. We present a deep HI mosaic completed as part of spectral line commissioning of the Karoo Array Telescope (KAT-7), and identify infrared counterparts from the WISE extended source catalog to study neutral atomic gas content and star formation within the cluster. We detect 37 cluster members out to a radius of ~0.9 Mpc with M_HI > 5x10^7 M_Sun. Of these, 35 are new HI detections, 27 do not have previous spectroscopic redshift measurements, and one is the Compton thick Seyfert II, NGC 3281, which we detect in HI absorption. The HI galaxies lie beyond the X-ray emitting region 200 kpc from the cluster center and have experienced ram pressure stripping out to at least 600 kpc. At larger radii, they are distributed asymmetrically suggesting accretion from surrounding filaments. Combining HI with optical redshifts, we perform a detailed dynamical analysis of the internal substructure, identify large infalling groups, and present the first compilation of the large scale distribution of HI, and star forming galaxies within the cluster. We find that elliptical galaxy NGC 3268 is at the center of the oldest substructure and argue that NGC 3258 and its companion population are more recent arrivals. Through the presence of HI and on-going star formation, we rank substructures with respect to their relative time since accretion onto Antlia.

HI observations of the Magellanic-type spiral NGC 3109, obtained with the seven dish Karoo Array Telescope (KAT-7), are used to analyze its mass distribution. Our results are compared to what is obtained using VLA data. KAT-7 is the precursor of the SKA pathfinder MeerKAT, which is under construction. The short baselines and low system temperature of the telescope make it sensitive to large scale low surface brightness emission. The new observations with KAT-7 allow the measurement of the rotation curve of NGC 3109 out to 32', doubling the angular extent of existing measurements. A total HI mass of 4.6 x 10^8 Msol is derived, 40% more than what was detected by the VLA observations. The observationally motivated pseudo-isothermal dark matter (DM) halo model can reproduce very well the observed rotation curve but the cosmologically motivated NFW DM model gives a much poorer fit to the data. While having a more accurate gas distribution has reduced the discrepancy between the observed RC and the MOdified Newtonian Dynamics (MOND) models, this is done at the expense of having to use unrealistic mass-to-light ratios for the stellar disk and/or very large values for the MOND universal constant a0. Different distances or HI contents cannot reconcile MOND with the observed kinematics, in view of the small errors on those two quantities. As for many slowly rotating gas-rich galaxies studied recently, the present result for NGC 3109 continues to pose a serious challenge to the MOND theory.

The relationship between the clustering of dark matter and that of luminous matter is often described using the bias parameter. Here, we provide a new method to probe the bias of intermediate to high-redshift radio continuum sources for which no redshift information is available. We matched radio sources from the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey data to their optical counterparts in the Sloan Digital Sky Survey (SDSS) to obtain photometric redshifts for the matched radio sources. We then use the publicly available semi-empirical simulation of extragalactic radio continuum sources (S3) to infer the redshift distribution for all FIRST sources and estimate the redshift distribution of unmatched sources by subtracting the matched distribution from the distribution of all sources. We infer that the majority of unmatched sources are at higher redshifts than the optically matched sources and demonstrate how the angular scales of the angular two-point correlation function can be used to probe different redshift ranges. We compare the angular clustering of radio sources with that expected for dark matter and estimate the bias of different samples.

We investigate the clustering of HI-selected galaxies in the ALFALFA survey and compare results with those obtained for HIPASS. Measurements of the angular correlation function and the inferred 3D-clustering are compared with results from direct spatial-correlation measurements. We are able to measure clustering on smaller angular scales and for galaxies with lower HI masses than was previously possible. We calculate the expected clustering of dark matter using the redshift distributions of HIPASS and ALFALFA and show that the ALFALFA sample is somewhat more anti-biased with respect to dark matter than the HIPASS sample.

We use the Millennium simulation to probe the correlation between cluster velocities and their shapes and the consequences for measurements of the kinetic Sunyaev-Zeldovich (kSZ) effect. Halos are generally prolate ellipsoids with orientations that are correlated with those of nearby halos. We measure the mean streaming velocities of halos along the lines that separate them, demonstrating that the peculiar velocities and the long axes of halos tend to be somewhat aligned, especially for the most massive halos. Since the kSZ effect is proportional to the line-of-sight velocity and the optical depth of the cluster, the alignment results in a strong enhancement of the kSZ signature in clusters moving along the line of sight. This effect has not been taken into account in many analyses of kSZ signatures.