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 magnetic field plays an important role in the evolution of the Galaxy, but its small-scale behaviour is still poorly known. It is also unknown whether it permeates the halo of the Galaxy or not. By using observations of pulsars in the halo globular cluster 47 Tucanae, we probed the Galactic magnetic field at arcsecond scales for the first time and discovered an unexpected large gradient in the component of the magnetic field parallel to the line of sight. This gradient is aligned with a direction perpendicular to the Galactic disk and could be explained by magnetic fields amplified to some 60 {\mu}G within the globular cluster. This scenario supports the existence of a magnetized outflow that extends from the Galactic disk to the halo and interacts with the studied globular cluster.

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.

Theoretical models suggest that intermediate mass black holes (IMBHs) may form and reside in the centers of globular clusters. IMBHs are still elusive to observations, but the accelerations of pulsars may bring along a unique fingerprint of their presence. In this work, we focus on the pulsars in the globular cluster M62. Using the new distance of M62 obtained from Gaia observations, we find that the measured pulsars' accelerations suggest a central excess of mass in the range [1200, 6000]$M_{\odot}$, corresponding to [0.2, 1] percent of the current total mass of the cluster. Our analysis can not unambiguously discriminate between an IMBH or a system of stellar mass dark remnants of comparable total mass.

Globular clusters may host intermediate mass black holes (IMBHs) at their centres. Here we propose a new method for their identification using millisecond pulsars (MSPs) as probes. We show that measuring the first (jerk) and second (jounce) derivatives of the accelerations of an ensemble of MSPs will let us infer the presence of an IMBH in a globular cluster better than measuring the sole accelerations. We test this concept by simulating a set of star clusters with and without a central IMBH to extract the distributions of the stellar jerks and jounces. We then apply this technique to the ensemble of MSPs in the Galactic globular cluster 47 Tucanae. Current timing observations are insufficient to constrain the presence of an IMBH and can only be used to pose upper limits on its mass. But, with few more years of observations it will be possible to test for the presence of a central IMBH with mass smaller than $\sim$ 1000 M$_{\odot}$. We conclude that jerks and jounces help significantly in reducing the upper limit of the mass of IMBHs in Galactic globular clusters.

Despite considerations of mass loss from stellar evolution suggesting otherwise, the content of gas in globular clusters seems poor and hence its measurement very elusive. One way of constraining the presence of ionized gas in a globular cluster is through its dispersive effects on the radiation of the millisecond pulsars included in the cluster. This effect led Freire et al. in 2001 to the first detection of any kind of gas in a globular cluster in the case of 47 Tucanae. By exploiting the results of 12 additional years of timing, as well as the observation of new millisecond pulsars in 47 Tucanae, we revisited this measurement: we first used the entire set of available timing parameters in order to measure the dynamical properties of the cluster and the three-dimensional position of the pulsars. Then we applied and tested various gas distribution models: assuming a constant gas density, we confirmed the detection of ionized gas with a number density of $n= 0.23\pm 0.05$ cm$^{-3}$, larger than the previous determination (at 2$\sigma$ uncertainty). Models predicting a decreasing density or following the stellar distribution density are highly disfavoured. We are also able to investigate the presence of an intermediate mass black hole in the centre of the cluster, showing that is not required by the available data, with an upper limit for the mass at $\sim 4000$ M$_{\odot}$.

The origin of the Nuclear Star Cluster in the centre of our Galaxy is still unknown. One possibility is that it formed after the disruption of stellar clusters that spiralled into the Galactic Centre due to dynamical friction. We trace the formation of the Nuclear Star Cluster around the central black hole, using state-of-the-art N-body simulations, and follow the dynamics of the neutron stars born in the clusters. We then estimate the number of Millisecond Pulsars (MSPs) that are released in the Nuclear Star Cluster, during its formation. The assembly and tidal dismemberment of globular clusters lead to a population of MSPs distributed over a radius of about 20 pc, with a peak near 3 pc. No clustering is found on the sub-parsec scale. We simulate the detectability of this population with future radio telescopes like the MeerKAT radio telescope and SKA1, and find that about of order ten MSPs can be observed over this large volume, with a paucity of MSPs within the central parsec. This helps discriminating this scenario from the in-situ formation model for the Nuclear Star Cluster that would predict an over abundance of MSPs closer to the black hole. We then discuss the potential contribution of our MSP population to the gamma-ray excess at the Galactic Centre.