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 ( ${\sim}2.8\,\mbox{K Jy}^{-1}$ ) low-system temperature ( ${\sim}18\,\mbox{K at }20\,\mbox{cm}$ ) radio array that currently operates at 580–1 670 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 $\mbox{J}0737{-}3039\mbox{A}$ , pulse profiles from 34 millisecond pulsars (MSPs) from a single 2.5-h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR $\mbox{J}0540{-}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 MSPs confirm that MeerKAT delivers exceptional timing. PSR $\mbox{J}2241{-}5236$ exhibits a jitter limit of $<4\,\mbox{ns h}^{-1}$ whilst timing of PSR $\mbox{J}1909{-}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 1 000 pulsars per day and the future deployment of S-band (1 750–3 500 MHz) receivers will further enhance its capabilities.

The initial results from timing observations of PSR J1141–6545, a relativistic pulsar-white dwarf binary system, are presented. Predictions from the timing baseline hint at the most stringent test of gravity by an asymmetric binary yet. The timing precision has been hindered by the dramatic variations of the pulse profile due to geodetic precession, a pulsar glitch and red timing noise. Methods to overcome such timing irregularities are briefly presented along with preliminary results from the test of the General Theory of Relativity (GR) from this pulsar.

The radial velocity of a pulsar induces the Doppler effect on its intrinsic spin properties. In particular, it can generate a contribution to the frequency second derivative. We estimated this effect for each of the International Pulsar Timing Array pulsars. We also assessed the possibility of measuring the frequency second derivative in the observational data.

I describe the concept of a pulsar timing array and give broad overview of the construction of a pulsar timing array, methods for high-precision timing and noise characterization, and algorithms for gravitational wave detection and source characterization. I then provide an overview of worldwide pulsar timing programs and the scale and sensitivity of the pulsar timing array efforts, with particular attention to the International Pulsar Timing Array (IPTA). I discuss the most recent results from pulsar timing arrays, emphasizing the gravitational wave detection efforts in particular. Finally, I describe the anticipated future growth in participants, telescopes, pulsars, and sensitivity of the IPTA, highlighting the transformational advances that it will enable over the next decade.

We perform a Bayesian analysis of pulsar-timing residuals from the NANOGrav pulsar-timing array to search for a specific form of stochastic narrow-band signal produced by oscillating gravitational potential (Gravitational Potential Background) in the Galactic halo. Such oscillations arise in models of warm dark matter composed of an ultralight massive scalar field (m ≃ 10−23 eV). The propagation of an electromagnetic signal from a pulsar through the time-dependent spacetime will leave an imprint in the pulsar timing, much like a gravitational wave. From the physical point of view, this is the classical Sachs-Wolfe effect. A distinctive feature of the pulsar-timing residuals due to GBP produced by a variable scalar field is that the amplitude of the TOA residuals should be independent of the pulsar location in the sky. In the monochromatic approximation, the stringent upper limit (95% C.L.) on the variable gravitational potential amplitude is found to be (Ψc<1.14 × 10−15), corresponding to the characteristic strain hc = 2$\sqrt{3}$Ψc < 4 × 10−15 at f=1.75 × 10−8 Hz. In the narrow-band approximation, the upper limit of this background energy density is ΩGPB < 1.27 × 10−9 at f=1.75 × 10−8 Hz. These limits are an order of magnitude higher than the expected signal amplitude assuming all Galactic dark matter is made of such scalar particles. The applied analysis of the pulsar-timing residuals can be used to search for any narrow-band stochastic signals with different correlation properties. As a by-product, parameters of the red noise present in four NANOGrav pulsars were found.

Simultaneous timing of several pulsars distributed over the sky, so called Pulsar Timing Array (PTA), is used for a variety of metrological and astronomical applications. Three examples of PTA application are presented: link between celestial reference frames, ensemble pulsar time scale and detection of gravitational waves.