Type: Article
Publication Date: 2013-10-31
Citations: 4
DOI: https://doi.org/10.1088/0004-637x/778/1/31
We study the pulsar timing properties and the data analysis methods of glitch recoveries. In some cases one first fits the times of arrival (TOAs) to obtain the "time-averaged" frequency ν and its first derivative , and then fits models to them. However, our simulations show that ν and obtained in this way are systematically biased unless the time intervals between the nearby data points of TOAs are smaller than about 104 s, which is much shorter than typical observation intervals. Alternatively, glitch parameters can be obtained by fitting the phases directly with relatively smaller biases, but the initial recovery time scale is usually chosen by eye, which may introduce a strong bias. We also construct a phenomenological model by assuming a pulsar spin-down law of with G(t) = 1 + κe−t/τ for a glitch recovery, where H0 is a constant and κ and τ are the glitch parameters to be found. This model can reproduce the observed data of slow glitches from B1822−09 and the giant classical glitch of B2334+61, with κ < 0 or κ > 0, respectively. We then use this model to simulate TOA data and test several fitting procedures for a glitch recovery. The best procedure is (1) to use a very high order polynomial (e.g., to 50th order) to precisely describe the phase, (2) obtain ν(t) and from the polynomial, then (3) obtain the glitch parameters from ν(t) or . Finally, the uncertainty in the starting time t0 of a classical glitch causes uncertainties in some glitch parameters, but less so for a slow glitch, t0, of which can be determined from data.