Pulsar glitch detection with a hidden Markov model

Abstract

Pulsar timing experiments typically generate a phase-connected timing solution from a sequence of times-of-arrival (TOAs) by absolute pulse numbering, i.e. by fitting an integer number of pulses between TOAs in order to minimize the residuals with respect to a parametrized phase model. In this observing mode, rotational glitches are discovered, when the residuals of the no-glitch phase model diverge after some epoch, and glitch parameters are refined by Bayesian follow-up. Here an alternative, complementary approach is presented which tracks the pulse frequency $f$ and its time derivative $df/dt$ with a hidden Markov model (HMM), whose dynamics include stochastic spin wandering (timing noise) and impulsive jumps in $f$ and $df/dt$ (glitches). The HMM tracks spin wandering explicitly, as a specific realization of a discrete-time Markov chain. It discovers glitches by comparing the Bayes factor for glitch and no-glitch models. It ingests standard TOAs for convenience and, being fully automated, allows performance bounds to be calculated quickly via Monte Carlo simulations. Practical, user-oriented plots are presented of the false alarm probability and detection threshold (e.g. minimum resolvable glitch size) versus observational scheduling parameters (e.g. TOA uncertainty, mean delay between TOAs) and glitch parameters (e.g. transient and permanent jump sizes, exponential recovery time-scale). The HMM is also applied to $\sim 1$ yr of real data bracketing the 2016 December 12 glitch in PSR J0835-4510 as a proof of principle. It detects the known glitch and confirms that no other glitch exists in the same data with size $>10^{-7}f$.

Accepted for publication in The Astrophysical Journal

References in corpus (19)

• Neutron Star Observations: Prognosis for Equation of State Constraints
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• Models of Pulsar Glitches
• Detection of 107 glitches in 36 southern pulsars
• Coherent Emission Mechanisms in Astrophysical Plasmas
• Avalanche dynamics of radio pulsar glitches
• Predicting the Starquakes in PSR J0537-6910
• A hydrodynamical trigger mechanism for pulsar glitches
• 30 Glitches in slow pulsars
• The UTMOST pulsar timing programme II: Timing noise across the pulsar population
• Directed searches for continuous gravitational waves from binary systems: parameter-space metrics and optimal Scorpius X-1 sensitivity
• Hidden Markov model tracking of continuous gravitational waves from a binary neutron star with wandering spin. II. Binary orbital phase tracking
• SOAP: A generalised application of the Viterbi algorithm to searches for continuous gravitational-wave signals
• Wide-band Profile Domain Pulsar Timing Analysis
• Empirically extending the range of validity of parameter-space metrics for all-sky searches for gravitational-wave pulsars
• Generative pulsar timing analysis
• Generating quasiperiodic pulsar glitches using a state-dependent Poisson process
• Diagnostics of timing noise in middle aged pulsars
• One year of monitoring the Vela pulsar using a Phased Array Feed