We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 yr. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS), with the characteristic Hellings-Downs (HD) spatial correlation. Our analysis is performed with two independent pipelines, ENTERPRISE, and TEMPONEST+FORTYTWO, which produce consistent results. A search for a CRS with simultaneous estimation of its spatial correlations yields spectral properties compatible with theoretical GWB predictions, but does not result in the required measurement of the HD correlation, as required for GWB detection. Further Bayesian model comparison between different types of CRSs, including a GWB, finds the most favoured model to be the common uncorrelated red noise described by a power law with $A = 5.13_{-2.73}^{+4.20} \times 10^{-15}$ and $\gamma = 3.78_{-0.59}^{+0.69}$ (95 per cent credible regions). Fixing the spectral index to γ = 13/3 as expected from the GWB by circular, inspiralling supermassive black hole binaries results in an amplitude of $A =2.95_{-0.72}^{+0.89} \times 10^{-15}$. We implement three different models, BAYESEPHEM, LINIMOSS, and EPHEMGP, to address possible Solar system ephemeris (SSE) systematics and conclude that our results may only marginally depend on these effects. This work builds on the methods and models from the studies on the EPTA DR1. We show that under the same analysis framework the results remain consistent after the data set extension.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form hc = A(f/1 yr−1)α, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times 10^{-15}$ and spectral index α = −0.5 ± 0.5, where the uncertainties represqent 95% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of α = −2/3, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times 10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.

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.

Radio pulsars provide some of the most important constraints for our understanding of matter at supranuclear densities. So far, these constraints are mostly given by precision mass measurements of neutron stars (NS). By combining single measurements of the two most massive pulsars, J0348+0432 and J0740+6620, the resulting lower limit of 1.98 M⊙ (99% confidence) of the maximum NS mass, excludes a large number of equations of state (EOSs). Further EOS constraints, complementary to other methods, are likely to come from the measurement of the moment of inertia (MOI) of binary pulsars in relativistic orbits. The Double Pulsar, PSR J0737−3039A/B, is the most promising system for the first measurement of the MOI via pulsar timing. Reviewing this method, based in particular on the first MeerKAT observations of the Double Pulsar, we provide well-founded projections into the future by simulating timing observations with MeerKAT and the SKA. For the first time, we account for the spin-down mass loss in the analysis. Our results suggest that an MOI measurement with 11% accuracy (68% confidence) is possible by 2030. If by 2030 the EOS is sufficiently well known, however, we find that the Double Pulsar will allow for a 7% test of Lense-Thirring precession, or alternatively provide a ∼3σ-measurement of the next-to-leading order gravitational wave damping in GR. Finally, we demonstrate that potential new discoveries of double NS systems with orbital periods shorter than that of the Double Pulsar promise significant improvements in these measurements and the constraints on NS matter.