No Arabic abstract
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 u/ u 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.
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$.
We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613$-$0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank TOAs recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of ${-1.6(39) times 10^{-19},text{Hz} text{s}^{-1}}$ are measured at the epoch of MJD 50888(30). After PSR B1821$-$24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of ${Delta u/ u=2.5(1) times 10^{-12}}$, which is several times smaller than the previous smallest glitch. PSR J0613$-$0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of $sim 50$% that another glitch will be observed in a timing array pulsar within 10 years.
We have observed a large glitch in the Crab pulsar (PSR B0531+21). The glitch occurred around MJD 58064 (2017 November 8) when the pulsar underwent an increase in the rotation rate of $Delta u = 1.530 times 10^{-5}$ Hz, corresponding to a fractional increase of $Delta u / u = 0.516 times 10^{-6}$ making this event the largest glitch ever observed in this source. Due to our high-cadence and long-dwell time observations of the Crab pulsar we are able to partially resolve a fraction of the total spin-up of the star. This delayed spin-up occurred over a timescale of $sim$1.7 days and is similar to the behaviour seen in the 1989 and 1996 large Crab pulsar glitches. The spin-down rate also increased at the glitch epoch by $Delta dot{ u} / dot{ u} = 7 times 10^{-3}$. In addition to being the largest such event observed in the Crab, the glitch occurred after the longest period of glitch inactivity since at least 1984 and we discuss a possible relationship between glitch size and waiting time. No changes to the shape of the pulse profile were observed near the glitch epoch at 610 MHz or 1520 MHz, nor did we identify any changes in the X-ray flux from the pulsar. The long-term recovery from the glitch continues to progress as $dot{ u}$ slowly rises towards pre-glitch values. In line with other large Crab glitches, we expect there to be a persistent change to $dot{ u}$. We continue to monitor the long-term recovery with frequent, high quality observations.
Glitches are sudden jumps in the spin frequency of pulsars believed to originate in the superfluid interior of neutron stars. Superfluid flow in a model neutron star is simulated by solving the equations of motion of a two-component superfluid consisting of a viscous proton-electron plasma and an inviscid neutron condensate in a spherical Couette geometry. We examine the response of our model neutron star to glitches induced in three different ways: by instantaneous changes of the spin frequency of the inner and outer boundaries, and by instantaneous recoupling of the fluid components in the bulk. All simulations are performed with strong and weak mutual friction. It is found that the maximum size of a glitch that originates in the bulk decreases as the mutual friction strengthens. It is also found that mutual friction determines the fraction of the frequency jump which is later recovered, a quantity known as the healing parameter. These behaviours may explain some of the diversity in observed glitch recoveries.
Since its discovery as a pulsar in 2000, the central compact object (CCO) 1E 1207.4-5209 in the supernova remnant PKS 1209-51/52 had been a stable 0.424 s rotator with an extremely small spin-down rate and weak (Bs ~ 9E10 G) surface dipole magnetic field. In 2016 we observed a glitch from 1E 1207.4-5209 of at least Delta f/f = (2.8+/-0.4)E-9, which is typical in size for the general pulsar population. However, glitch activity is closely correlated with spin-down rate fdot, and pulsars with fdot as small as that of 1E 1207.4-5209 are never seen to glitch. Unlike in glitches of ordinary pulsars, there may have been a large increase in fdot as well. The thermal X-ray spectrum of 1E 1207.4-5209, with its unique cyclotron absorption lines that measure the surface magnetic field strength, did not show any measurable change after the glitch, which rules out a major disruption in the dipole field as a cause or result of the glitch. A leading theory of the origin and evolution of CCOs, involving prompt burial of the magnetic field by fall-back of supernova ejecta, might hold the explanation for the glitch.