Search

The Large Magellanic Cloud (LMC) presents a unique environment for pulsar population studies due to its distinct star formation characteristics and proximity to the Milky Way. As part of the TRAPUM (TRAnsients and PUlsars with MeerKAT) Large Survey Project, we are using the core array of the MeerKAT radio telescope (MeerKAT) to conduct a targeted search of the LMC for radio pulsars at L-band frequencies, 856-1712$\,$MHz. The excellent sensitivity of MeerKAT, coupled with a 2-hour integration time, makes the survey 3 times more sensitive than previous LMC radio pulsar surveys. We report the results from the initial four survey pointings which has resulted in the discovery of seven new radio pulsars, increasing the LMC radio pulsar population by 30 per cent. The pulse periods of these new pulsars range from 278 to 1690$\,$ms, and the highest dispersion measure is 254.20$\,$pc$\,$cm$^{-3}$. We searched for, but did not find any significant pulsed radio emission in a beam centred on the SN$\,$1987A remnant, establishing an upper limit of 6.3$\,{\mu}$Jy on its minimum flux density at 1400$\,$MHz., Comment: 12 pages, 4 figures, 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We report new radio timing solutions from a three-year observing campaign conducted with the MeerKAT and Murriyang telescopes for nine Small Magellanic Cloud pulsars, increasing the number of characterised rotation-powered extragalactic pulsars by 40 per cent. We can infer from our determined parameters that the pulsars are seemingly all isolated, that six are ordinary pulsars, and that three of the recent MeerKAT discoveries have a young characteristic age of under 100 kyr and have undergone a spin-up glitch. Two of the sources, PSRs J0040$-$7337 and J0048$-$7317, are energetic young pulsars with spin-down luminosities of the order of 10$^{36}$ erg s$^{-1}$. They both experienced a large glitch, with a change in frequency of about 30 $\mu$Hz, and a frequency derivative change of order $-10^{-14}$ Hz s$^{-1}$. These glitches, the inferred glitch rate, and the properties of these pulsars (including potentially high inter-glitch braking indices) suggest these neutron stars might be Vela-like repeating glitchers and should be closely monitored in the future. The position and energetics of PSR J0048$-$7317 confirm it is powering a new Pulsar Wind Nebula (PWN) detected as a radio continuum source; and similarly the association of PSR J0040$-$7337 with the PWN of Supernova Remnant (SNR) DEM S5 (for which we present a new Chandra image) is strengthened. Finally, PSR J0040$-$7335 is also contained within the same SNR but is a chance superposition. It has also been seen to glitch with a change of frequency of $10^{-2}$ $\mu$Hz. This work more than doubles the characterised population of SMC radio pulsars., Comment: 20 pages, 13 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The Sculptor Galaxy (NGC 253), located in the Southern Hemisphere, far off the Galactic Plane, has a relatively high star-formation rate of about 7 M$_{\odot}$ yr$^{-1}$ and hosts a young and bright stellar population, including several super star clusters and supernova remnants. It is also the first galaxy, apart from the Milky Way Galaxy to be associated with two giant magnetar flares. As such, it is a potential host of pulsars and/or fast radio bursts in the nearby Universe. The instantaneous sensitivity and multibeam sky coverage offered by MeerKAT therefore make it a favourable target. We searched for pulsars, radio-emitting magnetars, and fast radio bursts in NGC 253 as part of the TRAPUM large survey project with MeerKAT. We did not find any pulsars during a four-hour observation, and derive a flux density limit of 4.4 $\mu$Jy at 1400 MHz, limiting the pseudo-luminosity of the brightest putative pulsar in this galaxy to 54 Jy kpc$^2$. Assuming universality of pulsar populations between galaxies, we estimate that detecting a pulsar as bright as this limit requires NGC 253 to contain a pulsar population of $\gtrsim$20 000. We also did not detect any single pulses and our single pulse search flux density limit is 62 mJy at 1284\,MHz. Our search is sensitive enough to have detected any fast radio bursts and radio emission similar to the brighter pulses seen from the magnetar SGR J1935+2154 if they had occurred during our observation., Comment: 6 pages, 3 figures, accepted for publication in MNRAS

The eclipses seen in the radio emission of some pulsars can be invaluable to study the properties of the material from the companion stripped away by the pulsar. We present a study of six consecutive eclipses of PSR J0024-7204O in the globular cluster 47 Tucanae as seen by the MeerKAT radio telescope in the UHF (544-1088 MHz) band. A high scintillation state boosted the signal during one of the orbits and allowed a detailed study of the eclipse properties. We measure significant dispersion measure (DM) variations and detect strong scattering that seems to be the dominating mechanism of the eclipses at these frequencies. A complete drop in the linear polarization together with a small increase in the rotation measure suggests the presence of a magnetic field of $\sim 2$ mG. The study of multiple eclipses allowed us to measure difference in the lengths of the eclipses and DM differences of $\sim 0.01$ pc cm$^{-3}$ in consecutive orbits. One orbit in particular shows a delay in recovery of the linear polarization and a visible delay in the arrival of the pulses caused by a stronger scattering event. We suggest that these are caused by a higher variance of density fluctuations during the event., Comment: 14 pages, 8 figures, 9 tables. Accepted for publication on Monthly Notices of the Royal Astronomical Society

Long-period radio transients are an emerging class of extreme astrophysical events of which only three are known. These objects emit highly polarised, coherent pulses of typically a few tens of seconds duration and minutes to hour-long periods. While magnetic white dwarfs and magnetars, either isolated or in binary systems, have been invoked to explain these objects, a consensus has not emerged. Here we report on the discovery of ASKAP J193505.1+214841.0 (henceforth ASKAPJ1935+2148) with a period of 53.8 minutes exhibiting three distinct emission states - a bright pulse state with highly linearly polarised pulses with widths of 10-50 seconds; a weak pulse state which is about 26 times fainter than the bright state with highly circularly polarised pulses of widths of approximately 370 milliseconds; and a quiescent or quenched state with no pulses. The first two states have been observed to progressively evolve over the course of 8 months with the quenched state interspersed between them suggesting physical changes in the region producing the emission. A constraint on the radius of the source for the observed period rules out a magnetic white dwarf origin. Unlike other long-period sources, ASKAPJ1935+2148 is the first to exhibit drastic variations in emission modes reminiscent of neutron stars. However, its radio properties challenge our current understanding of neutron star emission and evolution., Comment: Published in Nature Astronomy

Millisecond pulsars are subject to accelerations in globular clusters that manifest themselves in both the first and second spin period time derivatives, and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 years of pulsar timing observations of five millisecond radio pulsars in the core of the core-collapse globular cluster NGC\,6752 with the Parkes (Murriyang) and MeerKAT radio telescopes, that have allowed us to measure the proper motions, positions and first and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of non-luminous mass of at least $2.56\times10^3M_\odot$ is resent in the core of NGC\,6752. On the other hand, our studies highly disfavour the presence of an intermediate mass black hole at the center of the cluster, with a mass equal to or greater than $\sim3000M_\odot$., Comment: Accepted for publication in ApJ

The sensitivity of the MeerKAT radio interferometer is an opportunity to probe deeper into the population of rare and faint extragalactic pulsars. The TRAPUM (TRAnsients and PUlsars with MeerKAT) collaboration has conducted a radio-domain search for accelerated pulsars and transients in the Small Magellanic Cloud (SMC). This partially targeted survey, performed at L-band (856-1712 MHz) with the core array of the MeerKAT telescope in 2-h integrations, is twice as sensitive as the latest SMC radio pulsar survey. We report the discovery of seven new SMC pulsars, doubling this galaxy's radio pulsar population and increasing the total extragalactic population by nearly a quarter. We also carried out a search for accelerated millisecond pulsars in the SMC Globular Cluster NGC 121 using the full array of MeerKAT. This improved the previous upper limit on pulsed radio emission from this cluster by a factor of six. Our discoveries reveal the first radio pulsar-PWN systems in the SMC, with only one such system previously known outside our galaxy (the "Crab pulsar twin" in the Large Magellanic Cloud, PSR J0540$-$6919). We associate the 59 ms pulsar discovery PSR J0040$-$7337, now the fastest spinning radio pulsar in the SMC, with the bow-shock Pulsar Wind Nebula (PWN) of Supernova Remnant DEM S5. We also present a new young pulsar with a 79 ms period, PSR J0048$-$7317, in a PWN recently discovered in a MeerKAT radio continuum image. Using the multi-beam capability of MeerKAT, we localised our pulsar discoveries, and two previous Murriyang discoveries, to a positional uncertainty of a few arcseconds., Comment: 32 pages, 14 figures, 10 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish MeerKAT radio telescope, we use an interferometric beamforming technique to tile the potential pulsar locations with coherent beams which we search for radio pulsations, above a signal-to-noise of 9, down to an average flux density upper limit of 30 $\mu$Jy. This limit is target-dependent due to the contribution of the sky and nebula to the system temperature. Coherent beams are arranged to overlap at their 50 per cent power radius, so the sensitivity to pulsars is not degraded by more than this amount, though realistically averages around 65 per cent if every location in the beam is considered. We report the discovery of two new pulsars; PSR J1831$-$0941 is an adolescent pulsar likely to be the plerionic engine of the candidate PWN G20.0+0.0, and PSR J1818$-$1502 appears to be an old and faint pulsar that we serendipitously discovered near the centre of a SNR already hosting a compact central object. The survey holds importance for better understanding of neutron star birth rates and the energetics of young pulsars., Comment: 17 pages, 6 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We report here on the timing of 597 pulsars over the last four years with the MeerKAT telescope. We provide Times-of-Arrival, pulsar ephemeris files and per-epoch measurements of the flux density, dispersion measure (DM) and rotation measure (RM) for each pulsar. In addition we use a Gaussian process to model the timing residuals to measure the spin frequency derivative at each epoch. We also report the detection of 11 glitches in 9 individual pulsars. We find significant DM and RM variations in 87 and 76 pulsars respectively. We find that the DM variations scale approximately linearly with DM, which is broadly in agreement with models of the ionised interstellar medium. The observed RM variations seem largely independent of DM, which may suggest that the RM variations are dominated by variations in the interstellar magnetic field on the line of sight, rather than varying electron density. We also find that normal pulsars have around 5 times greater amplitude of DM variability compared to millisecond pulsars, and surmise that this is due to the known difference in their velocity distributions., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. 12 pages 7 figures

We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for nine discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748$-$2446ao which is an eccentric ($e = 0.32$) wide-orbit (orbital period $P_{\rm b} = 57.55$ d) system. We were able to measure the rate of advance of periastron ($\dot{\omega}$) for this system allowing us to determine a total mass of $3.17 \pm \, 0.02\, \rm M_{\odot}$. With a minimum companion mass ($M_{\rm c}$) of $\sim 0.8\, \rm M_{\odot}$, PSR J1748$-$2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar ($P = 2.27$ ms) and the longest orbital period measured for any known DNS system. PSR J1748$-$2446ap has the second highest eccentricity for any recycled pulsar ($e \sim 0.905$) and for this system we can measure the total mass ($1.997 \pm 0.006\, \rm M_{\odot}$) and also estimate the individual pulsar and companion masses. PSR J1748$-$2446ar is an eclipsing redback (minimum $M_{\rm c} \sim 0.34\, \rm M_{\odot}$) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies., Comment: 23 pages, 11 figures, 5 tables, published in A&A

Using MeerKAT, we have discovered three new millisecond pulsars (MSPs) in the bulge globular cluster M62: M62H, M62I, and M62J. All three are in binary systems, which means all ten known pulsars in the cluster are in binaries. M62H has a planetary-mass companion with a median mass $M_{\rm c,med} \sim 3$ M$_{\rm J}$ and a mean density of $\rho \sim 11$ g cm$^{-3}$. M62I has an orbital period of 0.51 days and a $M_{\rm c,med} \sim 0.15$ M$_{\odot}$. Neither of these low-mass systems exhibit eclipses. M62J has only been detected in the two UHF band (816 MHz) observations with a flux density $S_{816} = 0.08$ mJy. The non-detection in the L-band (1284 MHz) indicates it has a relatively steep spectrum ($\beta < -3.1$). We also present 23-yr-long timing solutions obtained using data from the Parkes "Murriyang", Effelsberg and MeerKAT telescopes for the six previously known pulsars. For all these pulsars, we measured the second spin-period derivatives and the rate of change of orbital period caused by the gravitational field of the cluster, and their proper motions. From these measurements, we conclude that the pulsars' maximum accelerations are consistent with the maximum cluster acceleration assuming a core-collapsed mass distribution. Studies of the eclipses of the redback M62B and the black widow M62E at four and two different frequency bands, respectively, reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies. The presence of only binary MSPs in this cluster challenges models which suggest that the MSP population of core-collapsed clusters should be dominated by isolated MSPs., Comment: 22 pages, 8 figures, 6 tables. Accepted for publication in MNRAS

Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838$-$2527, J0955$-$3947 and J2333$-$5526, confirming their redback nature. PSR J0838$-$2827 remained undetected for two years after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-year timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333$-$5526 may have a neutron star mass in excess of 2 M$_{\odot}$., Comment: 19 pages, 7 figures, accepted for publication in MNRAS

We present new radio continuum images and a source catalogue from the MeerKAT survey in the direction of the Small Magellanic Cloud (SMC). The observations, at a central frequency of 1.3 GHz across a bandwidth of 0.8 GHz, encompass a field of view ~7 x 7 degrees and result in images with resolution of 8 arcsec. The median broad-band Stokes I image Root Mean Squared noise value is ~11 microJy/beam. The catalogue produced from these images contains 108,330 point sources and 517 compact extended sources. We also describe a UHF (544-1088 MHz) single pointing observation. We report the detection of a new confirmed Supernova Remnant (SNR) (MCSNR J0100-7211) with an X-ray magnetar at its centre and 10 new SNR candidates. This is in addition to the detection of 21 previously confirmed SNRs and two previously noted SNR candidates. Our new SNR candidates have typical surface brightness an order of magnitude below those previously known, and on the whole they are larger. The high sensitivity of the MeerKAT survey also enabled us to detect the bright end of the SMC Planetary Nebulae (PNe) sample - point-like radio emission is associated with 38 of 102 optically known PNe, of which 19 are new detections. Lastly, we present the detection of three foreground radio stars amidst 11 circularly polarised sources, and a few examples of morphologically interesting background radio galaxies from which the radio ring galaxy ESO 029-G034 may represent a new type of radio object., Comment: 31 pages, 27 figures

We report on a radiopolarimetric observation of the Saturn-Jupiter Great Conjunction of 2020 using the MeerKAT L-band system, initially carried out for science verification purposes, which yielded a serendipitous discovery of a pulsar. The radiation belts of Jupiter are very bright and time variable: coupled with the sensitivity of MeerKAT, this necessitated development of dynamic imaging techniques, reported on in this work. We present a deep radio "movie" revealing Jupiter's rotating magnetosphere, a radio detection of Callisto, and numerous background radio galaxies. We also detect a bright radio transient in close vicinity to Saturn, lasting approximately 45 minutes. Follow-up deep imaging observations confirmed this as a faint compact variable radio source, and yielded detections of pulsed emission by the commensal MeerTRAP search engine, establishing the object's nature as a radio emitting neutron star, designated PSR J2009-2026. A further observation combining deep imaging with the PTUSE pulsar backend measured detailed dynamic spectra for the object. While qualitatively consistent with scintillation, the magnitude of the magnification events and the characteristic timescales are odd. We are tentatively designating this object a pulsar with anomalous refraction recurring on odd timescales (PARROT). As part of this investigation, we present a pipeline for detection of variable sources in imaging data, with dynamic spectra and lightcurves as the products, and compare dynamic spectra obtained from visibility data with those yielded by PTUSE. We discuss MeerKAT's capabilities and prospects for detecting more of such transients and variables., Comment: 21 pages, 17 figures, submitted to MNRAS

We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251\deg $\le l \le$ 358\deg and 2\deg $\le l \le$ 61\deg at $|b| \le 1.5\deg $). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $\mu$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance., Comment: Accepted for publication in MNRAS. The data release is live and links can be found in the Data Availability Statement in the paper

We present full Stokes MeerKAT L band (856--1712\,MHz) observations of \chg{36} high latitude supernova remnants. Sensitive, high dynamic range images show a wealth of structure. G15.1$-$1.6 appears to be an HII region rather than an SNR. G30.7$-$2.0 consists of three background extragalactic sources which appear to form an arc when imaged with much lower resolution. At least half of the remnants in the sample contain ``blowouts'', or ``ears'' showing these to be a common feature. Analysis of the polarimetric data reveals details of the magnetic field structure in the emitting regions of the remnants as well as magnetized thermal plasma in front of polarized emission. The chance alignment of G327.6+14.6 with a background AGN with very extended polarized jets allows testing for the presence of Faraday effects in the interior of the remnant. Scant evidence of Faraday rotating material is found in the interior of this remnant., Comment: 41 pages, 58 figures

The double pulsar system, PSR J0737$-$3039A/B, consists of two neutron stars bound together in a highly relativistic orbit that is viewed nearly edge-on from the Earth. This alignment results in brief radio eclipses of the fast-rotating pulsar A when it passes behind the toroidal magnetosphere of the slow-rotating pulsar B. The morphology of these eclipses is strongly dependent on the geometric orientation and rotation phase of pulsar B, and their time-evolution can be used to constrain the geodetic precession rate of the pulsar. We demonstrate a Bayesian inference framework for modelling eclipse light-curves obtained with MeerKAT between 2019-2023. Using a hierarchical inference approach, we obtained a precession rate of $\Omega_{\rm SO}^{\rm B} = {5.16^{\circ}}^{+0.32^{\circ}}_{-0.34^{\circ}}$ yr$^{-1}$ for pulsar B, consistent with predictions from General Relativity to a relative uncertainty of 6.5%. This updated measurement provides a 6.1% test of relativistic spin-orbit coupling in the strong-field regime. We show that a simultaneous fit to all of our observed eclipses can in principle return a $\sim$1.5% test of spin-orbit coupling. However, systematic effects introduced by the current geometric orientation of pulsar B along with inconsistencies between the observed and predicted eclipse light curves result in difficult to quantify uncertainties. Assuming the validity of General Relativity, we definitively show that the spin-axis of pulsar B is misaligned from the total angular momentum vector by $40.6^{\circ} \pm 0.1^{\circ}$ and that the orbit of the system is inclined by approximately $90.5^{\circ}$ from the direction of our line of sight. Our measured geometry for pulsar B suggests the largely empty emission cone contains an elongated horseshoe shaped beam centered on the magnetic axis, and that it may not be re-detected as a radio pulsar until early-2035., Comment: Abridged abstract. 13 pages, 9 figures and 2 tables

We present the results of observations aimed at discovering and studying pulsars in the core-collapsed globular cluster (GC) NGC 6522 performed by the MeerTIME and TRAPUM Large Survey Project with the MeerKAT telescope. We have discovered two new isolated pulsars bringing the total number of known pulsars in the cluster to six. PSR J1803$-$3002E is a mildly recycled pulsar with spin period of 17.9 ms while pulsar PSR J1803$-$3002F is a slow pulsar with spin period of 148.1 ms. The presence of isolated and slow pulsars is expected in NGC 6522 and confirms the predictions of previous theories for clusters at this stage in evolution. We further present a tentative timing solution for the millisecond pulsar (MSP) PSR J1803$-$3002C combining older observations taken with the Parkes 64m radio telescope, Murriyang. This solution implies a relatively small characteristic age of the pulsar in contrast with the old age of the GC. The presence of a slow pulsar and an apparently young MSP, both rare in GCs, suggests that their formation might be linked to the evolutionary stage of the cluster., Comment: 8 pages, 5 figues, 3 tables. Accepted for publication in Astronomy and Astrophysics

The MMGPS-L is the most sensitive pulsar survey in the Southern Hemisphere. We present a follow-up study of one of these new discoveries, PSR J1208-5936, a 28.71-ms recycled pulsar in a double neutron star system with an orbital period of Pb=0.632 days and an eccentricity of e=0.348. Through timing of almost one year of observations, we detected the relativistic advance of periastron (0.918(1) deg/yr), resulting in a total system mass of Mt=2.586(5) Mo. We also achieved low-significance constraints on the amplitude of the Einstein delay and Shapiro delay, in turn yielding constraints on the pulsar mass (Mp=1.26(+0.13/-0.25) Mo), the companion mass (Mc=1.32(+0.25/-0.13) Mo, and the inclination angle (i=57(12) degrees). This system is highly eccentric compared to other Galactic field double neutron stars with similar periods, possibly hinting at a larger-than-usual supernova kick during the formation of the second-born neutron star. The binary will merge within 7.2(2) Gyr due to the emission of gravitational waves. With the improved sensitivity of the MMGPS-L, we updated the Milky Way neutron star merger rate to be 25(+19/-9) Myr$^{-1}$ within 90% credible intervals, which is lower than previous studies based on known Galactic binaries owing to the lack of further detections despite the highly sensitive nature of the survey. This implies a local cosmic neutron star merger rate of 293(+222/-103} Gpc/yr, consistent with LIGO and Virgo O3 observations. With this, we predict the observation of 10(+8/-4) neutron star merger events during the LIGO-Virgo-KAGRA O4 run. We predict the uncertainties on the component masses and the inclination angle will be reduced to 5x10$^{-3}$ Mo and 0.4 degrees after two decades of timing, and that in at least a decade from now the detection of the shift in Pb and the sky proper motion will serve to make an independent constraint of the distance to the system.

We present a comprehensive study of the rotational and emission properties of PSR J0738$-$4042 using a combination of observations taken by the Deep Space Network, Hartebeesthoek, Parkes (Murriyang) and Molonglo observatories between 1972 and 2023. Our timing of the pulsar is motivated by previously reported profile/spin-down events that occurred in September 2005 and December 2015, which result in an anomalously large braking index of $n = 23300 \pm 1800$. Using a Gaussian process regression framework, we develop continuous models for the evolution of the pulsar spin-down rate ($\dot{\nu}$) and profile shape. We find that the pulse profile variations are similar regardless of radio observing frequency and polarisation. Small-scale differences can be ascribed to changes in the interstellar medium along the line of sight and frequency-dependent changes in magnetospheric radio emission height. No new correlated spin-down or profile events were identified in our extended dataset. However, we found that the disappearance of a bright emission component in the leading edge of archival profiles between 1981-1988 was not associated with a substantial change in $\dot{\nu}$. This marks a notable departure from the previous profile/spin-down events in this pulsar. We discuss the challenges these observations pose for physical models and conclude that interactions between the pulsar and in-falling asteroids or a form of magnetospheric state-switching with a long periodicity are plausible explanations., Comment: 14 pages, 8 figures. Fixed typo on page 5. Accepted for publication in MNRAS

PSR J1933$-$6211 is a 3.5-ms pulsar in a 12.8-d orbit with a white dwarf (WD). Its high proper motion and low dispersion measure result in such significant interstellar scintillation that high signal-to-noise detections require long observing durations or fortuitous timing. We turn to the sensitive MeerKAT telescope and, combined with historic Parkes data, leverage PSR J1933$-$6211's kinematic and relativistic effects to constrain its 3D orbital geometry and the component masses. We obtain precise proper motion and parallax estimates, and measure their effects as secular changes in the Keplerian orbital parameters: a variation in orbital period of $7(1) \times 10^{-13}$ s s$^{-1}$ and a change in projected semi-major axis of $1.60(5) \times 10^{-14}$ s s$^{-1}$. A self-consistent analysis of all kinematic and relativistic effects yields a distance of $1.6^{+0.2}_{-0.3}$ kpc, an orbital inclination, $i = 55(1)$ deg and a longitude of the ascending node, $\Omega = 255^{+8}_{-14}$ deg. The probability densities for $\Omega$ and $i$ and their symmetric counterparts, ($180-i$, $360-\Omega$), are seen to depend on the fiducial orbit used to measure the time of periastron passage. We investigate this unexpected dependence and rule out software-related causes using simulations. Nevertheless, we constrain the pulsar and WD masses to $1.4^{+0.3}_{-0.2}$ M$_\odot$ and $0.43(5)$ M$_\odot$ respectively. These strongly disfavour a helium-dominated WD. The orbital similarities between PSRs J1933$-$6211 and J1614$-$2230 suggest they underwent Case A Roche lobe overflow, an extended evolution while the companion star is still on the Main Sequence. However, with a mass of $\sim 1.4$ M$_\odot$, PSR J1933$-$6211 has not accreted significant matter. This highlights the low accretion efficiency of the spin-up process and suggests that observed neutron star masses are mostly a result of supernova physics., Comment: 16 pages, 7 figures. Abstract shortened to adhere to ArXiv limit

Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena. Performing such a survey using the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. We present here the overview and setup for the 3000 hour Max-Planck-Institut fuer Radioastronomie (MPIfR) MeerKAT Galactic Plane survey (MMGPS). The survey is unique by operating in a commensal mode, addressing key science objectives of the survey including the discovery of new pulsars and transients as well as studies of Galactic magnetism, the interstellar medium and star formation rates. We explain the strategy coupled with the necessary hardware and software infrastructure needed for data reduction in the imaging, spectral and time domains. We have so far discovered 78 new pulsars including 17 confirmed binary systems of which two are potential double neutron star systems. We have also developed an imaging pipeline sensitive to the order of a few tens of micro-Jansky with a spatial resolution of a few arcseconds. Further science operations with an in-house built S-Band receiver operating between 1.7-3.5 GHz are about to commence. Early spectral line commissioning observations conducted at S-Band, targeting transitions of the key molecular gas tracer CH at 3.3 GHz already illustrate the spectroscopic capabilities of this instrument. These results lay a strong foundation for future surveys with telescopes like the Square Kilometre Array (SKA)., Comment: 25 pages, 10 figures, Accepted in MNRAS

We present a sample of well-localised Fast Radio Bursts (FRBs) discovered by the MeerTRAP project at the MeerKAT telescope in South Africa. We discovered the three FRBs in single coherent tied-array beams and localised them to an area of ~1 arcmin$^2$. We investigate their burst properties, scattering, repetition rates, and localisations in a multi-wavelength context. FRB 20201211A shows hints of scatter broadening but is otherwise consistent with instrumental dispersion smearing. For FRB 20210202D, we discovered a faint post-cursor burst separated by ~200 ms, suggesting a distinct burst component or a repeat pulse. We attempt to associate the FRBs with host galaxy candidates. For FRB 20210408H, we tentatively (0.35 - 0.53 probability) identify a compatible host at a redshift ~0.5. Additionally, we analyse the MeerTRAP survey properties, such as the survey coverage, fluence completeness, and their implications for the FRB population. Based on the entire sample of 11 MeerTRAP FRBs discovered by the end of 2021, we estimate the FRB all-sky rates and their scaling with the fluence threshold. The inferred FRB all-sky rates at 1.28 GHz are $8.2_{-4.6}^{+8.0}$ and $2.1_{-1.1}^{+1.8} \times 10^3$ sky$^{-1}$ d$^{-1}$ above 0.66 and 3.44 Jy ms for the coherent and incoherent surveys, respectively. The scaling between the MeerTRAP rates is flatter than at higher fluences at the 1.4-$\sigma$ level. There seems to be a deficit of low-fluence FRBs, suggesting a break or turn-over in the rate versus fluence relation below 2 Jy ms. We speculate on cosmological or progenitor-intrinsic origins. The cumulative source counts within our surveys appear consistent with the Euclidean scaling., Comment: 23 pages, 10 figures, 3 tables. Accepted for publication in MNRAS. Supplementary data are available on Zenodo (https://doi.org/10.5281/zenodo.6047539) and via the Transient Name Server

We report on the subpulse modulation properties of 1198 pulsars using the Thousand-Pulsar-Array programme on MeerKAT. About 35% of the analysed pulsars exhibit drifting subpulses which are more pronounced towards the deathline, consistent with previous studies. We estimate that this common phenomenon is detectable in 60% of the overall pulsar population if high quality data were available for all. This large study reveals the evolution of drifting subpulses across the pulsar population in unprecedented detail. In particular, we find that the modulation period $P_3$ follows a V-shaped evolution with respect to the characteristic age $\tau_c$, such that the smallest $P_3$ values, corresponding to the Nyquist period $P_3>\sim2$, are found at $\tau_c>\sim10^{7.5}$ yr. The V-shaped evolution can be interpreted and reproduced if young pulsars possess aliased fast intrinsic $P_3$, which monotonically increase, ultimately achieving a slow unaliased $P_3$. Enhancement of irregularities in intrinsic subpulse modulation by aliasing in small $\tau_c$ pulsars would explain their observed less well defined $P_3$'s and weaker spectral features. Modelling these results as rotating subbeams, their circulation must slow down as the pulsar evolves. This is the opposite to that expected if circulation is driven by ExB drift. This can be resolved if the observed $P_3$ periodicity is due to a beat between an ExB system and the pulsar period. As a by-product, we identify the correct periods and spin-down rates for 12 pulsars, for which harmonically related values were reported in the literature., Comment: 18 pages, 13 figures in the main part of the paper. This version includes all online materials from the accepted MNRAS article

PSRJ1910-5959A (J1910A) is a binary millisecond pulsar in a 0.837 day circular orbit around a helium white dwarf (HeWD) companion. This pulsar is located 6.3 arcmin away from the centre of the globular cluster NGC6752. Given the large offset, the association of the pulsar to NGC6752 has been debated. We have made use of two decades of archival Parkes 64-m "Murriyang" telescope data and recently carried out observations with the MeerKAT telescope. We obtained Pulse times of arrival using standard data reduction techniques and analysed using Bayesian pulsar timing techniques. We analysed the pulsar's total intensity and polarisation profile, to study the interstellar scattering along the line of sight, and the pulsar's geometry by applying the rotating vector model. We obtain precise measurements of several post-Keplerian parameters: the range $r=0.202(6)T_\odot$ and shape s=0.999823(4) of the Shapiro delay, from which we infer the orbital inclination to be $88.9^{+0.15}_{-0.14}\deg$ and the masses of both the pulsar and the companion to be $1.55(7)M_{\odot}$ and $0.202(6)M_{\odot}$ respectively; a secular change in the orbital period $\dot{P}_{\rm b}=-53^{+7.4}_{-6.0}\times 10^{-15}$\,s\,s$^{-1}$ that proves the association to NGGC6752 and a secular change in the projected semi-major axis of the pulsar $\dot{x}= -40.7^{+7.3}_{-8.2}\times10^{-16}$\,s\,s$^{-1}$ that is likely caused by the spin-orbit interaction from a misaligned HeWD spin, at odds with the likely isolated binary evolution of the system. We also discuss some theoretical models for the structure and evolution of WDs in NS-WD binaries by using J1910A's companion as a test bed. J1910A is a rare system for which several parameters of both the pulsar and the HeWD companion can be accurately measured. As such, it is a test bed to discriminate between alternative models for HeWD structure and cooling., Comment: 18 pages, 12 figures. Abstract shortened and rephrased because of arXiv characters limits. Accepted for publication on A&A on January 9th 2023 Updated version after linguistic revision

The most massive globular cluster in our Galaxy, Omega Centauri, is an interesting target for pulsar searches, because of its multiple stellar populations and the intriguing possibility that it was once the nucleus of a galaxy that was absorbed into the Milky Way. The recent discoveries of pulsars in this globular cluster and their association with known X-ray sources was a hint that, given the large number of known X-ray sources, there is a much larger undiscovered pulsar population. We used the superior sensitivity of the MeerKAT radio telescope to search for pulsars in Omega Centauri. In this paper, we present some of the first results of this survey, including the discovery of 13 new pulsars; the total number of known pulsars in this cluster currently stands at 18. At least half of them are in binary systems and preliminary orbital constraints suggest that most of the binaries have light companions. We also discuss the ratio between isolated and binaries pulsars and how they were formed in this cluster.

More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sources identified as possible gamma-ray pulsar candidates by a Random Forest classification of unassociated sources from the 4FGL catalogue. Each source was observed for 10 minutes on two separate epochs using MeerKAT's L-band receiver (856-1712 MHz), with typical pulsed flux density sensitivities of $\sim$100$\,\mu$Jy. Nine new MSPs were discovered, eight of which are in binary systems, including two eclipsing redbacks and one system, PSR J1526$-$2744, that appears to have a white dwarf companion in an unusually compact 5 hr orbit. We obtained phase-connected timing solutions for two of these MSPs, enabling the detection of gamma-ray pulsations in the Fermi-LAT data. A follow-up search for continuous gravitational waves from PSR J1526$-$2744 in Advanced LIGO data using the resulting Fermi-LAT timing ephemeris yielded no detection, but sets an upper limit on the neutron star ellipticity of $2.45\times10^{-8}$. We also detected X-ray emission from the redback PSR J1803$-$6707 in data from the first eROSITA all-sky survey, likely due to emission from an intra-binary shock., Comment: 17 pages, 8 figures, accepted for publication in MNRAS

We present the largest single survey to date of average profiles of radio pulsars, observed and processed using the same telescope and data reduction software. Specifically, we present measurements for 1170 pulsars, observed by the Thousand Pulsar Array (TPA) programme at the 64-dish SARAO MeerKAT radio telescope, in a frequency band from 856 to 1712 MHz. We provide rotation measures (RM), dispersion measures, flux densities and polarization properties. The catalogue includes 254 new RMs that substantially increase the total number of known pulsar RMs. Our integration times typically span over 1000 individual rotations per source. We show that the radio (pseudo)luminosity has a strong, shallow dependence on the spin-down energy, proportional to $\dot{E}^{0.15\pm0.04}$, that contradicts some previous proposals of population synthesis studies. In addition, we find a significant correlation between the steepness of the observed flux density spectra and $\dot{E}$, and correlations of the fractional linear polarization with $\dot{E}$, the spectral index, and the pulse width, which we discuss in the context of what is known about pulsar radio emission and how pulsars evolve with time. On the whole, we do not see significant correlations with the estimated surface magnetic field strength, and the correlations with $\dot{E}$ are much stronger than those with the characteristic age. This finding lends support to the suggestion that magnetic dipole braking may not be the dominant factor for the evolution of pulsar rotation over the lifetimes of pulsars. A public data release of the high-fidelity time-averaged pulse profiles in full polarization accompanies our catalogue., Comment: Accepted for publication in MNRAS. 21 pages, 25 figures, 6 Tables

We present the polarization profiles of 22 pulsars in the globular cluster 47 Tucanae using observations from the MeerKAT radio telescope at UHF-band (544-1088 MHz) and report precise values of dispersion measure (DM) and rotation measure (RM). We use these measurements to investigate the presence of turbulence in electron density and magnetic fields. The structure function of DM shows a break at $\sim 30$ arcsec ($\sim 0.6$ pc at the distance of 47 Tucanae) that suggests the presence of turbulence in the gas in the cluster driven by the motion of wind-shedding stars. On the other hand, the structure function of RM does not show evidence of a break. This non-detection could be explained either by the limited number of pulsars or by the effects of the intervening gas in the Galaxy along the line of sight. Future pulsar discoveries in the cluster could help confirm the presence and localise the turbulence., Comment: Accepted for publication in MNRAS, 14 pages, 6 figures

The TRAPUM collaboration has used the MeerKAT telescope to conduct a search for pulsed radio emission from the young Small Magellanic Cloud pulsar J0058-7218 located in the supernova remnant IKT 16, following its discovery in X-rays with XMM-Newton. We report no significant detection of dispersed, pulsed radio emission from this source in three 2-hour L-band observations using the core dishes of MeerKAT, setting an upper limit of 7.0 {\mu}Jy on its mean flux density at 1284 MHz. This is nearly 7 times deeper than previous radio searches for this pulsar in Parkes L-band observations. This suggests that the radio emission of PSR J0058-7218 is not beamed towards Earth or that PSR J0058-7218 is similar to a handful of Pulsar Wind Nebulae systems that have a very low radio efficiency, such as PSR B0540-6919, the Large Magellanic Cloud Crab pulsar analogue. We have also searched for bright, dispersed, single radio pulses and found no candidates above a fluence of 93 mJy ms at 1284 MHz., Comment: 7 pages, 2 figures, 1 table. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The Double Pulsar, PSR J0737-3039A/B, has offered a wealth of gravitational experiments in the strong-field regime, all of which GR has passed with flying colours. In particular, among current gravity experiments that test photon propagation, the Double Pulsar probes the strongest spacetime curvature. Observations with MeerKAT and, in future, the SKA can greatly improve the accuracy of current tests and facilitate tests of NLO contributions in both orbital motion and signal propagation. We present our timing analysis of new observations of PSR J0737-3039A, made using the MeerKAT telescope over the last 3 years. The increased timing precision offered by MeerKAT yields a 2 times better measurement of Shapiro delay parameter s and improved mass measurements compared to previous studies. In addition, our results provide an independent confirmation of the NLO signal propagation effects and already surpass the previous measurement from 16-yr data by a factor of 1.65. These effects include the retardation effect due to the movement of B and the deflection of the signal by the gravitational field of B. We also investigate novel effects which are expected. For instance, we search for potential profile variations near superior conjunctions caused by shifts of the line-of-sight due to latitudinal signal deflection and find insignificant evidence with our current data. With simulations, we find that the latitudinal deflection delay is unlikely to be measured with timing because of its correlation with Shapiro delay. Furthermore, although it is currently not possible to detect the expected lensing correction to the Shapiro delay, our simulations suggest that this effect may be measured with the full SKA. Finally, we provide an improved analytical description for the signal propagation in the Double Pulsar system that meets the timing precision expected from future instruments such as the full SKA., Comment: 13 pages, 13 figures, accepted for publication in A&A

We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1'$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy's measured properties are consistent with other FRB hosts with secure associations., Comment: 14 pages, 9 figures, 3 tables, accepted for publication in MNRAS

We report 4 new pulsars discovered in the core-collapsed globular cluster (GC) NGC 6624 by the TRAPUM Large Survey Project with the MeerKAT telescope. All of the new pulsars found are isolated. PSR J1823$-$3021I and PSR J1823$-$3021K are millisecond pulsars with period of respectively 4.319 ms and 2.768 ms. PSR J1823$-$3021J is mildly recycled with a period of 20.899 ms, and PSR J1823$-$3022 is a long period pulsar with a period of 2.497 s. The pulsars J1823$-$3021I, J1823$-$3021J, and J1823$-$3021K have position and dispersion measure (DM) compatible with being members of the GC and are therefore associated with NGC 6624. Pulsar J1823$-$3022 is the only pulsar bright enough to be re-detected in archival observations of the cluster. This allowed the determination of a timing solution that spans over two decades. It is not possible at the moment to claim the association of pulsar J1823$-$3022 with the GC given the long period and large offset in position ($\sim 3$ arcminutes) and DM (with a fractional difference of 11 percent compared the average of the pulsars in NGC 6624). The discoveries made use of the beamforming capability of the TRAPUM backend to generate multiple beams in the same field of view which allows sensitive searches to be performed over a few half-light radii from the cluster center and can simultaneously localise the discoveries. The discoveries reflect the properties expected for pulsars in core-collapsed GCs., Comment: Accepted for publication on Monthly Notices of the Royal Astronomical Society. 11 pages, 6 figures

MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsars (MSPs) to high precision (< 1 $\mu$s) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing Array (IPTA). In order to plan for the remainder of the programme and to use the allocated time most efficiently, we have conducted an initial census with the MeerKAT "L-band" receiver of 189 MSPs visible to MeerKAT and here present their dispersion measures, polarization profiles, polarization fractions, rotation measures, flux density measurements, spectral indices, and timing potential. As all of these observations are taken with the same instrument (which uses coherent dedispersion, interferometric polarization calibration techniques, and a uniform flux scale), they present an excellent resource for population studies. We used wideband pulse portraits as timing standards for each MSP and demonstrated that the MeerTime Pulsar Timing Array (MPTA) can already contribute significantly to the IPTA as it currently achieves better than 1 $\mu$s timing accuracy on 89 MSPs (observed with fortnightly cadence). By the conclusion of the initial five-year MeerTime programme in July 2024, the MPTA will be extremely significant in global efforts to detect the gravitational wave background with a contribution to the detection statistic comparable to other long-standing timing programmes., Comment: Accepted to PASA. 27 figures. Data to be made available under the DOI 10.5281/zenodo.5347875 at the time of publication

Using the MeerKAT radio telescope, a series of observations have been conducted to time the known pulsars and search for new pulsars in the globular cluster NGC 6440. As a result, two pulsars have been discovered, NGC 6440G and NGC 6440H, one of which is isolated and the other a non-eclipsing (at frequencies above 962 MHz) "Black Widow", with a very low mass companion (M$_{\rm c}$ > 0.006 M$_{\odot}$). It joins the other binary pulsars discovered so far in this cluster which all have low companion masses (M$_{\rm c}$ < 0.30 M$_{\odot}$). We present the results of long-term timing solutions obtained using data from both Green Bank and MeerKAT telescopes for these two new pulsars and an analysis of the pulsars NGC 6440C and NGC 6440D. For the isolated pulsar NGC 6440C, we searched for planets using a Markov Chain Monte Carlo technique. We find evidence for significant unmodelled variations but they cannot be well modelled as planets nor as part of a power-law red-noise process. Studies of the eclipses of the "Redback" pulsar NGC 6440D at two different frequency bands reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies (962-1283 MHz)., Comment: 14 pages, 10 figures, 4 tables. Accepted for publication in MNRAS

We report the discovery of 13 new pulsars in the globular cluster NGC 1851 by the TRAPUM Large Survey Project using the MeerKAT radio telescope. The discoveries consist of six isolated millisecond pulsars (MSPs) and seven binary pulsars, of which six are MSPs and one is mildly recycled. For all the pulsars, we present the basic kinematic, astrometric, and orbital parameters, where applicable, as well as their polarimetric properties, when these are measurable. Two of the binary MSPs (PSR J0514-4002D and PSR J0514-4002E) are in wide and extremely eccentric (e > 0.7) orbits with a heavy white dwarf and a neutron star as their companion, respectively. With these discoveries, NGC 1851 is now tied with M28 as the cluster with the third largest number of known pulsars (14). Its pulsar population shows remarkable similarities with that of M28, Terzan 5 and other clusters with comparable structural parameters. The newly-found pulsars are all located in the innermost regions of NGC 1851 and will likely enable, among other things, detailed studies of the cluster structure and dynamics., Comment: 12 pages, 6 figures, 3 tables. Accepted for publication in Astronomy & Astrophysics

PSR J0955$-$6150 is a member of a class of eccentric MSP+He WD systems (eMSPs), whose binary evolution is poorly understood and believed to be different to that of traditional MSP+He WD systems. Measuring the masses of the stars in this system is important for testing hypotheses for the formation of eMSPs. We have carried out observations of this pulsar with the Parkes and MeerKAT radio telescopes. Our observations reveal a strong frequency evolution of this pulsar's intensity, with a spectral index ($\alpha$) of $-3.13(2)$. The sensitivity of MeerKAT has resulted in a $>10$-fold improvement in the timing precision compared to older Parkes observations. Combined with the 8-year timing baseline, it has allowed precise measurements of a proper motion and three orbital "post-Keplerian" parameters: the rate of advance of periastron, $\dot{\omega} = 0.00152(1) \, {\rm deg} \, yr^{-1}$ and the orthometric Shapiro delay parameters, $h_3 = 0.89(7) \, \mu$s and $\varsigma = 0.88(2)$. Assuming general relativity, we obtain $M_{p} = 1.71(2) \, M_{\odot}$ for the mass of the pulsar and $M_{c} = 0.254(2) \, M_{\odot}$ for the mass of the companion; the orbital inclination is 83.2(4) degrees. We find that the spin axis has a misalignment relative to the orbital angular momentum of $> 4.8$ degrees at 99% CI. While the value of $M_{\rm p}$ falls within the wide range observed in eMSPs, $M_{\rm c}$ is significantly smaller than expected, allowing several formation hypotheses being ruled out. $M_{\rm c}$ is also significantly different from the expected value for an ideal low mass X-ray binary evolution scenario. The putative misalignment between the spin axis of the pulsar and the orbital angular momentum suggests that the unknown process that created the orbital eccentricity of the binary was also capable of changing its orbital orientation, an important evidence for understanding the origin of eMSPs., Comment: 20 pages, 9 figures, 5 tables, accepted for publication by the A&A

The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total intensity and spectral index mosaics of the region from 20 pointings (144 hours on-target in total), covering 6.5 square degrees with an angular resolution of 4$"$,at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT's high sensitivity, exceptional $u,v$-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new non-thermal filament complexes, and enhanced views of the Radio Arc Bubble, Sgr A and Sgr B regions. This project is a SARAO public legacy survey, and the image products are made available with this article., Comment: 29 pages, 15 figures, 3 tables, accepted for publication in the The Astrophysical Journal. Replacement arXiv version with higher quality figures

The Magellanic Clouds are the only external galaxies known to host radio pulsars. The dispersion and rotation measures of pulsars in the Clouds can aid in understanding their structure, and studies of the pulsars themselves can point to potential differences between them and their Galactic counterparts. We use the high sensitivity of the MeerKAT telescope to observe 17 pulsars in the Small and Large Magellanic Clouds in addition to five foreground (Galactic) pulsars. We provide polarization profiles for 18 of these pulsars, improved measurements of their dispersion and rotation measures, and derive the mean parallel magnetic field along the lines of sight. The results are broadly in agreement with expectations for the structure and strength of the magnetic field in the Large and Small Magellanic Clouds. The Magellanic Cloud pulsars have profiles which are narrower than expected from the period-width relationship and we show this is due to selection effects in pulsar surveys rather than any intrinsic difference between the population of Galactic and Magellanic objects., Comment: Accepted for publication in MNRAS

MeerKAT's large number of antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L-band (900-1670 MHz) observations of 115 galaxy clusters, observed for $\sim$6-10 hours each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at $\sim$8" resolution, and enhanced spectral and polarisation image cubes at $\sim$8" and 15" resolutions. Typical sensitivities for the full-resolution MGCLS image products are $\sim$3-5 {\mu}Jy/beam. The basic cubes are full-field and span 4 deg^2. The enhanced products consist of the inner 1.44 deg^2 field of view, corrected for the primary beam. The survey is fully sensitive to structures up to $\sim$10' scales and the wide bandwidth allows spectral and Faraday rotation mapping. HI mapping at 209 kHz resolution can be done at $0

We present pulse width measurements for a sample of radio pulsars observed with the MeerKAT telescope as part of the Thousand-Pulsar-Array (TPA) programme in the MeerTime project. For a centre frequency of 1284 MHz, we obtain 762 $W_{10}$ measurements across the total bandwidth of 775 MHz, where $W_{10}$ is the width at the 10% level of the pulse peak. We also measure about 400 $W_{10}$ values in each of the four or eight frequency sub-bands. Assuming, the width is a function of the rotation period P, this relationship can be described with a power law with power law index $\mu=-0.29\pm 0.03$. However, using orthogonal distance regression, we determine a steeper power law with $\mu=-0.63\pm 0.06$. A density plot of the period-width data reveals such a fit to align well with the contours of highest density. Building on a previous population synthesis model, we obtain population-based estimates of the obliquity of the magnetic axis with respect to the rotation axis for our pulsars. Investigating the width changes over frequency, we unambiguously identify a group of pulsars that have width broadening at higher frequencies. The measured width changes show a monotonic behaviour with frequency for the whole TPA pulsar population, whether the pulses are becoming narrower or broader with increasing frequency. We exclude a sensitivity bias, scattering and noticeable differences in the pulse component numbers as explanations for these width changes, and attempt an explanation using a qualitative model of five contributing Gaussian pulse components with flux density spectra that depend on their rotational phase., Comment: Accepted for publication in MNRAS. 20 pages, 19 figures, 7 tables

We have measured the scattering timescale, $\tau$, and the scattering spectral index, $\alpha$, for 84 single-component pulsars. Observations were carried out with the MeerKAT telescope as part of the Thousand-Pulsar-Array programme in the MeerTime project at frequencies between 0.895 and 1.670 GHz. Our results give a distribution of values for $\alpha$ (defined in terms of $\tau$ and frequency $\nu$ as $\tau\propto\nu^{-\alpha}$) for which, upon fitting a Gaussian, we obtain a mean and standard deviation of $\langle\alpha\rangle = 4.0 \pm 0.6$. This is due to our identification of possible causes of inaccurate measurement of $\tau$, which, if not filtered out of modelling results, tend to lead to underestimation of $\alpha$. The pulsars in our sample have large dispersion measures and are therefore likely to be distant. We find that a model using an isotropic scatter broadening function is consistent with the data, likely due to the averaging effect of multiple scattering screens along the line of sight. Our sample of scattering parameters provides a strong data set upon which we can build to test more complex and time-dependent scattering phenomena, such as extreme scattering events., Comment: Accepted for publication in MNRAS. 14 pages, 14 figures, 2 tables

Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT optical telescope in the days and months after the April 28 event. We did not detect any significant single radio pulses down to fluence limits between 25 mJy ms and 18 Jy ms. Some observing epochs overlapped with times when X-ray bursts were detected. Radio images made on four days using the MeerKAT telescope revealed no point-like persistent or transient emission at the location of the magnetar. No transient or persistent optical emission was detected over seven days. Using the multi-colour MeerLICHT images combined with relations between DM, NH and reddening we constrain the distance to SGR J1935+2154, to be between 1.5 and 6.5 kpc. The upper limit is consistent with some other distance indicators and suggests that the April 28 burst is closer to two orders of magnitude less energetic than the least energetic FRBs. The lack of single-pulse radio detections shows that the single pulses detected over a range of fluences are either rare, or highly clustered, or both. It may also indicate that the magnetar lies somewhere between being radio-quiet and radio-loud in terms of its ability to produce radio emission efficiently., Comment: Accepted for publication in MNRAS. Corresponding author B. W. Stappers

We have used the central 44 antennas of the new 64-dish MeerKAT radio telescope array to conduct a deep search for new pulsars in the core of nine globular clusters. This has led to the discovery of eight new millisecond pulsars in six different clusters. Two new binaries, 47 Tuc ac and 47 Tuc ad, are eclipsing "spiders", featuring compact orbits ($\lesssim 0.32$ days), very low-mass companions and regular occultations of their pulsed emission. The other three new binary pulsars (NGC 6624G, M62G, and Ter 5 an) are in wider ($> 0.7$ days) orbits, with companions that are likely to be white dwarfs or neutron stars. NGC 6624G has a large eccentricity of $e\simeq 0.38$, which enabled us to detect the rate of advance of periastron. This suggests that the system is massive, with a total mass of $M{\rm tot} = 2.65 \pm 0.07$ M$_{\odot}$. Likewise, for Ter 5 an, with $e \simeq 0.0066$, we obtain $M{\rm tot}= 2.97 \pm 0.52$ M$_{\odot}$. The other three new discoveries (NGC 6522D, NGC 6624H and NGC 6752F) are faint isolated pulsars. Finally, we have used the whole MeerKAT array and synthesized 288 beams, covering an area of $\sim2$ arcmin in radius around the center of NGC 6624. This has allowed us to localize many of the pulsars in the cluster, demonstrating the beamforming capabilities of the TRAPUM software backend and paving the way for the upcoming MeerKAT globular cluster pulsar survey., Comment: 20 pages, 10 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We describe the ongoing Relativistic Binary programme (RelBin), a part of the MeerTime large survey project with the MeerKAT radio telescope. RelBin is primarily focused on observations of relativistic effects in binary pulsars to enable measurements of neutron star masses and tests of theories of gravity. We selected 25 pulsars as an initial high priority list of targets based on their characteristics and observational history with other telescopes. In this paper, we provide an outline of the programme, present polarisation calibrated pulse profiles for all selected pulsars as a reference catalogue along with updated dispersion measures. We report Faraday rotation measures for 24 pulsars, twelve of which have been measured for the first time. More than a third of our selected pulsars show a flat position angle swing confirming earlier observations. We demonstrate the ability of the Rotating Vector Model (RVM), fitted here to seven binary pulsars, including the Double Pulsar (PSR J0737$-$3039A), to obtain information about the orbital inclination angle. We present a high time resolution light curve of the eclipse of PSR J0737$-$3039A by the companion's magnetosphere, a high-phase resolution position angle swing for PSR J1141$-$6545, an improved detection of the Shapiro delay of PSR J1811$-$2405, and pulse scattering measurements for PSRs J1227$-$6208, J1757$-$1854, and J1811$-$1736. Finally, we demonstrate that timing observations with MeerKAT improve on existing data sets by a factor of, typically, 2-3, sometimes by an order of magnitude., Comment: 21 pages 16 figures, published in MNRAS (replaced earlier submission after small changes added in proofs)

Using the state-of-the-art SKA precursor, the MeerKAT radio telescope, we explore the limits to precision pulsar timing of millisecond pulsars achievable due to pulse stochasticity (jitter). We report new jitter measurements in 15 of the 29 pulsars in our sample and find that the levels of jitter can vary dramatically between them. For some, like the 2.2~ms pulsar PSR J2241--5236, we measure an implied jitter of just $\sim$ 4~ns/hr, while others like the 3.9~ms PSR J0636--3044 are limited to $\sim$ 100 ns/hr. While it is well known that jitter plays a central role to limiting the precision measurements of arrival times for high signal-to-noise ratio observations, its role in the measurement of dispersion measure (DM) has not been reported, particularly in broad-band observations. Using the exceptional sensitivity of MeerKAT, we explored this on the bright millisecond pulsar PSR J0437--4715 by exploring the DM of literally every pulse. We found that the derived single pulse DMs vary by typically 0.0085 cm$^{-3}$ pc from the mean, and that the best DM estimate is limited by the differential pulse jitter across the band. We postulate that all millisecond pulsars will have their own limit on DM precision which can only be overcome with longer integrations. Using high-time resolution filterbank data of 9 $\mu$s, we also present a statistical analysis of single pulse phenomenology. Finally, we discuss optimization strategies for the MeerKAT pulsar timing program and its role in the context of the International Pulsar Timing Array (IPTA)., Comment: 16 pages, 9 figures

The Thousand Pulsar Array (TPA) project currently monitors about 500 pulsars with the sensitive MeerKAT radio telescope by using subarrays to observe multiple sources simultaneously. Here we define the adopted observing strategy, which guarantees that each target is observed long enough to obtain a high fidelity pulse profile, thereby reaching a sufficient precision of a simple pulse shape parameter. This precision is estimated from the contribution of the system noise of the telescope, and the pulse-to-pulse variability of each pulsar, which we quantify under some simplifying assumptions. We test the assumptions and choice of model parameters using data from the MeerKAT 64-dish array, Lovell and Parkes telescopes. We demonstrate that the observing times derived from our method produce high fidelity pulse profiles that meet the needs of the TPA in studying pulse shape variability and pulsar timing. Our method can also be used to compare strategies for observing large numbers of pulsars with telescopes capable of forming multiple subarray configurations. We find that using two 32-dish MeerKAT subarrays is the most efficient strategy for the TPA project. We also find that the ability to observe in different array configurations will become increasingly important for large observing programmes using the Square Kilometre Array telescope., Comment: 13 pages, 8 figures. Accepted for publication in MNRAS

We present observations of 35 high spin-down energy radio pulsars using the MeerKAT telescope. Polarisation profiles and associated parameters are also presented. We derive the geometry for a selection of pulsars which show interpulse emission. We point out that, in several cases, these radio pulsars should also be seen in $\gamma$-rays but that improved radio timing is required to aid the high-energy detection. We discuss the relationship between the width of the radio profile and its high-energy detectability. Finally, we reflect on the correlation between the spin-down energy and the radio polarisation fraction and the implications this may have for $\gamma$-ray emission., Comment: 11 pages, 4 figures, accepted for publication by the MNRAS

The millisecond pulsar J1823-3021A is a very active giant pulse emitter in the globular cluster NGC 6624. New observations with the MeerKAT radio telescope have revealed 14350 giant pulses over 5 hours of integration time, with an average wait time of about 1 second between giant pulses. The giant pulses occur in phases compatible with the ordinary radio emission, follow a power-law distribution with an index of -2.63 $\pm$ 0.02 and contribute 4 percent of the total integrated flux. The spectral index of the giant pulses follows a Gaussian distribution centered around -1.9 with a standard deviation of 0.6 and is on average flatter than the integrated emission, which has a spectral index of -2.81 $\pm$ 0.02. The waiting times between the GPs are accurately described by a Poissonian distribution, suggesting that the time of occurrence of a GP is independent from the times of occurrence of other GPs. 76 GPs show multiple peaks within the same rotation, a rate that is also compatible with the mutual independence of the GP times of occurrence. We studied the polarization properties of the giant pulses finding, on average, linear polarization only at the 1 percent level and circular polarization at the 3 percent level, similar to the polarization percentages of the total integrated emission. In 4 cases it was possible to measure the RM of the GPs which are highly variable and, in two cases, is inconsistent with the mean RM of the total integrated pulsar signal., Comment: Accepted for publication in MNRAS. 8 Pages, 11 Figures

We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time., Comment: Accepted for publication in MNRAS