Search

We use data from the MeerTime project on the MeerKAT telescope to ask whether the radio emission properties of millisecond pulsars (MSPs) and slowly rotating, younger pulsars (SPs) are similar or different. We show that the flux density spectra of both populations are similarly steep, and the widths of MSP pulsar profiles obey the same dependence on the rotational period as slow pulsars. We also show that the polarization of MSPs has similar properties to slow pulsars. The commonly used pseudo-luminosity of pulsars, defined as the product of the flux density and the distance squared, is not appropriate for drawing conclusions about the relative intrinsic radio luminosity of SPs and MSPs. We show that it is possible to scale the pseudo-luminosity to account for the pulse duty cycle and the solid angle of the radio beam, in such a way that MSPs and SPs do not show clear differences in intrinsic luminosity. The data, therefore, support common emission physics between the two populations in spite of orders of magnitude difference in their period derivatives and inferred, surface, dipole magnetic field strengths., Comment: 9 pages, accepted for publication in MNRAS

The S-shaped swing of the linear polarization position angle (PPA) observed in many pulsars can be interpreted by the rotating vector model (RVM). However, efforts to fit the RVM for a large sample of pulsars observed with the MeerKAT telescope as a part of the Thousand-Pulsar-Array (TPA) programme, only succeeded for about half the cases. High time-resolution studies suggest that the failed cases arise due to the presence of orthogonal polarization modes, or highly disordered distribution of PPA points. One such example is PSR~J1645-0317. Recently it has been shown that the RVM can be recovered in this pulsar by using only time samples which are greater than 80% linearly polarized. In this work we test this novel approach on the brightest 249 pulsars from the TPA sample, of which 177 yield sufficient highly polarized samples to be amenable to our method. Remarkably, only 9 of these pulsars (5%) now fail to fit the RVM as opposed to 59% from the original analysis. This result favours the paradigm that the underlying mechanism is coherent curvature radiation., Comment: Delayed interminably by MNRAS

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

A number of pulsars are known to have profile evolution on timescales of months, often correlated with spin-down rate changes. Here, we present the first result from 3 years of monitoring observations from MeerKAT as part of the Thousand Pulsar Array programme. This programme obtains high-fidelity pulse profiles for $\sim$ 500 pulsars, which enabled the detection of subtle changes in seven sources not previously known to exhibit long-term profile evolution. A 2D Gaussian convolution is used to highlight correlated emission variability in both the pulse phase and observing epoch direction. Simulations show that for one additional source the observed profile variability is likely to originate from stochastic single-pulse shape variability (jitter). We find that it is common for long-term profile variability to be associated with changes in polarization fractions, but not with polarisation position angle (PA) changes. PA changes are expected if emission height changes or precession is responsible for the profile variability. PSR J1741$-$3927 is the only pulsar in our sample that shows correlated PA variability, and this is associated with orthogonal polarization mode activity. For the six %the rest, without correlated PA variability, other pulsars limits on possible emission height changes and impact angle changes are derived. These limits are consistent with the small changes in the total intensity profile shape. None of the sources show detectable spin-down variability correlated with the emission changes, which are thought to be driven by magnetospheric current fluctuations. Therefore the absence of correlated spin-down rate variability allows upper limits to be placed on changes in the magnetospheric charge density., Comment: 19 pages, 12 figures, 3 Tables. Accepted for publication in MNRAS

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

The population of radio pulsars is observed to demonstrate certain polarization properties not explained by the conventional picture of pulsar polarization, namely frequency evolution of polarization, deviations of the linear polarization angle from a curve of geometric origins and the presence of features in the circular polarization. We present the partial-coherence model as a way to explain the co-occurrence of these features and to provide an origin for circular polarization in radio pulsar profiles. We describe the mathematics of the model and demonstrate how it can explain these observed features, both on a population level and for the idiosyncrasies of individual pulsars. The partial coherence model can account for complex polarization behaviour, enabling improved access to information about pulsar geometries. We discuss the scientific implications of this for our understanding of pulsar radio emission and propagation., Comment: 15 pages, 10 figures

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

The radio polarization properties of the pulsar population are only superficially captured by the conventional picture of pulsar radio emission. We study the broadband polarization of 271 young radio pulsars, focusing particularly on circular polarization, using high quality observations made with the Ultra-Wideband Low receiver on Murriyang, the Parkes radio telescope. We seek to encapsulate polarization behaviour on a population scale by defining broad categories for frequency- and phase-dependent polarization evolution, studying the co-occurrences of these categorizations and comparing them with average polarization measurements and spin-down energy ($\dot{E}$). This work shows that deviations of the linear polarization position angle (PA) from the rotating vector model (RVM) are linked to the presence of circular polarization features and to frequency evolution of the polarization. Polarization fraction, circular polarization contribution and profile complexity all evolve with $\dot{E}$ across the population, with the profiles of high-$\dot{E}$ pulsars being simple and highly linearly polarized. The relationship between polarization fraction and circular contribution is also seen to evolve such that highly polarized profiles show less variation in circular contribution with frequency than less strongly polarized profiles. This evolution is seen both across the population and across frequency for individual sources. Understanding pulsar radio polarization requires detailed study of individual sources and collective understanding of population-level trends. For the former, we provide visualizations of their phase- and frequency-resolved polarization parameters. For the latter, we have highlighted the importance of including the impact of circular polarization and of $\dot{E}$., Comment: 20 pages, 7 figures, supplementary material available at the MNRAS website or by contacting the lead author

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

We present the first 2.5 years of data from the MeerKAT Pulsar Timing Array (MPTA), part of MeerTime, a MeerKAT Large Survey Project. The MPTA aims to precisely measure pulse arrival times from an ensemble of 88 pulsars visible from the Southern Hemisphere, with the goal of contributing to the search, detection and study of nanohertz-frequency gravitational waves as part of the International Pulsar Timing Array. This project makes use of the MeerKAT telescope, and operates with a typical observing cadence of two weeks using the L-band receiver that records data from 856-1712 MHz. We provide a comprehensive description of the observing system, software, and pipelines used and developed for the MeerTime project. The data products made available as part of this data release are from the 78 pulsars that had at least $30$ observations between the start of the MeerTime programme in February 2019 and October 2021. These include both sub-banded and band-averaged arrival times, as well as the initial timing ephemerides, noise models, and the frequency-dependent standard templates (portraits) used to derive pulse arrival times. After accounting for detected noise processes in the data, the frequency-averaged residuals of $67$ of the pulsars achieved a root-mean-square residual precision of $< 1 \mu \rm{s}$. We also present a novel recovery of the clock correction waveform solely from pulsar timing residuals, and an exploration into preliminary findings of interest to the international pulsar timing community. The arrival times, standards and full Stokes parameter calibrated pulsar timing archives are publicly available., Comment: 18 pages, 6 figures

In spite of the rich phenomenology of the polarization properties of radio pulsars, the rotating vector model (RVM) created 50 years ago remains the best method to determine the beam geometry of a pulsar. We apply the RVM to a sample of 854 radio pulsars observed with the MeerKAT telescope in order to draw conclusions about the population of pulsars as a whole. The main results are that (i) the geometrical interpretation of the position angle traverse is valid in the majority of the population, (ii) the pulsars for which the RVM fails tend to have a high fraction of circular polarization compared to linear polarization, (iii) emission heights obtained through both geometrical and relativistic methods show that the majority of pulsars must have emission heights less than 1000~km independent of spin period, (iv) orthogonal mode jumps are seen in the position angle traverse in about one third of the population. All these results are weakly dependent on the pulsar spin-down energy., Comment: accepted by 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

The LOFAR Tied-Array All-Sky Survey (LOTAAS) is the most sensitive untargeted radio pulsar survey performed at low radio frequencies (119--151\,MHz) to date and has discovered 76 new radio pulsars, among which the 23.5-s pulsar J0250+5854, up until recently the slowest-spinning radio pulsar known. Here, we report on the timing solutions of 35 pulsars discovered by LOTAAS, which include a nulling pulsar and a mildly recycled pulsar, and thereby complete the full timing analysis of the LOTAAS pulsar discoveries. We give an overview of the findings from the full LOTAAS sample of 76 pulsars, discussing their pulse profiles, radio spectra and timing parameters. We found that the pulse profiles of some of the pulsars show profile variations in time or frequency and while some pulsars show signs of scattering, a large majority display no pulse broadening. The LOTAAS discoveries have on average steeper radio spectra and have longer spin periods ($1.4\times$) as well as lower spin-down rates ($3.1\times$) compared to the known pulsar population. We discuss the cause of these differences, and attribute them to a combination of selection effects of the LOTAAS survey as well as previous pulsar surveys, though can not rule out that older pulsars tend to have steeper radio spectra., Comment: Accepted to Astronomy & Astrophysics

We present the detection of 107 pulsars with interstellar scintillation arcs at 856--1712\,MHz, observed with the MeerKAT Thousand Pulsar Array Programme. Scintillation arcs appear to be ubiquitous in clean, high S/N observations, their detection mainly limited by short observing durations and coarse frequency channel resolution. This led the survey to be sensitive to nearby, lightly scattered pulsars with high effective velocity -- from a large proper motion, a screen nearby the pulsar, or a screen near the Earth. We measure the arc curvatures in all of our sources, which can be used to give an estimate of screen distances in pulsars with known proper motion, or an estimate of the proper motion. The short scintillation timescale in J1731$-$4744 implies a scattering screen within 12\,pc of the source, strongly suggesting the association between this pulsar and the supernova remnant RCW 114. We measure multiple parabolic arcs of 5 pulsars, all of which are weakly scintillating with high proper motion. Additionally, several sources show hints of inverted arclets suggesting scattering from anisotropic screens. Building on this work, further targeted MeerKAT observations of many of these pulsars will improve understanding of our local scattering environment and the origins of scintillation; annual scintillation curves would lead to robust screen distance measurements, and the evolution of arclets in time and frequency can constrain models of scintillation., Comment: 12 pages, 5 figures, 2 tables. Accepted in MNRAS

This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwidths, forming several 100s of beams, and processing the data in the SKA real-time pulsar search pipelines. A brief account of the different architectures of the accelerators, primarily the new generation Field Programmable Gate Array-based accelerators, showing their critical roles to achieve high-performance computing and in handling the enormous data volume problems of the SKA is presented here. It also presents the power-performance efficiency of this emerging technology and presents potential future scenarios., Comment: Accepted for publication in the JoAA special issue on SKA (2023)

With the ever-increasing sensitivity and timing baselines of modern radio telescopes, a growing number of pulsars are being shown to exhibit transitions in their rotational and radio emission properties. In many of these cases, the two are correlated with pulsars assuming a unique spin-down rate ($\dot{\nu}$) for each of their specific emission states. In this work we revisit 17 radio pulsars previously shown to exhibit spin-down rate variations. Using a Gaussian process regression (GPR) method to model the timing residuals and the evolution of the profile shape, we confirm the transitions already observed and reveal new transitions in 8 years of extended monitoring with greater time resolution and enhanced observing bandwidth. We confirm that 7 of these sources show emission-correlated $\dot{\nu}$ transitions ($\Delta \dot{\nu}$) and we characterise this correlation for one additional pulsar, PSR B1642$-$03. We demonstrate that GPR is able to reveal extremely subtle profile variations given sufficient data quality. We also corroborate the dependence of $\Delta \dot{\nu}$ amplitude on $\dot{\nu}$ and pulsar characteristic age. Linking $\Delta \dot{\nu}$ to changes in the global magnetospheric charge density $\Delta \rho$, we speculate that $\dot{\nu}$ transitions associated with large $\Delta \rho$ values may be exhibiting detectable profile changes with improved data quality, in cases where they have not previously been observed., Comment: 20 pages, 16 figures, 3 tables

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

Discovering radio pulsars, particularly millisecond pulsars (MSPs), is important for a range of astrophysical applications, such as testing theories of gravity or probing the magneto-ionic interstellar medium. We aim to discover pulsars that may have been missed in previous pulsar searches by leveraging known pulsar observables (primarily polarisation) in the sensitive, low-frequency radio images from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS), and have commenced the Targeted search, using LoTSS images, for polarised pulsars (TULIPP) survey. For this survey, we identified linearly and circularly polarised point sources with flux densities brighter than 2 mJy in LoTSS images at a centre frequency of 144 MHz with a 48 MHz bandwidth. Over 40 known pulsars, half of which are MSPs, were detected as polarised sources in the LoTSS images and excluded from the survey. We have obtained beam-formed LOFAR observations of 30 candidates, which were searched for pulsations using coherent de-dispersion. Here, we present the results of the first year of the TULIPP survey. We discovered two pulsars, PSRs J1049+5822 and J1602+3901, with rotational periods of P=0.73 s and 3.7 ms, respectively. We also detected a further five known pulsars (two slowly-rotating pulsars and three MSPs) for which accurate sky positions were not available to allow a unique cross-match with LoTSS sources. This targeted survey presents a relatively efficient method by which pulsars, particularly MSPs, may be discovered using the flexible observing modes of sensitive radio telescopes such as the Square Kilometre Array and its pathfinders/precursors, particularly since wide-area all-sky surveys using coherent de-dispersion are currently computationally infeasible., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in A&A

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 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

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array., Comment: 24 pages, 6 figures. Updated to match published version

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 present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this object, with the detections at 1250 MHz (FAST) and 57 MHz (NenuFAR) being the highest- and lowest-frequency published respectively to date. We measure a flux density of $4\pm2$ $\mu$Jy at 1250 MHz and an exceptionally steep spectral index of $-3.5^{+0.2}_{-1.5}$, with a turnover below $\sim$95 MHz. In conjunction with observations of this pulsar with the GBT and the LOFAR Core, we show that the intrinsic profile width increases drastically towards higher frequencies, contrary to the predictions of conventional radius-to-frequency mapping. We examine polarimetric data from FAST and the LOFAR Core and conclude that its polar cap radio emission is produced at an absolute height of several hundreds of kilometres around 1.5 GHz, similar to other rotation-powered pulsars across the population. Its beam is significantly underfilled at lower frequencies, or it narrows because of the disappearance of conal outriders. Finally, the results for PSR J0250+5854 and other slowly spinning rotation-powered pulsars are contrasted with the radio-detected magnetars. We conclude that magnetars have intrinsically wider radio beams than the slow rotation-powered pulsars, and that consequently the latter's lower beaming fraction is what makes objects such as PSR J0250+5854 so scarce., Comment: 14 pages, 6 figures

Since their discovery more than 50 years ago, broadband radio studies of pulsars have generated a wealth of information about the underlying physics of radio emission. In order to gain some further insights into this elusive emission mechanism, we performed a multi-frequency study of two very well-known pulsars, PSR~B0919+06 and PSR~B1859+07. These pulsars show peculiar radio emission properties whereby the emission shifts to an earlier rotation phase before returning to the nominal emission phase in a few tens of pulsar rotations (also known as `swooshes'). We confirm the previous claim that the emission during the swoosh is not necessarily absent at low frequencies and the single pulses during a swoosh show varied behaviour at 220~MHz. We also confirm that in PSR~B0919+06, the pulses during the swoosh show a chromatic dependence of the maximum offset from the normal emission phase with the offset following a consistent relationship with observing frequency. We also observe that the flux density spectrum of the radio profile during the swoosh is inverted compared to the normal emission. For PSR~B1859+07, we have discovered a new mode of emission in the pulsar that is potentially quasi-periodic with a different periodicity than is seen in its swooshes. We invoke an emission model previously proposed in the literature and show that this simple model can explain the macroscopic observed characteristics in both pulsars. We also argue that pulsars that exhibit similar variability on short timescales may have the same underlying emission mechanism., Comment: 13 pages, 13 figures, 1 table, accepted for publication in MNRAS

PSR J0540$-$6919 is the second-most energetic radio pulsar known and resides in the Large Magellanic Cloud. Like the Crab pulsar it is observed to emit giant radio pulses (GPs). We used the newly-commissioned PTUSE instrument on the MeerKAT radio telescope to search for GPs across three observations. In a total integration time of 5.7 hrs we detected 865 pulses above our 7$\sigma$ threshold. With full polarisation information for a subset of the data, we estimated the Faraday rotation measure, $\rm{RM}=-245.8 \pm 1.0$ rad m$^{-2}$ toward the pulsar. The brightest of these pulses is $\sim$ 60% linearly polarised but the pulse-to-pulse variability in the polarisation fraction is significant. We find that the cumulative GP flux distribution follows a power law distribution with index $-2.75 \pm 0.02$. Although the detected GPs make up only $\sim$ 10% of the mean flux, their average pulse shape is indistinguishable from the integrated pulse profile, and we postulate that there is no underlying emission. The pulses are scattered at L-band frequencies with the brightest pulse exhibiting a scattering time-scale of $\tau = 0.92 \pm 0.02$ ms at 1.2 GHz. We find several of the giants display very narrow-band "flux knots" similar to those seen in many Fast Radio Bursts, which we assert cannot be due to scintillation or plasma lensing. The GP time-of-arrival distribution is found to be Poissonian on all but the shortest time-scales where we find four GPs in six rotations, which if GPs are statistically independent is expected to occur in only 1 of 7000 observations equivalent to our data., Comment: 15 pages, 10 figures

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

The low-frequency radio spectra of the hotspots within powerful radio galaxies can provide valuable information about the physical processes operating at the site of the jet termination. These processes are responsible for the dissipation of jet kinetic energy, particle acceleration, and magnetic-field generation. Here we report new observations of the powerful radio galaxy Cygnus A using the Low Frequency Array (LOFAR) between 109 and 183 MHz, at an angular resolution of ~3.5 arcsec. The radio emission of the lobes is found to have a complex spectral index distribution, with a spectral steepening found towards the centre of the source. For the first time, a turnover in the radio spectrum of the two main hotspots of Cygnus A has been directly observed. By combining our LOFAR imaging with data from the Very Large Array at higher frequencies, we show that the very rapid turnover in the hotspot spectra cannot be explained by a low-energy cut-off in the electron energy distribution, as has been previously suggested. Thermal (free-free) absorption or synchrotron self absorption models are able to describe the low-frequency spectral shape of the hotspots, however, as with previous studies, we find that the implied model parameters are unlikely, and interpreting the spectra of the hotspots remains problematic., Comment: 9 pages, 3 figures, accepted for publication in MNRAS

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

Binary supermassive black hole (SMBH) systems result from galaxy mergers, and will eventually coalesce due to gravitational wave (GW) emission if the binary separation can be reduced to $\lesssim 0.1$ pc by other mechanisms. Here, we explore a gravitational self-lensing binary SMBH model for the sharp (duration $\sim 1$ hr), quasi-regular X-ray flares -- dubbed quasi-periodic eruptions -- recently observed from two low mass active galactic nuclei: GSN 069 and RX J1301.9+2747. In our model, the binary is observed $\sim$edge-on, such that each SMBH gravitationally lenses light from the accretion disc surrounding the other SMBH twice per orbital period. The model can reproduce the flare spacings if the current eccentricity of RX J1301.9+2747 is $\epsilon_0 \gtrsim 0.16$, implying a merger within $\sim 1000$ yrs. However, we cannot reproduce the observed flare profiles with our current calculations. Model flares with the correct amplitude are $\sim 2/5$ the observed duration, and model flares with the correct duration are $\sim 2/5$ the observed amplitude. Our modelling yields three distinct behaviours of self-lensing binary systems that can be searched for in current and future X-ray and optical time-domain surveys: i) periodic lensing flares, ii) partial eclipses (caused by occultation of the background mini-disc by the foreground mini-disc), and iii) partial eclipses with a very sharp in-eclipse lensing flare. Discovery of such features would constitute very strong evidence for the presence of a supermassive binary, and monitoring of the flare spacings will provide a measurement of periastron precession., Comment: Accepted for publication in MNRAS

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

FRB~121102 was the first Fast Radio Burst (FRB) that was shown to repeat. Since its discovery in 2012, more than two hundred bursts have been detected from the source. These bursts exhibit a diverse range of spectral and temporal characteristics and many questions about their origin and form remain unanswered. Here, we present a sample of radio bursts from FRB 121102 detected using the Lovell telescope at Jodrell Bank Observatory. We show four examples of bursts that show peculiar spectro-temporal characteristics and compare them with properties of bursts of FRB~121102 detected at other observatories. We report on a precursor burst that is separated by just 17~ms from the main burst, the shortest reported separation between two individual bursts to date. We also provide access to data for all the detections of FRB~121102 in this campaign., Comment: 3 pages, 1 Figure, to be published in RNAAS

The population of young, non-recycled pulsars with spin down energies Edot >10^35 erg/s is sampled predominantly at gamma-ray and radio wavelengths. A total of 137 such pulsars are known, with partial overlap between the sources detectable in radio and gamma-rays. We use a very small set of assumptions in an attempt to test whether the observed pulsar sample can be explained by a single underlying population of neutron stars. For radio emission we assume a canonical conal beam with a fixed emission height of 300~km across all spin periods and a luminosity law which depends on Edot^{0.25}. For gamma-ray emission we assume the outer-gap model and a luminosity law which depends on Edot^{0.5}. We synthesise a population of fast-spinning pulsars with a birth rate of one per 100 years. We find that this simple model can reproduce most characteristics of the observed population with two caveats. The first is a deficit of gamma-ray pulsars at the highest Edot which we surmise to be an observational selection effect due to the difficulties of finding gamma-ray pulsars in the presence of glitches without prior knowledge from radio frequencies. The second is a deficit of radio pulsars with interpulse emission, which may be related to radio emission physics. We discuss the implications of these findings., Comment: Accepted by MNRAS

Pulsar radio emission and its polarization are observed to evolve with frequency. This frequency dependence is key to the emission mechanism and the structure of the radio beam. With the new Ultra-Wideband receiver (UWL) on the Parkes radio telescope we are able, for the first time, to observe how pulsar profiles evolve over a broad continuous bandwidth of 700-4000 MHz. We describe here a technique for processing broadband polarimetric observations to establish a meaningful alignment and visualize the data across the band. We apply this to observations of PSRs J1056-6258 and J1359-6038, chosen due to previously unresolved questions about the frequency evolution of their emission. Application of our technique reveals that it is possible to align the polarization position angle (PA) across a broad frequency range when constrained to applying only corrections for dispersion and Faraday rotation to do so. However, this does not correspond to aligned intensity profiles for these two sources. We find that it is possible to convert these misalignments into emission height range estimates that are consistent with published and simulated values, suggesting that they can be attributed to relativistic effects in the magnetosphere. We discuss this work in the context of the radio beam structure and prepare the ground for a wider study of pulsar emission using broadband polarimetric data., Comment: Accepted for publication in MNRAS. 13 pages, 10 figures, 2 tables

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., Comment: 20 pages, 16 Figures, 4 Tables, accepted for publication in PASA

We present the data analysis pipeline, commissioning observations and initial results from the GREENBURST fast radio burst (FRB) detection system on the Robert C. Byrd Green Bank Telescope (GBT) previously described by Surnis et al. which uses the 21~cm receiver observing commensally with other projects. The pipeline makes use of a state-of-the-art deep learning classifier to winnow down the very large number of false positive single-pulse candidates that mostly result from radio frequency interference. In our observations totalling 156.5 days so far, we have detected individual pulses from 20 known radio pulsars which provide an excellent verification of the system performance. We also demonstrate, through blind injection analyses, that our pipeline is complete down to a signal-to-noise threshold of 12. Depending on the observing mode, this translates to peak flux sensitivities in the range 0.14--0.89~Jy. Although no FRBs have been detected to date, we have used our results to update the analysis of Lawrence et al. to constrain the FRB all-sky rate to be $1140^{+200}_{-180}$ per day above a peak flux density of 1~Jy. We also constrain the source count index $\alpha=0.83\pm0.06$ which indicates that the source count distribution is substantially flatter than expected from a Euclidean distribution of standard candles (where $\alpha=1.5$). We discuss this result in the context of the FRB redshift and luminosity distributions. Finally, we make predictions for detection rates with GREENBURST, as well as other ongoing and planned FRB experiments., Comment: 9 pages, 7 figures, submitted to MNRAS

The discovery that at least some Fast Radio Bursts (FRBs) repeat has ruled out cataclysmic events as the progenitors of these particular bursts. FRB~121102 is the most well-studied repeating FRB but despite extensive monitoring of the source, no underlying pattern in the repetition has previously been identified. Here, we present the results from a radio monitoring campaign of FRB~121102 using the 76-m Lovell telescope. Using the pulses detected in the Lovell data along with pulses from the literature, we report a detection of periodic behaviour of the source over the span of five years of data. We predict that the source is currently `off' and that it should turn `on' for the approximate MJD range $59002-59089$ (2020-06-02 to 2020-08-28). This result, along with the recent detection of periodicity from another repeating FRB, highlights the need for long-term monitoring of repeating FRBs at a high cadence. Using simulations, we show that one needs at least 100 hours of telescope time to follow-up repeating FRBs at a cadence of 0.5--3 days to detect periodicities in the range of 10--150 days. If the period is real, it shows that repeating FRBs can have a large range in their activity periods that might be difficult to reconcile with neutron star precession models., Comment: 9 pages, 6 Figures, 3 Tables. More detections added to the periodicity search. Updated period and activity period. Complete tables (Table 1 and Table 3) have been submitted as online supplementary materials

Fast Radio Bursts (FRBs) are bright, extragalactic radio pulses whose origins are still unknown. Until recently, most FRBs have been detected at frequencies greater than 1 GHz with a few exceptions at 800 MHz. The recent discoveries of FRBs at 400 MHz from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope has opened up possibilities for new insights about the progenitors while many other low frequency surveys in the past have failed to find any FRBs. Here, we present results from a FRB survey recently conducted at the Jodrell Bank Observatory at 332 MHz with the 76-m Lovell telescope for a total of 58 days. We did not detect any FRBs in the survey and report a 90$\%$ upper limit of 5500 FRBs per day per sky for a Euclidean Universe above a fluence threshold of 46 Jy ms. We discuss the possibility of absorption as the main cause of non-detections in low frequency (< 800 MHz) searches and invoke different absorption models to explain the same. We find that Induced Compton Scattering alone cannot account for absorption of radio emission and that our simulations favour a combination of Induced Compton Scattering and Free-Free Absorption to explain the non-detections. For a free-free absorption scenario, our constraints on the electron density are consistent with those expected in the post-shock region of the ionized ejecta in Super-Luminous SuperNovae (SLSNe)., Comment: 12 pages, 9 Figures, 2 Tables, Second revision submitted to MNRAS

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., Comment: Accepted by MNRAS

The frequency widening of pulsar profiles is commonly attributed to lower frequencies being produced at greater heights above the surface of the pulsar; so-called radius-to-frequency mapping. The observer's view of pulsar emission is a 1D cut through a 3D magnetosphere: we can only see that emission which points along our line of sight. However, by comparing the frequency evolution of many single pulses positioned at different phases, we can build up an understanding of the shape of the active emission region. We use single pulses observed with the Giant Metrewave Radio Telescope to investigate the emission region of PSR J1136+1551 and test radius-to-frequency mapping. Assuming that emission is produced tangential to the magnetic field lines and that each emission frequency corresponds to a single height, we simulate the single pulse profile evolution resulting from the canonical conal beam model and a fan beam model. Comparing the results of these simulations with the observations, we conclude that the emission region of PSR J1136+1551 is better described by the fan beam model. The diversity of profile widening behaviour observed for the single pulses can be explained by orthogonally polarized modes propagating along differing frequency-dependent paths in the magnetosphere., Comment: Accepted for publication in MNRAS. 17 pages, 14 figures, 1 table

We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees), as well as three incoherent beams (covering 67 square degrees) are formed in each survey pointing. For each ofthe 222 beams, total intensity is recorded at 491.52 us time resolution. Each observation integrates for 1 hr and covers 2592 channels from 119 to 151 MHz. This instrumental setup allows LOTAAS to reach a detection threshold of 1 to 5 mJy for periodic emission. Thus far, the LOTAAS survey has resulted in the discovery of 73 radio pulsars. Among these are two mildly recycled binary millisecond pulsars (P = 13 and 33 ms), as well as the slowest-spinning radio pulsar currently known (P = 23.5 s). The survey has thus far detected 311 known pulsars, with spin periods ranging from 4 ms to 5.0 s and dispersion measures from 3.0 to 217 pc/cc. Known pulsars are detected at flux densities consistent with literature values. We find that the LOTAAS pulsar discoveries have, on average, longer spin periods than the known pulsar population. This may reflect different selection biases between LOTAAS and previous surveys, though it is also possible that slower-spinning pulsars preferentially have steeper radio spectra. LOTAAS is the deepest all-sky pulsar survey using a digital aperture array; we discuss some of the lessons learned that can inform the approach for similar surveys using future radio telescopes such as the Square Kilometre Array., Comment: 22 pages, 9 figures, 3 tables. Accepted for publication in A&A

It is now known that the emission from radio pulsars can vary over a wide range of timescales, from fractions of seconds to decades. However, it is not yet known if long- and short-term emission variability are caused by the same physical processes. It has been observed that long-term emission variability is often correlated with rotational changes in the pulsar. We do not yet know if the same is true of short-term emission variability, as the rotational changes involved cannot be directly measured over such short timescales. To remedy this, we propose a continuous pulsar monitoring technique that permits the statistical detection of any rotational changes in nulling and mode-changing pulsars with certain properties. Using a simulation, we explore the range of pulsar properties over which such an experiment would be possible., Comment: 11 pages, 8 figures, 2 tables