Researcher profile

B. W. Stappers

B. W. Stappers contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
37works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2026arXiv

Evaluating the effectiveness of radio frequency interference removal algorithms for single pulse searches

Radio Frequency Interference (RFI), the presence of artificial and/or terrestrial signals in astronomical data, poses a great challenge to the search for pulsars and radio transients, such as Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs), by obscuring or distorting the signal of interest and resulting in large numbers of erroneous detections. RFI mitigation algorithms aim to remove this interference and improve the chance of detection of transients, but with the growing number of techniques, selecting the most appropriate method for a given survey can be problematic. The choice of method is particularly important in real-time searches planned for next-generation telescopes such as those of the SKAO, where there is no possibility to reprocess the data. In this paper, we explore the algorithm selection problem by injecting pulses into data which simulates several RFI environments. A set of these files is then cleaned using RFI mitigation algorithms and run through a single pulse search pipeline to analyse the recovery of the injected pulses. We examine the recovery of the injected single pulses with an emphasis on a number of cases spanning a range of pulse brightness, width and dispersion measure. The efficacy and side effects of a few popular RFI excision methods, namely IQRM, SKF, and ZDMF are evaluated.

preprint2026arXiv

Radio timing constraints on the orbital orientation and component masses of PSR J1455$-$3330

PSR J1455$-$3330 is a $\sim$7.98 ms pulsar in a $\sim$76.17 day nearly circular orbit with a white dwarf companion. In this work, we combine the available Lovell, Nançay decimetric Radio Telescope, Green Bank, and MeerKAT pulsar timing data spanning $\sim$ 30 years to measure the kinematic and relativistic effects of PSR J1455$-$3330 to constrain its 3D orbital geometry and component masses. We detect a relativistic Shapiro delay signal. We measure a significant orthometric amplitude $h_3 = 0.307^{+0.022}_{-0.026}$ $μ$s and an orthometric ratio $ς= 0.551^{+0.057}_{-0.054}$. We measure the change in projected semi-major axis $\dot{x} = -202.1^{+2.5}_{-2.7} \times10^{-16} \, \rm s\,s^{-1}$ with high significance, parallax, $\varpi$ = 1.11(6) mas, parallax derived distance 0.90(5) kpc, and a precise total proper motion magnitude of 12.432(2) mas yr$^{-1}$. A self-consistent analysis of all kinematic and relativistic effects, assuming general relativity, yields two solutions: (1) a pulsar mass $M_{\rm p} = 1.39^{+0.38}_{-0.18}\, \rm M_{\odot}$, a companion mass $M_{\rm c} = 0.293^{+0.056}_{-0.026}$ $\rm M_{\odot}$, an orbital inclination, $i = 63(2)^{\circ}$, and longitude of the ascending node, $Ω= 212(12)^{\circ}$ or (2) a pulsar mass $M_{\rm p} = 1.53^{+1.10}_{-0.22} \, \rm M_{\odot}$, a companion mass $M_{\rm c} = 0.309^{+0.163}_{-0.026}\, \rm M_{\odot}$, an orbital inclination, $i = 123(4)^{\circ}$, and longitude of the ascending node, $Ω= 334(12)^{\circ}$. All uncertainties represent the 68.27$\%$ credibility region. These results strongly favour a helium-dominated white dwarf companion.

preprint2022arXiv

21 new long-term variables in the GX 339-4 field: two years of MeerKAT monitoring

We present 21 new long-term variable radio sources found commensally in two years of weekly MeerKAT monitoring of the low-mass X-ray binary GX 339-4. The new sources vary on time scales of weeks to months and have a variety of light curve shapes and spectral index properties. Three of the new variable sources are coincident with multi-wavelength counterparts; and one of these is coincident with an optical source in deep MeerLICHT images. For most sources, we cannot eliminate refractive scintillation of active galactic nuclei as the cause of the variability. These new variable sources represent $2.2\pm0.5$ per cent of the unresolved sources in the field, which is consistent with the 1-2 per cent variability found in past radio variability surveys. However, we expect to find short-term variable sources in the field as well as these 21 new long-term variable sources. We present the radio light curves and spectral index variability of the new variable sources, as well as the absolute astrometry and matches to coincident sources at other wavelengths.

preprint2022arXiv

Detection of quasi-periodic micro-structure in three millisecond pulsars with the Large European Array for Pulsars

We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology.

preprint2022arXiv

Discoveries and Timing of Pulsars in NGC 6440

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) &#34;Black Widow&#34;, 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 &#34;Redback&#34; pulsar NGC 6440D at two different frequency bands reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies (962-1283 MHz).

preprint2022arXiv

First discoveries and localisations of Fast Radio Bursts with MeerTRAP: a real-time, commensal MeerKAT survey

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&#39;$ 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&#39;s measured properties are consistent with other FRB hosts with secure associations.

preprint2022arXiv

Four pulsar discoveries in NGC 6624 by TRAPUM using MeerKAT

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.

preprint2022arXiv

Irradiated but not eclipsed, the case of PSR J0610-2100

We report on radio timing observations of the black widow binary pulsar J0610$-$2100 and optical observations of its binary companion. The radio timing observations extend the timing baseline to 16yr and reveal a marginal detection of the orbital period derivative, but they show no significant evidence of orbital variations such as those seen in other black widow pulsars. Furthermore, no eclipses are seen in the observations at observing frequencies ranging from 310 to 2700MHz. The optical VRI light curves were modulated with the orbital period, reaching maximum brightness of $V=26.8$, $R=25.4$, and $I=23.8$ at superior conjunction of the companion, confirming irradiation of the companion by the pulsar. Modelling the light curves indicates that the companion is likely not filling its Roche lobe, while having a moderate inclination ($i > 54\degr$). We find an unusually low temperature and a low irradiation for the irradiated hemisphere of the companion. We investigate the absence of radio eclipses in PSR J0610-2100 and in other black widow systems in relation to their binary, pulsar, and companion properties. We also discuss the suitability of PSR J0610-2100 for pulsar timing array observations aimed at detecting nano-Hertz gravitational waves.

preprint2022arXiv

Long term radio and X-ray evolution of the magnetar Swift J1818.0-1607

We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift~J1818.0-1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of $-$1.37$\times$10$^{-11}$ Hz s$^{-1}$. The magnitude of the variation in the frequency derivative reduces with time before converging on the mean value. We were able to identify four states in the spin-frequency derivative that were quantified by the amount of modulation about the mean value and the transition between these states seem to be correlated with the change in the radio emission of the magnetar while no correlation is seen in the average radio profile variability on a shorter timescale (days). The 0.5--12 keV X-ray flux shows a monotonic decrease that can be attributed to thermal emission from a hot-spot on the surface of the neutron star that is reducing in size. Such decrease is consistent with what is seen in other magnetars. The potential correlation between the radio emission mode and the behaviour of the spin-down rate hints to a global change in the magnetopshere of the magnetar akin to the correlation seen in a subset of mode-changing radio pulsars and suggests a physical link between the two sub-populations.

preprint2022arXiv

Long-term rotational and emission variability of 17 radio pulsars

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ν$) 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ν$ transitions ($Δ\dotν$) 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 $Δ\dotν$ amplitude on $\dotν$ and pulsar characteristic age. Linking $Δ\dotν$ to changes in the global magnetospheric charge density $Δρ$, we speculate that $\dotν$ transitions associated with large $Δρ$ values may be exhibiting detectable profile changes with improved data quality, in cases where they have not previously been observed.

preprint2022arXiv

MeerTRAP: Twelve Galactic fast transients detected in a real-time, commensal MeerKAT survey

MeerTRAP is a real-time untargeted search project using the MeerKAT telescope to find single pulses from fast radio transients and pulsars. It is performed commensally with the MeerKAT large survey projects (LSPs), using data from up to 64 of MeerKAT&#39;s 13.96~m dishes to form hundreds of coherent beams on sky, each of which is processed in real time to search for millisecond-duration pulses. We present the first twelve Galactic sources discovered by MeerTRAP, with DMs in the range of 33--381~pc~cm$^{-3}$. One source may be Galactic or extragalactic depending on the Galactic electron density model assumed. Follow-up observations performed with the MeerKAT, Lovell, and Parkes radio telescopes have detected repeat pulses from seven of the twelve sources. Pulse periods have been determined for four sources. Another four sources could be localised to the arcsecond-level using a novel implementation of the tied-array beam localisation method.

preprint2022arXiv

Modelling annual scintillation arc variations in PSR J1643-1224 using the Large European Array for Pulsars

In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised interstellar medium (IISM), and Earth. We observe a clear parabolic scintillation arc which varies in curvature throughout the year. The distribution of power in the secondary spectra are inconsistent with a single scattering screen which is fully 1D, or entirely isotropic. We fit the observed arc curvature variations with two models; an isotropic scattering screen, and a model with two independent 1D screens. We measure the distance to the scattering screen to be in the range 114-223 pc, depending on the model, consistent with the known distance of the foreground large-diameter HII region Sh 2-27 (112+/-17 pc), suggesting that it is the dominant source of scattering. We obtain only weak constraints on the pulsar&#39;s orbital inclination and angle of periastron, since the scintillation pattern is not very sensitive to the pulsar&#39;s motion, since the screen is much closer to the Earth than the pulsar. More measurements of this kind - where scattering screens can be associated with foreground objects - will help to inform the origins and distribution of scattering screens within our galaxy.

preprint2022arXiv

Multi-colour optical light curves of the companion star to the millisecond pulsar PSR J2051-0827

We present simultaneous, multi-colour optical light curves of the companion star to the black-widow pulsar PSR J2051-0827, obtained approximately 10 years apart using ULTRACAM and HiPERCAM, respectively. The ULTRACAM light curves confirm the previously reported asymmetry in which the leading hemisphere of the companion star appears to be brighter than the trailing hemisphere. The HiPERCAM light curves, however, do not show this asymmetry, demonstrating that whatever mechanism is responsible for it varies on timescales of a decade or less. We fit the symmetrical HiPERCAM light curves with a direct-heating model to derive the system parameters, finding an orbital inclination of $55.9^{+4.8}_{-4.1}$ degrees, in good agreement with radio-eclipse constraints. We find that approximately half of the pulsar&#39;s spin-down energy is converted to optical luminosity, resulting in temperatures ranging from approximately $5150^{+190}_{-190}$ K on the day side to $2750^{+130}_{-150}$ K on the night side of the companion star. The companion star is close to filling its Roche lobe ($f_{\rm RL} =0.88^{+0.02}_{-0.02}$) and has a mass of $0.039^{+0.010}_{-0.011}$ M$_{\odot}$, giving a mean density of $20.24^{+0.59}_{-0.44}$ g cm$^{-3}$ and an apsidal motion constant in the range $0.0036 < k_2 < 0.0047$. The companion mass and mean density values are consistent with those of brown dwarfs, but the apsidal motion constant implies a significantly more centrally-condensed internal structure than is typical for such objects.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo&#39;s third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

News and views regarding PSR J1757-1854, a highly-relativistic binary pulsar

We provide an update on the ongoing monitoring and study of the highly-relativistic double neutron star binary system PSR J1757-1854, a 21.5-ms pulsar in a highly eccentric, 4.4-hour orbit. The extreme nature of this pulsar&#39;s orbit allows it to probe a parameter space largely unexplored by other relativistic binary pulsars. For example, it displays one of the highest gravitational wave (GW) luminosities of any known binary pulsar, as well as the highest rate or orbital decay due to GW damping. PSR J1757-1854 is also notable in that it is an excellent candidate for exploring new tests of General Relativity and other gravitational theories, with possible measurements of both Lense-Thirring precession and relativistic orbital deformation (through the post-Keplerian parameter $δ_θ$) anticipated within the next 3-5 years. Here we present a summary of the latest interim results from the ongoing monitoring of this pulsar as part of an international, multi-telescope campaign. This includes an update of the pulsar&#39;s long-term timing and post-Keplerian parameters, new constraints on the pulsar&#39;s proper motion and corresponding Shklovskii kinematic correction, and new limits on the pulsar&#39;s geodetic precession as determined by monitoring for secular changes in the pulse profile. We also highlight prospects for future work, including an updated timeline on new relativistic tests following the introduction of MeerKAT observations.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars&#39; mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2022arXiv

Searching for pulsars associated with polarised point sources using LOFAR: Initial discoveries from the TULIPP project

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.

preprint2022arXiv

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

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

preprint2022arXiv

TRAPUM discovery of thirteen new pulsars in NGC 1851 using MeerKAT

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.

preprint2021arXiv

A comparative analysis of pulse time-of-arrival creation methods

Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived from total intensity pulse profiles of pulsars via cross-correlation with an idealised `1D&#39; template of that profile. While a number of results have been presented in the literature relying on the ever increasing sensitivity of such pulsar timing experiments, there is no consensus on the most reliable methods for TOA creation and, more importantly, the associated TOA uncertainties for each scheme. In this article, we present a comprehensive comparison of TOA determination practices, focusing on the creation of timing templates, TOA determination methods and the most useful TOA bandwidth. The aim is both to present a possible approach towards TOA optimisation as well as the (partial) identification of an optimal TOA-creation scheme and the demonstration of optimisation differences between pulsars and data sets. We compare the values of data-derived template profiles as compared to analytic profiles and evaluate the three most commonly used template-matching methods. Finally, we study the relation between timing precision and TOA bandwidth to identify any potential breaks in that relationship. As a practical demonstration, we apply our selected methods to European Pulsar Timing Array data on the three test pulsars PSRs\ J0218+4232, J1713+0747 and J2145$-$0750.

preprint2021arXiv

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst.

preprint2021arXiv

Apertif, Phased Array Feeds for the Westerbork Synthesis Radio Telescope

We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics.

preprint2021arXiv

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age $\sim$1 Gy), and PSR J1954+2529, which likely belongs to a newly-emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a non-recycled neutron star with a 931-ms spin period in an eccentric ($e\,=\,0.114$) wide ($P_b\,=\,82.7\,$d) orbit with a companion of undetermined nature having a minimum mass of $\sim0.6\,M_{\odot}$. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSP) have twice the dispersion measure per unit spin period than the known population of MSP in the Plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, RRATS, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought.

preprint2020arXiv

Discovery and timing of pulsars in the globular cluster M13 with FAST

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

preprint2020arXiv

Giant pulses from J1823-3021A observed with the MeerKAT telescope

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.

preprint2020arXiv

Limits on Absorption from a 332-MHz survey for Fast Radio Bursts

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

preprint2020arXiv

Optimal periodicity searching: Revisiting the Fast Folding Algorithm for large-scale pulsar surveys

The Fast Folding Algorithm (FFA) is a phase-coherent search technique for periodic signals. It has rarely been used in radio pulsar searches, having been historically supplanted by the less computationally expensive Fast Fourier Transform (FFT) with incoherent harmonic summing (IHS). Here we derive from first principles that an FFA search closely approaches the theoretical optimum sensitivity to all periodic signals; it is analytically shown to be significantly more sensitive than the standard FFT+IHS method, regardless of pulse period and duty cycle. A portion of the pulsar phase space has thus been systematically under-explored for decades; pulsar surveys aiming to fully sample the pulsar population should include an FFA search as part of their data analysis. We have developed an FFA software package, riptide, fast enough to process radio observations on a large scale; riptide has already discovered sources undetectable using existing FFT+IHS implementations. Our sensitivity comparison between search techniques also shows that a more realistic radiometer equation is needed, which includes an additional term: the search efficiency. We derive the theoretical efficiencies of both the FFA and the FFT+IHS methods and discuss how excluding this term has consequences for pulsar population synthesis studies.

preprint2020arXiv

Possible periodic activity in the repeating FRB 121102

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&#39; and that it should turn `on&#39; 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.

preprint2020arXiv

Simultaneous multi-telescope observations of FRB 121102

We present 11 detections of FRB 121102 in ~3 hours of observations during its &#39;active&#39; 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&#39;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 &#39;precursors&#39; 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 &#39;inactive&#39; over a broad frequency range during this time.

preprint2020arXiv

Study of spider pulsar binary eclipses and discovery of an eclipse mechanism transition

We present a comparative study of the low-frequency eclipses of spider (compact, irradiating binary) PSRs B1957+20 and J1816+4510. Combining these data with those of three other eclipsing systems we study the frequency dependence of the eclipse duration. PSRs B1957+20 and J1816+4510 have similar orbital properties, but the companions to the pulsars have masses that differ by an order of magnitude. A dedicated campaign to simultaneously observe the pulsed and imaged continuum flux densities throughout the eclipses reveals many similarities between the excess material within the two binaries, irrespective of the companion star properties. The observations show that the pulsar fluxes are removed from the line of sight throughout the main body of the eclipses. For PSR J1816+4510 we present the first direct evidence of an eclipse mechanism that transitions from one that removes the pulsar flux from the line of sight to one that merely smears out pulsations, and claim that this is a consequence of scattering in a tail of material flowing behind the companion. Inferred mass loss rates from the companion stars are found to be $\dot{M}_{\text{C}} \sim 10^{-12}~M_\odot$~yr$^{-1}$ and $\dot{M}_{\text{C}} \sim 2 \times 10^{-13}~M_\odot$~yr$^{-1}$ for PSR B1957+20 and PSR J1816+4510, respectively; seemingly too low to evaporate the stars within Hubble time. Measurements of eclipse durations over a wide range of radio-frequencies show a significant dependence of eclipse duration on frequency for all pulsars, with wider eclipses at lower-frequencies. These results provide a marked improvement in the observational constraints available for theoretical studies of the eclipse mechanisms.

preprint2020arXiv

The High Time Resolution Universe Pulsar Survey -- XVI. Discovery and timing of 40 pulsars from the southern Galactic plane

We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls.

preprint2020arXiv

The LOFAR Tied-Array All-Sky Survey: Timing of 21 pulsars including the first binary pulsar discovered with LOFAR

We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known.

preprint2020arXiv

The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime

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.

preprint2020arXiv

The radio pulsar population of the Small Magellanic Cloud

We model the present day, observable, normal radio pulsar population of the Small Magellanic Cloud (SMC). The pulsars are generated with SeBa, a binary population synthesis code that evolves binaries and the constituent stellar objects up to remnant formation and beyond. We define radio pulsars by selecting neutron stars that satisfy a selection of criteria defined by Galactic pulsars, and apply the detection thresholds of previous and future SMC pulsar surveys.The number of synthesised and recovered pulsars are exceptionally sensitive to the assumed star formation history and applied radio luminosity model, but is not affected extensively by the assumed common envelope model, metallicity, and neutron star kick velocity distribution. We estimate that the SMC formed (1.6$\pm$0.3)$\times 10^4$ normal pulsars during the last 100 Myrs. We study which pulsars could have been observed by the Parkes multibeam survey of the SMC, by applying the survey&#39;s specific selection effects, and recover 4.0$\pm$0.8 synthetic pulsars.This is in agreement with their five observed pulsars. We also apply a proposed MeerKAT configuration for the upcoming SMC survey, and predict that the MeerKAT survey will detect 17.2$\pm$2.5 pulsars.

preprint2020arXiv

The SUrvey for Pulsars and Extragalactic Radio Bursts IV: Discovery and polarimetry of a 12.1-second radio pulsar

We report the discovery of PSR~J2251$-$3711, a radio pulsar with a spin period of 12.1 seconds, the second longest currently known. Its timing parameters imply a characteristic age of 15 Myr, a surface magnetic field of $1.3 \times 10^{13}$~G and a spin-down luminosity of $2.9 \times 10^{29}~\mathrm{erg~s}^{-1}$. Its dispersion measure of 12.12(1)~$\mathrm{pc}~\mathrm{cm}^{-3}$ leads to distance estimates of 0.5 and 1.3 kpc according to the NE2001 and YMW16 Galactic free electron density models, respectively. Some of its single pulses show an uninterrupted 180 degree sweep of the phase-resolved polarization position angle, with an S-shape reminiscent of the rotating vector model prediction. However, the fact that this sweep occurs at different phases from one pulse to another is remarkable and without straightforward explanation. Although PSR~J2251$-$3711 lies in the region of the $P-\dot{P}$ parameter space occupied by the X-ray Isolated Neutron Stars (XINS), there is no evidence for an X-ray counterpart in our Swift XRT observation; this places a 99\%-confidence upper bound on its unabsorbed bolometric thermal luminosity of $1.1 \times 10^{31}~(d / 1~\mathrm{kpc})^2~\mathrm{erg/s}$ for an assumed temperature of 85 eV, where $d$ is the distance to the pulsar. Further observations are needed to determine whether it is a rotation-powered pulsar with a true age of at least several Myr, or a much younger object such as an XINS or a recently cooled magnetar. Extreme specimens like PSR J2251$-$3711 help bridge populations in the so-called neutron star zoo in an attempt to understand their origins and evolution.

preprint2019arXiv

A LOFAR census of non-recycled pulsars: extending below 80 MHz

We present the results from the low-frequency (40--78 MHz) extension of the first LOFAR pulsar census of non-recycled pulsars. We have used the Low-Band Antennas of the LOFAR core stations to observe 87 pulsars out of 158 that have been detected previously with the High-Band Antennas. Forty-three pulsars have been detected and we present here their flux densities and flux-calibrated profiles. Seventeen of these pulsars have not been, to our knowledge, detected before at such low frequencies. We re-calculate the spectral indices using the new low-frequency flux density measurements from the LOFAR census and discuss the prospects of studying pulsars at the very low frequencies with the current and upcoming facilities, such as NenuFAR.

preprint2019arXiv

A pulsar-based timescale from the International Pulsar Timing Array

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.