Source author record

M. Bailes

M. Bailes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

86works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

86 published item(s)

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.

preprint2024arXiv

An insight into chromatic behaviour of jitter in pulsars and its modelling: A case study of PSR J0437$-$4715

Pulse-to-pulse profile shape variations introduce correlations in pulsar times of arrival (TOAs) across radio frequency measured at the same observational epoch. This leads to a broadband noise in excess of radiometer noise, which is termed pulse jitter noise. The presence of jitter noise limits the achievable timing precision and decreases the sensitivity of pulsar-timing data sets to signals of interest such as nanohertz-frequency gravitational waves. Current white noise models used in pulsar timing analyses attempt to account for this, assuming complete correlation of uncertainties through the arrival times collected in a unique observation and no frequency dependence of jitter (which corresponds to a rank-one covariance matrix). However, previous studies show that the brightest millisecond pulsar at decimetre wavelengths, PSR J0437$-$4715, shows decorrelation and frequency dependence of jitter noise. Here we present a detailed study of the decorrelation of jitter noise in PSR J0437$-$4715 and implement a new technique to model it. We show that the rate of decorrelation due to jitter can be expressed as a power-law in frequency. We analyse the covariance matrix associated with the jitter noise process and find that a higher-rank-approximation is essential to account for the decorrelation and to account for frequency dependence of jitter noise. We show that the use of this novel method significantly improves the estimation of other chromatic noise parameters such as dispersion measure variations. However, we find no significant improvement in errors and estimation of other timing model parameters suggesting that current methods are not biased for other parameters, for this pulsar due to this misspecification. We show that pulse energy variations show a similar decorrelation to the jitter noise, indicating a common origin for both observables.

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) "Black Widow", with a very low mass companion (M$_{\rm c}$ > 0.006 M$_{\odot}$). It joins the other binary pulsars discovered so far in this cluster which all have low companion masses (M$_{\rm c}$ < 0.30 M$_{\odot}$). We present the results of long-term timing solutions obtained using data from both Green Bank and MeerKAT telescopes for these two new pulsars and an analysis of the pulsars NGC 6440C and NGC 6440D. For the isolated pulsar NGC 6440C, we searched for planets using a Markov Chain Monte Carlo technique. We find evidence for significant unmodelled variations but they cannot be well modelled as planets nor as part of a power-law red-noise process. Studies of the eclipses of the "Redback" pulsar NGC 6440D at two different frequency bands reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies (962-1283 MHz).

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'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'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's long-term timing and post-Keplerian parameters, new constraints on the pulsar's proper motion and corresponding Shklovskii kinematic correction, and new limits on the pulsar'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' 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

Systematic upper limits on the size of missing pulsar glitches in the first UTMOST open data release

A systematic, semi-automated search for pulsar glitches in the first UTMOST public data release is presented. The search is carried out using a hidden Markov model which incorporates both glitches and timing noise into the model of the assumed phase evolution of the pulsar. Glitches are detected through Bayesian model selection between models with and without glitches present with minimal human intervention. Nine glitches are detected among seven objects, all of which have been previously reported. No new glitches were detected. Injection studies are used to place 90\% frequentist upper limits on the size of undetected glitches in each of the 282 objects searched. The mean upper limit obtained is $Δf^{90\%}/f = 1.9 \times 10^{-8}$, with a range of $4.1 \times 10^{-11} \leq Δf^{90\%}/f \leq 2.7 \times 10^{-7}$, assuming step events with no post-glitch recoveries. It is demonstrated that including glitch recovery has a mild effect, in most cases increasing the upper limit by a factor of $\lesssim 5$ conservatively assuming complete recovery on a timescale of $100\,\mathrm{d}$.

preprint2022arXiv

The eccentric millisecond pulsar, PSR J0955$-$6150 I: Pulse profile analysis, mass measurements and constraints on binary evolution

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 ($α$) 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ω = 0.00152(1) \, {\rm deg} \, yr^{-1}$ and the orthometric Shapiro delay parameters, $h_3 = 0.89(7) \, μ$s and $ς= 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.

preprint2022arXiv

The High Time Resolution Universe Pulsar Survey -- XVII. PSR J1325-6253, a low eccentricity double neutron star system from an ultra-stripped supernova

The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surprisingly low eccentricity of just $e = 0.064$. Through 1.4 yr of dedicated timing with the Parkes radio telescope we have been able to measure its rate of advance of periastron, $\dotω=0.138 \pm 0.002$ $\rm deg$ $\rm yr^{-1}$. If this induced $\dotω$ is solely due to general relativity then the total mass of the system is, $M_{\rm sys} = 2.57 \pm 0.06$ M$_{\odot}$. Assuming an edge-on orbit the minimum companion mass is constrained to be $M_\mathrm{c,min}>0.98$ M$_{\odot}$ which implies the pulsar mass is $M_\mathrm{p,max}<1.59 $ M$_{\odot}$. Its location in the $P$-$\dot{P}$ diagram suggests that, like other DNS systems, PSR J1325$-$6253 is a recycled pulsar and if its mass is similar to the known examples ($>1.3$ M$_\odot$), then the companion neutron star is probably less than $\sim1.25$ M$_\odot$ and the system is inclined at about $50^{\circ}$-$60^{\circ}$. The low eccentricity along with the wide orbit of the system strongly favours a formation scenario involving an ultra-stripped supernova explosion.

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

The MeerTime Pulsar Timing Array -- A Census of Emission Properties and Timing Potential

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 $μ$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 $μ$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.

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

Measurements of pulse jitter and single-pulse variability in millisecond pulsars using MeerKAT

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

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

Lense-Thirring frame dragging induced by a fast-rotating white dwarf in a binary pulsar system

Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141$-$6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observe a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of a Newtonian quadrupole moment and Lense-Thirring precession of the orbit resulting from rapid rotation of the WD. Lense-Thirring precession, an effect of relativistic frame-dragging, is a prediction of general relativity. This detection is consistent with the evolutionary scenario in which the WD accreted matter from the pulsar progenitor, spinning up the WD to a period $< 200$ seconds.

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 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 Parkes Pulsar Timing Array Project: Second data release

We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 $μ$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files.

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.

preprint2020arXiv

The SUrvey for Pulsars and Extragalactic Radio Bursts V: Recent Discoveries and Full Timing Solutions

The SUrvey for Pulsars and Extragalactic Radio Bursts ran from 2014 April to 2019 August, covering a large fraction of the southern hemisphere at mid- to high-galactic latitudes, and consisting of 9-minute pointings taken with the 20-cm multibeam receiver on the Parkes Radio Telescope. Data up to 2017 September 21 have been searched using standard Fourier techniques, single-pulse searches, and Fast Folding Algorithm searches. We present 19 new discoveries, bringing the total to 27 discoveries in the programme, and we report the results of follow-up timing observations at Parkes for 26 of these pulsars, including the millisecond pulsar PSR J1421-4409; the faint, highly-modulated, slow pulsar PSR J1646-1910; and the nulling pulsar PSR J1337-4441. We present new timing solutions for 23 pulsars, and we report flux densities, modulation indices, and polarization properties.

preprint2020arXiv

The Thousand-Pulsar-Array programme on MeerKAT I: Science objectives and first results

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.

preprint2020arXiv

Timing of young radio pulsars II. Braking indices and their interpretation

In Paper I of this series, we detected a significant value of the braking index ($n$) for 19 young, high-$\dot{E}$ radio pulsars using $\sim$ 10 years of timing observations from the 64-m Parkes radio telescope. Here we investigate this result in more detail using a Bayesian pulsar timing framework to model timing noise and to perform selection to distinguish between models containing exponential glitch recovery and braking index signatures. We show that consistent values of $n$ are maintained with the addition of substantial archival data, even in the presence of glitches. We provide strong arguments that our measurements are unlikely due to exponential recovery signals from unseen glitches even though glitches play a key role in the evolution of a pulsar's spin frequency. We conclude that, at least over decadal time scales, the value of $n$ can be significantly larger than the canonical 3 and discuss the implications for the evolution of pulsars.

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.

preprint2019arXiv

An ultra-wide bandwidth (704 to 4032 MHz) receiver for the Parkes radio telescope

We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.

preprint2019arXiv

The UTMOST Survey for Magnetars, Intermittent pulsars, RRATs and FRBs I: System description and overview

We describe the ongoing `Survey for Magnetars, Intermittent pulsars, Rotating radio transients and Fast radio bursts' (SMIRF), performed using the newly refurbished UTMOST telescope. SMIRF repeatedly sweeps the southern Galactic plane performing real-time periodicity and single-pulse searches, and is the first survey of its kind carried out with an interferometer. SMIRF is facilitated by a robotic scheduler which is capable of fully autonomous commensal operations. We report on the SMIRF observational parameters, the data analysis methods, the survey's sensitivities to pulsars, techniques to mitigate radio frequency interference and present some early survey results. UTMOST's wide field of view permits a full sweep of the Galactic plane to be performed every fortnight, two orders of magnitude faster than previous surveys. In the six months of operations from January to June 2018, we have performed $\sim 10$ sweeps of the Galactic plane with SMIRF. Notable blind re-detections include the magnetar PSR J1622$-$4950, the RRAT PSR J0941$-$3942 and the eclipsing pulsar PSR J1748$-$2446A. We also report the discovery of a new pulsar, PSR J1705$-$54. Our follow-up of this pulsar with the UTMOST and Parkes telescopes at an average flux limit of $\leq 20$ mJy and $\leq 0.16$ mJy respectively, categorizes this as an intermittent pulsar with a high nulling fraction of $< 0.002$

preprint2016arXiv

A Fast Radio Burst Host Galaxy

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $Ω_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

preprint2016arXiv

Are the distributions of Fast Radio Burst properties consistent with a cosmological population?

High time resolution radio surveys over the last few years have discovered a population of millisecond-duration transient bursts called Fast Radio Bursts (FRBs), which remain of unknown origin. FRBs exhibit dispersion consistent with propagation through a cold plasma and dispersion measures indicative of an origin at cosmological distances. In this paper we perform Monte Carlo simulations of a cosmological population of FRBs, based on assumptions consistent with observations of their energy distribution, their spatial density as a function of redshift and the properties of the interstellar and intergalactic media. We examine whether the dispersion measures, fluences, inferred redshifts, signal-to-noises and effective widths of known FRBs are consistent with a cosmological population. Statistical analyses indicate that at least 50 events at Parkes are required to distinguish between a constant co-moving FRB density, and a FRB density that evolves with redshift like the cosmological star formation rate density.

preprint2016arXiv

Fast Radio Transient searches with UTMOST at 843 MHz

We report the first radio interferometric search at 843 MHz for fast transients, particularly Fast Radio Bursts (FRBs). The recently recommissioned Swinburne University of Technology's digital backend for the Molonglo Observatory Synthesis Telescope array (the UTMOST) with its large collecting area (18,000 $\mathrm{m^2}$) and wide instantaneous field of view (7.80 $\mathrm{deg^2}$) is expected to be an efficient tool to detect FRBs. As an interferometer it will be capable of discerning whether the FRBs are truly a celestial population. We show that UTMOST at full design sensitivity can detect an event approximately every few days. We report on 2 preliminary FRB surveys at about 7% and 14% respectively of the array's final sensitivity. Several pulsars have been detected via single pulses and no FRBs were discovered with pulse widths ($W$), in the range 655.36 $μ$s $< W < 41.9$ ms and dispersion measures (DMs) in the range $100 < $DM$< 2000$ $\mathrm{pc\,cm^{-3}}$. This non-detection sets a 2$σ$ upper limit of the sky rate of not more than 1000 events $\mathrm{sky^{-1}}$ $\mathrm{day^{-1}}$ at 843 MHz down to a flux limit of 11 Jy for 1 ms FRBs. We show that this limit is consistent with previous survey limits at 1.4 GHz and 145 MHz and set a lower limit on the mean spectral index of FRBs of $α> -3.2$.

preprint2016arXiv

FRBCAT: The Fast Radio Burst Catalogue

Here we present a catalogue of known Fast Radio Burst (FRB) sources in the form of an online catalogue, FRBCAT. The catalogue includes information about the instrumentation used for the observations for each detected burst, the measured quantities from each observation, and model-dependent quantities derived from observed quantities. To aid in consistent comparisons of burst properties such as width and signal-to-noise ratios we have reprocessed all the bursts for which we have access to the raw data, with software which we make available. The originally derived properties are also listed for comparison. The catalogue is hosted online as a MySQL database which can also be downloaded in tabular or plain text format for off-line use. This database will be maintained for use by the community for studies of the FRB population as it grows.

preprint2016arXiv

The disturbance of a millisecond pulsar magnetosphere

Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability.

preprint2016arXiv

The magnetic field and turbulence of the cosmic web measured using a brilliant fast radio burst

Fast radio bursts (FRBs) are millisecond-duration events thought to originate beyond the Milky Way galaxy. Uncertainty surrounding the burst sources, and their propagation through intervening plasma, has limited their use as cosmological probes. We report on a mildly dispersed (dispersion measure 266.5+-0.1 pc cm^-3), exceptionally intense (120+-30 Jy), linearly polarized, scintillating burst (FRB 150807) that we directly localize to 9 arcmin^2. Based on a low Faraday rotation (12.0+-0.7 rad m^-2), we infer negligible magnetization in the circum-burst plasma and constrain the net magnetization of the cosmic web along this sightline to <21 nG, parallel to the line-of-sight. The burst scintillation suggests weak turbulence in the ionized intergalactic medium.

preprint2015arXiv

A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $λ^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

preprint2015arXiv

A survey of FRB fields: Limits on repeatability

Several theories exist to explain the source of the bright, millisecond duration pulses known as fast radio bursts (FRBs). If the progenitors of FRBs are non-cataclysmic, such as giant pulses from pulsars, pulsar-planet binaries, or magnetar flares, FRB emission may be seen to repeat. We have undertaken a survey of the fields of eight known FRBs from the High Time Resolution Universe survey to search for repeating pulses. Although no repeat pulses were detected the survey yielded the detection of a new FRB, described in Petroff et al. (2015a). From our observations we rule out periodic repeating sources with periods P $\leq$ 8.6 hours and rule out sources with periods 8.6 < P < 21 hours at the 90% confidence level. At P $\geq$ 21 hours our limits fall off as ~1/P. Dedicated and persistent observations of FRB source fields are needed to rule out repetition on longer timescales, a task well-suited to next generation wide-field transient detectors.

preprint2015arXiv

Five new Fast Radio Bursts from the HTRU high latitude survey: first evidence for two-component bursts

The detection of five new fast radio bursts (FRBs) found in the High Time Resolution Universe high latitude survey is presented. The rate implied is 6$^{+4}_{-3}\times~10^3$ (95%) FRBs sky$^{-1}$ day$^{-1}$ above a fluence of between 0.13 and 5.9 Jy ms for FRBs between 0.128 and 262 ms in duration. One of these FRBs has a clear two-component profile, each component is similar to the known population of single component FRBs and are separated by 2.4(4) ms. All the FRB components appear to be unresolved following deconvolution with a scattering tail and accounting for intra-channel smearing. The two-component FRB also has the highest dispersion measure (1629 pc cm$^{-3}$) of any FRB to-date. Many of the proposed models to explain FRBs use a single high energy event involving compact objects (such as neutron star mergers) and therefore cannot easily explain a two-component FRB. Models that are based on extreme versions of flaring, pulsing or orbital events however could produce multiple component profiles. The compatibility of these models and the FRB rate implied by these detections is discussed.

preprint2015arXiv

Gravitational waves from binary supermassive black holes missing in pulsar observations

Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.

preprint2015arXiv

Pulsar Observations of Extreme Scattering Events

Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the electron density of order 10s of cm^{-3}. Similar structures have been observed in measurements of pulsar intensity scintillation and group delay. Here we report observations of two ESEs showing increases in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). These allow us to make more complete models of the ESE, including an estimate of the "outer-scale" of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully turbulent on an AU scale. They provide some support for the idea that the structures are extended along the line of sight, such as would be the case for a scattering shell. The dispersion measurements also show a variety of AU scale structures which would not be called ESEs, yet involve electron density variations typical of ESEs and likely have the same origin.

preprint2015arXiv

The High Time Resolution Universe Pulsar Survey XII : Galactic plane acceleration search and the discovery of 60 pulsars

We present initial results from the low-latitude Galactic plane region of the High Time Resolution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the computational challenges arising from the processing of the terabyte-sized survey data. Two new radio interference mitigation techniques are introduced, as well as a partially-coherent segmented acceleration search algorithm which aims to increase our chances of discovering highly-relativistic short-orbit binary systems, covering a parameter space including potential pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio of data length over orbital period of ~0.1 results in the highest effectiveness for this search algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with periods less than 30ms. PSR J1101-6424 is a millisecond pulsar whose heavy white dwarf (WD) companion and short spin period of 5.1ms indicate a rare example of full-recycling via Case A Roche lobe overflow. PSR J1757-27 appears to be an isolated recycled pulsar with a relatively long spin period of 17ms. In addition, PSR J1244-6359 is a mildly-recycled binary system with a heavy WD companion, PSR J1755-25 has a significant orbital eccentricity of 0.09, and PSR J1759-24 is likely to be a long-orbit eclipsing binary with orbital period of the order of tens of years. Comparison of our newly-discovered pulsar sample to the known population suggests that they belong to an older population. Furthermore, we demonstrate that our current pulsar detection yield is as expected from population synthesis.

preprint2015arXiv

Timing analysis for 20 millisecond pulsars in the Parkes Pulsar Timing Array

We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with $M_p=1.44\pm0.07$ $M_\odot$ and $M_p=1.47\pm0.03$ $M_\odot$ respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16% level of precision, $D=156.79\pm0.25$ pc, one of the most fractionally precise distance measurements of any star to date.

preprint2015arXiv

Versatile Directional Searches for Gravitational Waves with Pulsar Timing Arrays

By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "phased-up" to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise.

preprint2014arXiv

A real-time fast radio burst: polarization detection and multiwavelength follow-up

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$σ$) circularly polarized on the leading edge with a 1-$σ$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

preprint2014arXiv

An absence of fast radio bursts at intermediate galactic latitudes

Fast radio bursts (FRBs) are an emerging class of bright, highly dispersed radio pulses. Recent work by Thornton et al. (2013) has revealed a population of FRBs in the High Time Resolution Universe (HTRU) survey at high Galactic latitudes. A variety of progenitors have been proposed including cataclysmic events at cosmological distances, Galactic flare stars, and terrestrial radio frequency interference. Here we report on a search for FRBs at intermediate Galactic latitudes ($-15^{\circ}$ $< b <$ 15$^{\circ}$) in data taken as part of the HTRU survey. No FRBs were discovered in this region. Several effects such as dispersion, scattering, sky temperature and scintillation decrease the sensitivity by more than 3$σ$ in $\sim$20\% of survey pointings. Including all of these effects, we exclude the hypothesis that FRBs are uniformly distributed on the sky with 99\% confidence. This low probability implies that additional factors -- not accounted for by standard Galactic models -- must be included to ease the discrepancy between the detection rates at high and low Galactic latitudes. A revised rate estimate or another strong and heretofore unknown selection effect in Galactic latitude would provide closer agreement between the surveys' detection rates. The dearth of detections at low Galactic latitude disfavors a Galactic origin for these bursts.

preprint2014arXiv

An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set

We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude $h_0$ for a range of GW frequencies. For example, at $10^{-8}$ Hz our analysis has excluded with $95\%$ confidence the presence of signals with $h_0\geqslant 1.7\times 10^{-14}$. Our new limits are about a factor of four more stringent than those of Yardley et al. (2010) based on an earlier PPTA data set and a factor of two better than those reported in the recent Arzoumanian et al. (2014) paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of $10^{9} M_{\odot}$ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of $4 \times 10^{-3} {\rm{Mpc}}^{-3} {\rm{Gyr}}^{-1}$ at $95\%$ confidence on the coalescence rate of nearby ($z \lesssim 0.1$) supermassive binary black holes in circular orbits with chirp masses of $10^{10}M_{\odot}$.

preprint2014arXiv

Limitations in timing precision due to single-pulse shape variability in millisecond pulsars

High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.

preprint2014arXiv

Pulsars in Globular Clusters with the SKA

Globular clusters are highly efficient radio pulsar factories. These pulsars can be used as precision probes of the clusters' structure, gas content, magnetic field, and formation history; some of them are also highly interesting in their own right because they probe exotic stellar evolution scenarios as well as the physics of dense matter, accretion, and gravity. Deep searches with SKA1-MID and SKA1-LOW will plausibly double to triple the known population. Such searches will only require one to a few tied-array beams, and can be done during early commissioning of the telescope - before an all-sky pulsar survey using hundreds to thousands of tied-array beams is feasible. With SKA2 it will be possible to observe most of the active radio pulsars within a large fraction of the Galactic globular clusters, an estimated population of 600 - 3700 observable pulsars (those beamed towards us). This rivals the total population of millisecond pulsars that can be found in the Galactic field; fully characterizing it will provide the best-possible physical laboratories as well as a rich dynamical history of the Galactic globular cluster system.

preprint2014arXiv

Searching for gravitational wave memory bursts with the Parkes Pulsar Timing Array

Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.

preprint2014arXiv

SPINN: a straightforward machine learning solution to the pulsar candidate selection problem

We describe SPINN (Straightforward Pulsar Identification using Neural Networks), a high-performance machine learning solution developed to process increasingly large data outputs from pulsar surveys. SPINN has been cross-validated on candidates from the southern High Time Resolution Universe (HTRU) survey and shown to identify every known pulsar found in the survey data while maintaining a false positive rate of 0.64%. Furthermore, it ranks 99% of pulsars among the top 0.11% of candidates, and 95% among the top 0.01%. In conjunction with the PEASOUP pipeline (Barr et al., in prep.), it has already discovered four new pulsars in a re-processing of the intermediate Galactic latitude area of HTRU, three of which have spin periods shorter than 5 milliseconds. SPINN's ability to reduce the amount of candidates to visually inspect by up to four orders of magnitude makes it a very promising tool for future large-scale pulsar surveys. In an effort to provide a common testing ground for pulsar candidate selection tools and stimulate interest in their development, we also make publicly available the set of candidates on which SPINN was cross-validated.

preprint2014arXiv

The High Time Resolution Universe Pulsar Survey X: Discovery of four millisecond pulsars and updated timing solutions of a further 12

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

preprint2014arXiv

The High Time Resolution Universe Survey - XI. Discovery of five recycled pulsars and the optical detectability of survey white dwarf companions

We present the discovery of a further five recycled pulsar systems in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods ranging from 2 ms to 66 ms, and four are in binary systems with orbital periods between 10.8 hours and 9.0 days. Three of these binary systems are particularly interesting; PSR J1227-6208 has a pulse period of 34.5 ms and the highest mass function of all pulsars with near-circular orbits. The circular orbit suggests that the companion is not another neutron star, so future timing experiments may reveal one of the heaviest white dwarfs ever found ($>$ 1.3 M$_\odot$). Timing observations of PSR J1431$-$4715 indicate that it is eclipsed by its companion which has a mass indicating it belongs to the redback class of eclipsing millisecond pulsars. PSR J1653-2054 has a companion with a minimum mass of only $0.08$ M$_\odot$, placing it among the class of pulsars with low-mass companions. Unlike the majority of such systems, however, no evidence of eclipses is seen at 1.4 GHz.

preprint2014arXiv

Timing, polarimetry and physics of the bright, nearby millisecond pulsar PSR J0437-4715 - a single-pulse perspective

Single pulses from radio pulsars contain a wealth of information about emission and propagation in the magnetosphere and insight into their timing properties. It was recently demonstrated that single-pulse emission is responsible for limiting the timing stability of the brightest of millisecond pulsars. We report on an analysis of more than a million single-pulses from PSR J0437-4715 and present various statistical properties such as the signal-to-noise ratio (S/N) distribution, timing and polarimetry of average profiles integrated from subpulses with chosen S/N cut-offs, modulation properties of the emission, phase-resolved statistics of the S/N, and two dimensional spherical histograms of the polarization vector orientation. The last of these indicates the presence of orthogonally polarised modes (OPMs). Combined with the dependence of the polarisation fraction on the S/N and polarimetry of the brightest pulses, the existence of OPMs constrains pulsar emission mechanisms and models for the plasma physics in the magnetosphere.

preprint2013arXiv

A Population of Fast Radio Bursts at Cosmological Distances

Searches for transient astrophysical sources often reveal unexpected classes of objects that are useful physical laboratories. In a recent survey for pulsars and fast transients we have uncovered four millisecond-duration radio transients all more than 40° from the Galactic plane. The bursts' properties indicate that they are of celestial rather than terrestrial origin. Host galaxy and intergalactic medium models suggest that they have cosmological redshifts of 0.5 to 1, and distances of up to 3 gigaparsecs. No temporally coincident x- or gamma-ray signature was identified in association with the bursts. Characterization of the source population and identification of host galaxies offers an opportunity to determine the baryonic content of the Universe.

preprint2013arXiv

Detection of fast transients with radio interferometric arrays

Next-generation radio arrays, including the Square Kilometre Array (SKA) and its pathfinders, will open up new avenues for exciting transient science at radio wavelengths. Their innovative designs, comprising a large number of small elements, pose several challenges in digital processing and optimal observing strategies. The Giant Metrewave Radio Telescope (GMRT) presents an excellent test-bed for developing and validating suitable observing modes and strategies for transient experiments with future arrays. Here we describe the first phase of the ongoing development of a transient detection system for GMRT that is planned to eventually function in a commensal mode with other observing programs. It capitalizes on the GMRT's interferometric and sub-array capabilities, and the versatility of a new software backend. We outline considerations in the plan and design of transient exploration programs with interferometric arrays, and describe a pilot survey that was undertaken to aid in the development of algorithms and associated analysis software. This survey was conducted at 325 and 610 MHz, and covered 360 deg$^2$ of the sky with short dwell times. It provides large volumes of real data that can be used to test the efficacies of various algorithms and observing strategies applicable for transient detection. We present examples that illustrate the methodologies of detecting short-duration transients, including the use of sub-arrays for higher resilience to spurious events of terrestrial origin, localization of candidate events via imaging and the use of a phased array for improved signal detection and confirmation. In addition to demonstrating applications of interferometric arrays for fast transient exploration, our efforts mark important steps in the roadmap toward SKA-era science.

preprint2013arXiv

Gravitational-wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be $Ω_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.

preprint2013arXiv

The High Time Resolution Universe Pulsar Survey - VII: discovery of five millisecond pulsars and the different luminosity properties of binary and isolated recycled pulsars

This paper presents the discovery and timing parameters for five millisecond pulsars (MSPs), four in binary systems with probable white dwarf companions and one isolated, found in ongoing processing of the High Time Resolution Universe Pulsar Survey (HTRU). We also present high quality polarimetric data on four of them. These further discoveries confirm the high potential of our survey in finding pulsars with very short spin periods. At least two of these five MSPs are excellent candidates to be included in the Pulsar Timing Array projects. Thanks to the wealth of MSP discoveries in the HTRU survey, we revisit the question of whether the luminosity distributions of isolated and binary MSPs are different. Using the Cordes and Lazio distance model and our new and catalogue flux density measurements, we find that 41 of the 42 most luminous MSPs in the Galactic disk are in binaries and a statistical analysis suggests that the luminosity functions differ with 99.9% significance. We conclude that the formation process that leads to solitary MSPs affects their luminosities, despite their period and period derivatives being similar to those of pulsars in binary systems.

preprint2013arXiv

The High Time Resolution Universe Pulsar Survey VIII: The Galactic millisecond pulsar population

We have used millisecond pulsars (MSPs) from the southern High Time Resolution Universe (HTRU) intermediate latitude survey area to simulate the distribution and total population of MSPs in the Galaxy. Our model makes use of the scale factor method, which estimates the ratio of the total number of MSPs in the Galaxy to the known sample. Using our best fit value for the z-height, z=500 pc, we find an underlying population of MSPs of 8.3(\pm 4.2)*10^4 sources down to a limiting luminosity of L_min=0.1 mJy kpc^2 and a luminosity distribution with a steep slope of d\log N/d\log L = -1.45(\pm 0.14). However, at the low end of the luminosity distribution, the uncertainties introduced by small number statistics are large. By omitting very low luminosity pulsars, we find a Galactic population above L_min=0.2 mJy kpc^2 of only 3.0(\pm 0.7)*10^4 MSPs. We have also simulated pulsars with periods shorter than any known MSP, and estimate the maximum number of sub-MSPs in the Galaxy to be 7.8(\pm 5.0)*10^4 pulsars at L=0.1 mJy kpc^2. In addition, we estimate that the high and low latitude parts of the southern HTRU survey will detect 68 and 42 MSPs respectively, including 78 new discoveries. Pulsar luminosity, and hence flux density, is an important input parameter in the model. Some of the published flux densities for the pulsars in our sample do not agree with the observed flux densities from our data set, and we have instead calculated average luminosities from archival data from the Parkes Telescope. We found many luminosities to be very different than their catalogue values, leading to very different population estimates. Large variations in flux density highlight the importance of including scintillation effects in MSP population studies.

preprint2013arXiv

The High Time Resolution Universe Survey - IX: Polarimetry of long-period pulsars

We present a polarimetric analysis of 49 long-period pulsars discovered as part of the High Time Resolution Universe (HTRU) southern survey. The sources exhibit the typical characteristics of "old" pulsars, with low fractional linear and circular polarisation and narrow, multicomponent profiles. Although the position angle swings are generally complex, for two of the analysed pulsars (J1622-3751 and J1710-2616) we obtained an indication of the geometry via the rotating vector model. We were able to determine a value of the rotation measure (RM) for 34 of the sources which, when combined with their dispersion measures (DM), yields an integrated magnetic field strength along the line of sight. With the data presented here, the total number of values of RM associated to pulsars discovered during the HTRU southern survey sums to 51. The RMs are not consistent with the hypothesis of a counter-clockwise direction of the Galactic magnetic field within an annulus included between 4 and 6 kpc from the Galactic centre. A partial agreement with a counter-clockwise sense of the Galactic magnetic field within the spiral arms is however found in the area of the Carina-Sagittarius arm.

preprint2012arXiv

A Shapiro delay detection in the binary system hosting the millisecond pulsar PSR J1910-5959A

PSR J1910-5959A is a binary pulsar with a helium white dwarf companion located about 6 arcmin from the center of the globular cluster NGC6752. Based on 12 years of observations at the Parkes radio telescope, the relativistic Shapiro delay has been detected in this system. We obtain a companion mass Mc = 0.180+/-0.018Msun (1sigma) implying that the pulsar mass lies in the range 1.1Msun <= Mp <= 1.5Msun. We compare our results with previous optical determinations of the companion mass, and examine prospects for using this new measurement for calibrating the mass-radius relation for helium white dwarfs and for investigating their evolution in a pulsar binary system. Finally we examine the set of binary systems hosting a millisecond pulsar and a low mass helium white dwarf for which the mass of both stars has been measured. We confirm that the correlation between the companion mass and the orbital period predicted by Tauris & Savonije reproduces the observed values but find that the predicted Mp - Pb correlation over-estimates the neutron star mass by about 0.5Msun in the orbital period range covered by the observations. Moreover, a few systems do not obey the observed Mp - Pb correlation. We discuss these results in the framework of the mechanisms that inhibit the accretion of matter by a neutron star during its evolution in a low-mass X-ray binary.

preprint2012arXiv

Development of a pulsar-based timescale

Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.

preprint2012arXiv

Measurement and correction of variations in interstellar dispersion in high-precision pulsar timing

Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales. In this paper we demonstrate that without consideration of wavelength- independent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem. Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603-7202 of ~2x10^{-3} cm^{-3}pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.

preprint2012arXiv

PSR J1141-6545: a powerful laboratory of GR and tensor-scalar theories of gravity

Pulsars in close binary systems have provided some of the most stringent tests of strong-field gravity to date. The pulsar--white-dwarf binary system J1141-6545 is specifically interesting due to its gravitational asymmetry which makes it one of the most powerful probes of tensor-scalar theories of gravity. We give an overview of current gravitational tests provided by the J1141-6545 binary system and comment on how anomalous accelerations, geodetic precession and timing instabilities may be prevented from limiting future tests of gravity to come from this system.

preprint2012arXiv

Radio emission evolution, polarimetry and multifrequency single pulse analysis of the radio magnetar PSR J1622-4950

Here we report on observations of the radio magnetar PSR J1622-4950 at frequencies from 1.4 to 17 GHz. We show that although its flux density is varying up to a factor of ~10 within a few days, it has on average decreased by a factor of 2 over the last 700 days. At the same time, timing analysis indicates a trend of decreasing spin-down rate over our entire data set, again of about a factor of 2 over 700 days, but also an erratic variability in the spin-down rate within this time span. Integrated pulse profiles are often close to 100 per cent linearly polarized, but large variations in both the profile shape and fractional polarization are regularly observed. Furthermore, the behaviour of the position angle of the linear polarization is very complex - offsets in both the absolute position angle and the phase of the position angle sweep are often seen and the occasional presence of orthogonal mode jumps further complicates the picture. However, model fitting indicates that the magnetic and rotation axes are close to aligned. Finally, a single pulse analysis has been carried out at four observing frequencies, demonstrating that the wide pulse profile is built up of narrow spikes of emission, with widths that scale inversely with observing frequency. All three of the known radio magnetars seem to have similar characteristics, with highly polarized emission, time-variable flux density and pulse profiles, and with spectral indices close to zero.

preprint2012arXiv

The High Time Resolution Universe Survey - V: Single-pulse energetics and modulation properties of 315 pulsars

We report on the pulse-to-pulse energy distributions and phase-resolved modulation properties for catalogued pulsars in the southern High Time Resolution Universe intermediate-latitude survey. We selected the 315 pulsars detected in a single-pulse search of this survey, allowing a large sample unbiased regarding any rotational parameters of neutron stars. We found that the energy distribution of many pulsars is well-described by a log-normal distribution, with few deviating from a small range in log-normal scale and location parameters. Some pulsars exhibited multiple energy states corresponding to mode changes, and implying that some observed "nulling" may actually be a mode-change effect. PSRJ1900-2600 was found to emit weakly in its previously-identified "null" state. We found evidence for another state-change effect in two pulsars, which show bimodality in their nulling time scales; that is, they switch between a continuous-emission state and a single-pulse-emitting state. Large modulation occurs in many pulsars across the full integrated profile, with increased sporadic bursts at leading and trailing sub-beam edges. Some of these high-energy outbursts may indicate the presence of "giant pulse" phenomena. We found no correlation with modulation and pulsar period, age, or other parameters. Finally, the deviation of integrated pulse energy from its average value was generally quite small, despite the significant phase-resolved modulation in some pulsars; we interpret this as tenuous evidence of energy regulation between distinct pulsar sub-beams.

preprint2012arXiv

The High Time Resolution Universe Survey VI: An Artificial Neural Network and Timing of 75 Pulsars

We present 75 pulsars discovered in the mid-latitude portion of the High Time Resolution Universe survey, 54 of which have full timing solutions. All the pulsars have spin periods greater than 100 ms, and none of those with timing solutions are in binaries. Two display particularly interesting behaviour; PSR J1054-5944 is found to be an intermittent pulsar, and PSR J1809-0119 has glitched twice since its discovery. In the second half of the paper we discuss the development and application of an artificial neural network in the data-processing pipeline for the survey. We discuss the tests that were used to generate scores and find that our neural network was able to reject over 99% of the candidates produced in the data processing, and able to blindly detect 85% of pulsars. We suggest that improvements to the accuracy should be possible if further care is taken when training an artificial neural network; for example ensuring that a representative sample of the pulsar population is used during the training process, or the use of different artificial neural networks for the detection of different types of pulsars.

preprint2012arXiv

The Parkes Pulsar Timing Array Project

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

preprint2011arXiv

17 and 24 GHz observations of southern pulsars

We present observations of PSRs J0437-4715, J0738-4042, J0835-4510, J0908-4913, J1048-5832, J1622-4950, J1644-4559, J1721-3532 and J1740-3015 at 17 GHz using the Parkes radio telescope. All 9 were detected at 17 GHz, additionally, we detected PSR J0835-4510 and J1622-4950 at 24 GHz. Polarisation profiles of each pulsar and the variation with frequency are discussed. In general, we find that the highly polarised edge components of young pulsars continue to dominate their profiles at 17 GHz. Older pulsars (>10^5 years) appear to be almost completely depolarised. Our detection of PSR J0437-4715 is the highest frequency observation of a millisecond pulsar to date, and implies a luminosity at 17 GHz of 14 μJy kpc^2, and a mean spectral index of 2.2. We find that the spectral index of the magnetar PSR J1622-4950 is flat between 1.4 and 24 GHz, similar to the other known radio magnetars XTE J1810-197 and 1E 1547.0-5408. The profile is similar to that at 3.1 GHz, and is highly linearly polarised. Analysis of the frequency evolution of the profile of PSR J0835-4510 show that the profile is made of four components that vary with frequency only in their amplitude. The width and separation of the components remains fixed and the spectral index of each component can be determined independently.

preprint2011arXiv

On detection of the stochastic gravitational-wave background using the Parkes pulsar timing array

We search for the signature of an isotropic stochastic gravitational-wave background in pulsar timing observations using a frequency-domain correlation technique. These observations, which span roughly 12 yr, were obtained with the 64-m Parkes radio telescope augmented by public domain observations from the Arecibo Observatory. A wide range of signal processing issues unique to pulsar timing and not previously presented in the literature are discussed. These include the effects of quadratic removal, irregular sampling, and variable errors which exacerbate the spectral leakage inherent in estimating the steep red spectrum of the gravitational-wave background. These observations are found to be consistent with the null hypothesis, that no gravitational-wave background is present, with 76 percent confidence. We show that the detection statistic is dominated by the contributions of only a few pulsars because of the inhomogeneity of this data set. The issues of detecting the signature of a gravitational-wave background with future observations are discussed.

preprint2011arXiv

The High Time Resolution Universe Pulsar Survey IV: Discovery and polarimetry of millisecond pulsars

We present the discovery of six millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) survey for pulsars and fast transients carried out with the Parkes radio telescope. All six are in binary systems with approximately circular orbits and are likely to have white dwarf companions. PSR J1017-7156 has a high flux density and a narrow pulse width, making it ideal for precision timing experiments. PSRs J1446-4701 and J1125-5825 are coincident with gamma-ray sources, and folding the high-energy photons with the radio timing ephemeris shows evidence of pulsed gamma-ray emission. PSR J1502-6752 has a spin period of 26.7 ms, and its low period derivative implies that it is a recycled pulsar. The orbital parameters indicate it has a very low mass function, and therefore a companion mass much lower than usually expected for such a mildly recycled pulsar. In addition we present polarisation profiles for all 12 MSPs discovered in the HTRU survey to date. Similar to previous observations of MSPs, we find that many have large widths and a wide range of linear and circular polarisation fractions. Their polarisation profiles can be highly complex, and although the observed position angles often do not obey the rotating vector model, we present several examples of those that do. We speculate that the emission heights of MSPs are a substantial fraction of the light cylinder radius in order to explain broad emission profiles, which then naturally leads to a large number of cases where emission from both poles is observed.

preprint2011arXiv

The High Time Resolution Universe Survey - III. Single-pulse searches and preliminary analysis

We present the search methods and initial results for transient radio signals in the High Time Resolution Universe (HTRU) Survey. The HTRU survey's single-pulse search, the software designed to perform the search, and a determination of the HTRU survey's sensitivity to single pulses are described. Initial processing of a small fraction of the survey has produced 11 discoveries, all of which are sparsely-emitting neutron stars, as well as provided confirmation of two previously unconfirmed neutron stars. Most of the newly discovered objects lie in regions surveyed previously, indicating both the improved sensitivity of the HTRU survey observing system and the dynamic nature of the radio sky. The cycles of active and null states in nulling pulsars, rotating radio transients (RRATs), and long-term intermittent pulsars are explored in the context of determining the relationship between these populations, and of the sensitivity of a search to the various radio-intermittent neutron star populations. This analysis supports the case that many RRATs are in fact high-null-fraction pulsars (i. e. with null fraction >~0.95), and indicates that intermittent pulsars appear distinct from nulling pulsars in their activity cycle timescales. We find that in the measured population, there is a deficit of pulsars with typical emission timescales greater than ~300 s that is not readily explained by selection effects. The HTRU low-latitude survey will be capable of addressing whether this deficit is physical. We predict that the HTRU survey will explore pulsars with a broad range of nulling fractions (up to and beyond 0.999), and at its completion is likely to increase the currently known RRATs by a factor of more than two.

preprint2011arXiv

The High Time Resolution Universe Survey II: Discovery of 5 Millisecond Pulsars

We present the discovery of 5 millisecond pulsars found in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods from ~2.3 to ~7.5 ms, and all are in binary systems with orbital periods ranging from ~0.3 to ~150 d. In four of these systems, the most likely companion is a white dwarf, with minimum masses of ~0.2 Solar Masses. The other pulsar, J1731-1847, has a very low mass companion and exhibits eclipses, and is thus a member of the "black widow" class of pulsar binaries. These eclipses have been observed in bands centred near frequencies of 700, 1400 and 3000 MHz, from which measurements have been made of the electron density in the eclipse region. These measurements have been used to examine some possible eclipse mechanisms. The eclipse and other properties of this source are used to perform a comparison with the other known eclipsing and "black widow" pulsars. These new discoveries occupy a short-period and high-dispersion measure (DM) region of parameter space, which we demonstrate is a direct consequence of the high time and frequency resolution of the HTRU survey. The large implied distances to our new discoveries makes observation of their companions unlikely with both current optical telescopes and the Fermi Gamma-ray Space Telescope. The extremely circular orbits make any advance of periastron measurements highly unlikely. No relativistic Shapiro delays are obvious in any of the systems, although the low flux densities would make their detection difficult unless the orbits were fortuitously edge-on.

preprint2011arXiv

The Parkes Observatory Pulsar Data Archive

The Parkes pulsar data archive currently provides access to 144044 data files obtained from observations carried out at the Parkes observatory since the year 1991. Around 10^5 files are from surveys of the sky, the remainder are observations of 775 individual pulsars and their corresponding calibration signals. Survey observations are included from the Parkes 70cm and the Swinburne Intermediate Latitude surveys. Individual pulsar observations are included from young pulsar timing projects, the Parkes Pulsar Timing Array and from the PULSE@Parkes outreach program. The data files and access methods are compatible with Virtual Observatory protocols. This paper describes the data currently stored in the archive and presents ways in which these data can be searched and downloaded.

preprint2011arXiv

Transformation of a Star into a Planet in a Millisecond Pulsar Binary

Millisecond pulsars are thought to be neutron stars that have been spun-up by accretion of matter from a binary companion. Although most are in binary systems, some 30% are solitary, and their origin is therefore mysterious. PSR J1719-1438, a 5.7 ms pulsar, was detected in a recent survey with the Parkes 64m radio telescope. We show that it is in a binary system with an orbital period of 2.2 h. Its companion's mass is near that of Jupiter, but its minimum density of 23 g cm$^{-3}$ suggests that it may be an ultra-low mass carbon white dwarf. This system may thus have once been an Ultra Compact Low-Mass X-ray Binary, where the companion narrowly avoided complete destruction.

preprint2010arXiv

A 6.5-GHz Multibeam Pulsar Survey

A survey of the Galactic plane in the region $-60\degree \leq l \leq 30\degree$, $|b| \leq 0.25\degree$ was carried out using the seven-beam Parkes Methanol Multibeam (MMB) receiver, which operates at a frequency of 6.5 GHz. Three pulsars were discovered, and 16 previously known pulsars detected. In this paper we present two previously-unpublished discoveries, both with extremely high dispersion measures, one of which is very close, in angular distance, to the Galactic centre. The survey data also contain the first known detection, at radio frequencies, of the radio magnetar PSR J1550-5418. Our survey observation was made 46 days prior to that previously published and places constraints on the beginning of pulsed radio emission from the source. The detection of only three previously undiscovered pulsars argues that there are few pulsars in the direction of the inner Galaxy whose flux density spectrum is governed by a flat power law. However, these pulsars would be likely to remain undetected at lower frequencies due to the large amount of scatter broadening which affects pulsars with high values of dispersion measure. Surveys with future telescopes at high observing frequencies will, therefore, play an important role in the discovery of pulsars at the Galactic centre. By simulating pulsar surveys of the Galaxy with Phase 1 SKA at frequencies of 1.4 GHz and 10 GHz, we find that high-frequency observations are the only way to discover and observe the Galactic-centre pulsar population.

preprint2010arXiv

DSPSR: Digital Signal Processing Software for Pulsar Astronomy

DSPSR is a high-performance, open-source, object-oriented, digital signal processing software library and application suite for use in radio pulsar astronomy. Written primarily in C++, the library implements an extensive range of modular algorithms that can optionally exploit both multiple-core processors and general-purpose graphics processing units. After over a decade of research and development, DSPSR is now stable and in widespread use in the community. This paper presents a detailed description of its functionality, justification of major design decisions, analysis of phase-coherent dispersion removal algorithms, and demonstration of performance on some contemporary microprocessor architectures.

preprint2010arXiv

Measuring the mass of solar system planets using pulsar timing

High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

preprint2010arXiv

The Commensal Real-time ASKAP Fast Transients (CRAFT) survey

We are developing a purely commensal survey experiment for fast (<5s) transient radio sources. Short-timescale transients are associated with the most energetic and brightest single events in the Universe. Our objective is to cover the enormous volume of transients parameter space made available by ASKAP, with an unprecedented combination of sensitivity and field of view. Fast timescale transients open new vistas on the physics of high brightness temperature emission, extreme states of matter and the physics of strong gravitational fields. In addition, the detection of extragalactic objects affords us an entirely new and extremely sensitive probe on the huge reservoir of baryons present in the IGM. We outline here our approach to the considerable challenge involved in detecting fast transients, particularly the development of hardware fast enough to dedisperse and search the ASKAP data stream at or near real-time rates. Through CRAFT, ASKAP will provide the testbed of many of the key technologies and survey modes proposed for high time resolution science with the SKA.

preprint2010arXiv

The High Time Resolution Universe Pulsar Survey I: System configuration and initial discoveries

We have embarked on a survey for pulsars and fast transients using the 13-beam Multibeam receiver on the Parkes radio telescope. Installation of a digital backend allows us to record 400 MHz of bandwidth for each beam, split into 1024 channels and sampled every 64 us. Limits of the receiver package restrict us to a 340 MHz observing band centred at 1352 MHz. The factor of eight improvement in frequency resolution over previous multibeam surveys allows us to probe deeper into the Galactic plane for short duration signals such as the pulses from millisecond pulsars. We plan to survey the entire southern sky in 42641 pointings, split into low, mid and high Galactic latitude regions, with integration times of 4200, 540 and 270 s respectively. Simulations suggest that we will discover 400 pulsars, of which 75 will be millisecond pulsars. With ~30% of the mid-latitude survey complete, we have re-detected 223 previously known pulsars and discovered 27 pulsars, 5 of which are millisecond pulsars. The newly discovered millisecond pulsars tend to have larger dispersion measures than those discovered in previous surveys, as expected from the improved time and frequency resolution of our instrument.

preprint2010arXiv

The Sensitivity of the Parkes Pulsar Timing Array to Individual Sources of Gravitational Waves

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

preprint2009arXiv

Precision southern hemisphere VLBI pulsar astrometry II: Measurement of seven parallaxes

Accurate measurement of pulsar distances via astrometry using very long baseline interferometry enables the improvement of Galactic electron density distribution models, improving distance estimates for the vast majority of pulsars for which parallax measurements are unavailable. However, pulsars at southern declinations have been under-represented in previous interferometric astrometry campaigns. In order to redress this imbalance, we have conducted a two-year astrometric campaign targeting eight southern pulsars with the Australian Long Baseline Array. The program summarized in this paper has resulted in the measurement of seven new pulsar parallaxes, with success on objects down to a mean flux density of 0.8 mJy at 1600 MHz. Our results highlight the substantial uncertainties that remain when utilizing free electron density models for individual pulsar distances. Until this study, PSR J0630-2834 was believed to convert 16% of its spin-down energy into x-rays, but our measured parallax distance of 332 (+52 -40) pc has revised this value to <1%. In contrast, PSR J0108-1431 was found to be almost a factor of two more distant than previously thought, making its conversion of spin-down energy to x-rays the most efficient known (>1%). The 8.5 second radio pulsar J2144-3933 was found to be closer than previously predicted, making its apparent 1400 MHz radio luminosity the lowest of any known pulsar (20 microJy kpc^2). We have examined the growing population of neutron stars with accurate parallaxes to determine the effect of distance errors on the underlying neutron star velocity distribution, and find that typical distance errors may be biasing the estimated mean pulsar velocity upwards by 5%, and are likely to exaggerate the distribution's high-velocity tail.

preprint2009arXiv

Status Update of the Parkes Pulsar Timing Array

The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a successful detection of gravitational waves within a decade, or alternatively to rule out all current predictions for gravitational wave backgrounds formed by supermassive black-hole mergers.

preprint2009arXiv

The international pulsar timing array project: using pulsars as a gravitational wave detector

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

preprint2009arXiv

The Millisecond Radio Sky: Transients from a Blind Single Pulse Search

We present the results of a search for transient radio bursts of between 0.125 and 32 millisecond duration in two archival pulsar surveys of intermediate galactic latitudes with the Parkes multibeam receiver. Fourteen new neutron stars have been discovered, seven of which belong to the recently identified "rotating radio transients" (RRATs) class. Here we describe our search methodology, and discuss the new detections in terms of how the RRAT population relates to the general population of pulsars. The new detections indicate (1) that the galactic z-distribution of RRATs in the surveys closely resembles the distribution of pulsars, with objects up to 0.86 kpc from the galactic plane; (2) where measurable, the RRAT pulse widths are similar to that of individual pulses from pulsars of similar period, implying a similar beaming fraction; and (3) our new detections span a variety of nulling fractions, and thus we postulate that the RRATs may simply be nulling pulsars that are only "on" for less than a pulse period. Finally, the newly discovered object PSR J0941-39 may represent a link between pulsars and RRATs. This bizarre object was discovered as an RRAT, but in follow-up observations often appeared as a bright (~10 mJy) pulsar with a low nulling fraction. It is obvious therefore that a neutron star can oscillate between being an RRAT and a pulsar. Crucially, the sites of the RRAT pulses are coincident with the pulsar's emission, implying that the two emission mechanisms are linked, and that RRATs are not just pulsars observed from different orientations.

preprint2009arXiv

Timing stability of millisecond pulsars and prospects for gravitational-wave detection

Analysis of high-precision timing observations of an array of approx. 20 millisecond pulsars (a so-called "timing array") may ultimately result in the detection of a stochastic gravitational-wave background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a five-year timescale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within five years to a decade are good.

preprint2008arXiv

Gravitational wave detection using pulsars: status of the Parkes Pulsar Timing Array project

The first direct detection of gravitational waves may be made through observations of pulsars. The principal aim of pulsar timing array projects being carried out worldwide is to detect ultra-low frequency gravitational waves (f ~ 10^-9 to 10^-8 Hz). Such waves are expected to be caused by coalescing supermassive binary black holes in the cores of merged galaxies. It is also possible that a detectable signal could have been produced in the inflationary era or by cosmic strings. In this paper we review the current status of the Parkes Pulsar Timing Array project (the only such project in the Southern hemisphere) and compare the pulsar timing technique with other forms of gravitational-wave detection such as ground- and space-based interferometer systems.

preprint2008arXiv

Precision timing of PSR J0437-4715: an accurate pulsar distance, a high pulsar mass and a limit on the variation of Newton's gravitational constant

Analysis of ten years of high-precision timing data on the millisecond pulsar PSR J0437-4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, Pbdot, determined at the 1.5% level of precision (Dk = 157.0 +/- 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 Solar System ephemerides. The precise measurement of Pbdot allows a limit on the variation of Newton's gravitational constant, |Gdot/G| < 23 x 10^{-12} 1/yr. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a(Sun)/c| < 1.5 x 10^{-18} 1/s at 95% confidence, and derive a pulsar mass, m(psr) = 1.76 +/- 0.20 M, one of the highest estimates so far obtained.

preprint2004arXiv

PSR J1022+1001: Profile Stability & Precision Timing

We present an investigation of the morphology and arrival times of integrated radio pulses from the binary millisecond pulsar PSR J1022+1001. This pulsar is renowned for its poor timing properties, which have been postulated to originate from variability in its average pulse profile. Although a sub-class of long-period pulsars are known to exhibit mode changes that give rise to very large deviations in their integrated profiles, this was the first millisecond pulsar thought to have an unstable mean profile. As part of a precision timing program at the Parkes radio telescope we observed this pulsar between January 2003 and March 2004 using a coherent de-dispersion system (CPSR2). A study of morphological variability during our brightest observations suggests that the pulse profile varies by at most a few percent, similar to the uncertainty in our calibration. Unlike previous authors, we find that this pulsar times extremely well. In five minute integrations of 64 MHz bands we obtain a weighted RMS residual of just 2.27 microseconds. The reduced chi-squared of our best fit is 1.43, which suggests that this pulsar can be timed to high accuracy with standard cross-correlation techniques. Combining relativistic constraints with the pulsar mass function and consideration of the Chandrasekhar mass limit on the white dwarf companion, we can constrain the inclination angle of the system to lie within the range 37 < i < 56 degrees. For reasonable pulsar masses, this suggests that the white dwarf is at least 0.9 solar masses. We also find evidence for secular evolution of the projected semi-major axis.

preprint2001arXiv

The Swinburne Intermediate Latitude Pulsar Survey

We have conducted a survey of intermediate Galactic latitudes using the 13-beam 21-cm multibeam receiver of the Parkes 64-m radio telescope. The survey covered the region enclosed by 5 deg < |b| < 15 deg and -100 deg < l < 50 deg with 4,702 processed pointings of 265 s each, for a total of 14.5 days of integration time. Thirteen 2x96-channel filterbanks provided 288 MHz of bandwidth at a centre frequency of 1374 MHz, one-bit sampled every 125 microsec and incurring ~DM/13.4 cm^-3 pc samples of dispersion smearing. The system was sensitive to slow and most millisecond pulsars in the region with flux densities greater than approximately 0.3--1.1 mJy. Offline analysis on the 64-node Swinburne workstation cluster resulted in the detection of 170 pulsars of which 69 were new discoveries. Eight of the new pulsars, by virtue of their small spin periods and period derivatives, may be recycled and have been reported elsewhere. The slow pulsars discovered are typical of those already known in the volume searched, being of intermediate to old age. Several pulsars experience pulse nulling and two display very regular drifting sub-pulses. We discuss the new discoveries and provide timing parameters for the 48 slow pulsars for which we have a phase-connnected solution.

preprint1996arXiv

Timing Measurements and Their Implications for Four Binary Millisecond Pulsars

We present timing observations of four millisecond pulsars, using data obtained over three years at the ATNF Parkes and NRAL Jodrell Bank radio telescopes. Astrometric, spin, and binary parameters are updated, and substantially improved for three pulsars, PSRs J0613-0200, J1045-4509 and J1643-1224. We have measured the time variation of the projected semi-major axis of the PSR J0437-4715 orbit due to its proper motion, and use it to constrain the inclination of the orbit and the mass of the companion. Some evidence is found for changes in the dispersion measures of PSRs J1045-4509 and J1643-1224. Limits are placed on the existence of planetary mass companions, ruling out companions with masses and orbits similar to the terrestrial planets of the solar system for eight pulsars.