Pulsar data analysis with PSRCHIVE

PSRCHIVE is an open-source, object-oriented, scientific data analysis software library and application suite for pulsar astronomy. It implements an extensive range of general-purpose algorithms for use in data calibration and integration, statistical analysis and modeling, and visualisation. These are utilised by a variety of applications specialised for tasks such as pulsar timing, polarimetry, radio frequency interference mitigation, and pulse variability studies. This paper presents a general overview of PSRCHIVE functionality with some focus on the integrated interfaces developed for the core applications.Comment: 21 pages, 5 figures; tutorial presented at IPTA 2010 meeting in Leiden merged with talk presented at 2011 pulsar conference in Beijing; includes further research and development on algorithms for RFI mitigation and TOA bias correctio

Long Term Variability of a Black Widow's Eclipses -- A Decade of PSR J2051−-0827

In this paper we report on $\sim10$ years of observations of PSR J2051$-$0827, at radio frequencies in the range 110--4032 MHz. We investigate the eclipse phenomena of this black widow pulsar using model fits of increased dispersion and scattering of the pulsed radio emission as it traverses the eclipse medium. These model fits reveal variability in dispersion features on timescales as short as the orbital period, and previously unknown trends on timescales of months--years. No clear patterns are found between the low-frequency eclipse widths, orbital period variations and trends in the intra-binary material density. Using polarisation calibrated observations we present the first available limits on the strength of magnetic fields within the eclipse region of this system; the average line of sight field is constrained to be $10^{-4}$ G $\lesssim B_{||} \lesssim 10^2$ G, while for the case of a field directed near-perpendicular to the line of sight we find $B_{\perp} \lesssim 0.3$ G. Depolarisation of the linearly polarised pulses during the eclipse is detected and attributed to rapid rotation measure fluctuations of $\sigma_{\text{RM}} \gtrsim 100$ rad m$^{-2}$ along, or across, the line of sights averaged over during a sub-integration. The results are considered in the context of eclipse mechanisms, and we find scattering and/or cyclotron absorption provide the most promising explanation, while dispersion smearing is conclusively ruled out. Finally, we estimate the mass loss rate from the companion to be $\dot{M}_{\text{C}} \sim 10^{-12} M_\odot$ yr$^{-1}$, suggesting that the companion will not be fully evaporated on any reasonable timescale

Detection of a glitch in the pulsar J1709-4429

We report the detection of a glitch event in the pulsar J1709$-$4429 (also known as B1706$-$44) during regular monitoring observations with the Molonglo Observatory Synthesis Telescope (UTMOST). The glitch was found during timing operations, in which we regularly observe over 400 pulsars with up to daily cadence, while commensally searching for Rotating Radio Transients, pulsars, and FRBs. With a fractional size of $\Delta\nu/\nu \approx 52.4 \times10^{-9}$, the glitch reported here is by far the smallest known for this pulsar, attesting to the efficacy of glitch searches with high cadence using UTMOST.Comment: 3 pages, 1 figur

The UTMOST pulsar timing programme II:Timing noise across the pulsar population

While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesis Telescope (MOST). We reanalysed the spin and spin-down changes associated with nine previously reported pulsar glitches, report the discovery of three new glitches and four unusual glitch-like events in the rotational evolution of PSR J1825$-$0935. We develop a refined Bayesian framework for determining how red noise strength scales with pulsar spin frequency ($\nu$) and spin-down frequency ($\dot{\nu}$), which we apply to a sample of 280 non-recycled pulsars. With this new method and a simple power-law scaling relation, we show that red noise strength scales across the non-recycled pulsar population as $\nu^{a} |\dot{\nu}|^{b}$, where $a = -0.84^{+0.47}_{-0.49}$ and $b = 0.97^{+0.16}_{-0.19}$. This method can be easily adapted to utilise more complex, astrophysically motivated red noise models. Lastly, we highlight our timing of the double neutron star PSR J0737$-$3039, and the rediscovery of a bright radio pulsar originally found during the first Molonglo pulsar surveys with an incorrectly catalogued position.Comment: Accepted by MNRAS. 28 pages, 8 figures, 8 table

The gravitational-wave background null hypothesis: Characterizing noise in millisecond pulsar arrival times with the Parkes Pulsar Timing Array

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $\gamma=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission.Comment: 18 pages, 10 figures. Accepted for publication in ApJ