<HTML>
<TITLE>African Skies 4 - Radio Pulsar Glitch Studies</TITLE>
<BODY BGCOLOR="WHITE">
<center><h1>Radio Pulsar Glitch Studies</h1>
<h3>J.O. Urama<SUP>1</SUP></h3>
<EM>Hartebeesthoek Radio Astronomy Observatory, South Africa<br>
<a href="mailto"johnson@bootes.hartrao.ac.za">
johnson@bootes.hartrao.ac.za</a></EM>
<P>
<h3>P.N. Okeke</h3>
<EM>Space Research Centre and Department of Physics
& Astronomy, University of Nigeria, Nsukka</center></em>
<P>
<small><sup>1</sup> On leave from the Department of Physics &
Astronomy, University of Nigeria, Nsukka</small>
<p>
<B>Abstract</B>. Timing observations of glitches in radio
pulsars is currently one of the best probes of the neutron star
interior. Such observations have led to a number of models describing
the various components of the neutron star and their mutual coupling.
Here, we briefly summarize some of our current work in this regard.
<P>
<B>Sommaire</B>. Les observations des changements soudains de fr&#233;quence
des pulsars radio sont de nos jours un des meilleurs sondages de
l'int&#233;rieur d'une &#233;toile &#224; neutrons. De telles observations ont
conduit &#224; plusieurs mod&#232;les qui d&#233;crivent les diff&#233;rentes
composantes de l'&#233;toile &#224; neutrons et leur couplage mutuel. Nous
r&#233;sumons ici bri&#232;vement nos recherches actuelles &#224; cet &#233;gard.
<P>
<H2><A NAME="SECTION00010000000000000000">
1 Introduction</A>
</H2>
<P>
Radio pulsars were discovered as sequences of rapid (~ 1Hz) and
regular radio pulses in 1967, by Jocelyn Bell and Anthony Hewish (Hewish
<I>et al.</I>, 1968). Currently about 1000 pulsars are known and the galactic
population is estimated at around 10<SUP>5</SUP>.
<P>
Shorly after their discovery, the pulses were firmly associated with
highly magnetized, rapidly rotating, neutron stars. Neutron stars are
one of the final evolutionary stages of stars formed by the collapse of
a massive star in a supernova explosion. Containing typically 1.4 solar
masses within a diameter of about 15 km, they are also the densest known
form of directly-observable matter in the Universe.
<P>
Although neutron stars have now been recognized in a variety of stellar
systems such as X-ray binaries, radio pulsars are by far the most common
observable manifestation of the fascinating neutron star. The emitted
radio pulses allow us to directly measure the rotation of the underlying star.
This rotation rate is, in comparison with most other astronomical measurements,
exceptionally stable and easy to measure with high accuracy. Glitches are
perturbations to this regular behaviour which give us a window into the
interior of the neutron star, as described below.
<P>
<H2><A NAME="SECTION00020000000000000000">
2 Pulsar Glitches and the Neutron Star Interior</A>
</H2>
<b>2.1 Glitches</b>
<P>
Pulsars slow down steadily, due to the loss of kinetic (rotational)
energy via the emission of high-energy particles and radiation. In other
words, the pulsar spin rate, <IMG
 WIDTH="20" HEIGHT="29" ALIGN="MIDDLE" BORDER="0"
 SRC="img16.gif"
 ALT="$\Omega,$"> has a negative derivative
(<!-- MATH: $\dot{\Omega}$ -->
<IMG  WIDTH="16" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
 SRC="img17.gif"
 ALT="$\dot{\Omega}$">). 
<!-- MATH: $\dot{\Omega}$ -->
<IMG  WIDTH="16" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
 SRC="img17.gif"
 ALT="$\dot{\Omega}$"> is easily observable and is
traditionally used to determine the age of the pulsar.
<P>
Outside the steady decrease in <IMG
 WIDTH="20" HEIGHT="29" ALIGN="MIDDLE" BORDER="0"
 SRC="img16.gif"
 ALT="$\Omega,$"> all pulsars exhibit to some
extent rotational irregularities known as timing noise, which is
observable as random wandering in either the pulse phase or its
frequency. Sudden changes in the rotational frequency are known as
glitches, and may or may not be related to timing noise. Glitches
involve relative increases in the pulsar rotational frequency of up to 1
part in 10<SUP>6</SUP>. Since the rotation rate can be measured to an accuracy
of up to 10<SUP>-12</SUP>, such events are easy to recognize.
<P>
<TABLE>
<small><center>Table 1: Period (P), period-derivative (<IMG
 WIDTH="17" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
 SRC="img18.gif"
 ALT="$\dot{P}$">), and the characteristic age (<IMG
 WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img19.gif"
 ALT="$\tau$">)
of PSRs 1046-58 and 1737-30 as compared to those of the Vela pulsar. 
The parameters were obtained from the Tables of Taylor 
<I> et al.</I> (1993).</small></center>
<p>
<center><table cellpadding=2 border=1>
<TR>
<TD>
<TR><TD ALIGN="CENTER">PSR</TD>
<TD ALIGN="CENTER">P (s)</TD>
<TD ALIGN="CENTER"><IMG
 WIDTH="17" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
 SRC="img18.gif"
 ALT="$\dot{P}$"></TD>
<TD ALIGN="CENTER">Log <IMG
 WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img19.gif"
 ALT="$\tau$">(yr)</TD>
</TR>
<TR><TD ALIGN="CENTER">B1046-58</TD>
<TD ALIGN="CENTER">0.1236526</TD>
<TD ALIGN="CENTER">96.93</TD>
<TD ALIGN="CENTER">4.31</TD>
</TR>
<TR><TD ALIGN="CENTER">B1737-30</TD>
<TD ALIGN="CENTER">0.6066433</TD>
<TD ALIGN="CENTER">465.67</TD>
<TD ALIGN="CENTER">4.31</TD>
</TR>
<TR><TD ALIGN="CENTER">Vela</TD>
<TD ALIGN="CENTER">0.0892991</TD>
<TD ALIGN="CENTER">126.68</TD>
<TD ALIGN="CENTER">4.05</TD></TR>
</TABLE></center>
<p>
While the slow-down is believed to be caused by processes above the
surface of the neutron star, pulsars are (unusual amongst stars) almost
all isolated bodies, and the glitch is generally attributed to internal
processes. The few intensive observations of post-glitch behaviour
reveal recovery over a range of time scales. These observations have
been used to study the pulsar interior, in particular the various
components of the interior and the manner in which they couple with each
other.
<P>
To date about eighty glitches have been noted, involving over thirty pulsars.
The extremely young PSR 0531+21 (Crab) and the adolescent
PSRs 0833-45 (Vela) and 1737-30 are responsible for almost half of
these glitches. Evidence is growing in support of the view that glitch
characteristics such as frequency of occurrence and rate of post-glitch
recovery are dependent on the age of the pulsar. Unfortunately, only in
the Crab and Vela pulsars have multiple glitches been monitored
intensively. Similarly, intense observations from a larger sample of
pulsars is required to distinguish between competing theories of the
neutron star interior.
<P>
<center><img src="fig1.gif">
<p>
<small><b>Figure 1:</b>
Integrated pulse profiles showing: (a) 72,750
integrations of 1046-58 at 18 cm; and <br>
(b) 18,000 integrations of 1737-30 at 13 cm.</small></center>
<P>
<b><A NAME="SECTION00022000000000000000">
2.2 Glitch Monitoring Project</b>
<P>
Pulsar observations using the Hartebeesthoek Radio Astronomy Observatory
(HartRAO) 26-m radio telescope commenced in 1984. Observations are made
at 1.7 and 2.3 GHz. In 1988, software was installed to analyse observations
in near-real time and then initiate continuous observations should a glitch
be detected. The Observatory also has software that allows short
observations of a range of astronomical objects to be scheduled
automatically at selected times. These developments allow for a rapid
and automatic response to glitches, and the Vela pulsar became the first
target of a new ``glitch detection'' project. The project has been quite
successful in providing good observational coverage of immediate postglitch
behaviour of five glitches in the Vela pulsar (Flanagan, 1995). Recently, two
more young pulsars (1046-58 and 1737-30) have been included in the list of
frequently-observed pulsars, although the low signal-to-noise ratio of these
observations prevented on-line glitch-detection with automatically-scheduled 
follow-up observations.
<P>
<center><img src="urama.jpg">
<p>
<small><b>Figure 2:</b>
The rotational parameters of PRS <I>B</I>1046-58 showing the glitch
that occured on MJD 50791.5 &#177 0.5. <br>
(a) shows the residual frequency after subtracting the pre-glitch model 
while (b) is the frequency derivative.<br> 
In panel (a) the uncertainties are much smaller than the
size of the symbols.</small></center>
<p>
PSRs 1737-30 and 1046-58 are of similar age to Vela (Table 1). The
former is the most frequently-glitching pulsar known, although its
post-glitch behaviour has not yet been observed as well as that of Vela.
No glitch observations of the more southern PSR 1046-58 appear to have
been published; because of its age and the similarity of its spin-down
rate to that of Vela, we expect it also to undergo frequent glitches.
PSR 1046-58 is observed twice daily at approximately 8-hourly intervals
while PSR 1737-30 is observed three times daily at approximately 5-hourly
intervals. A single integration consists of 5000 samples of PSR
1737-30 and 1089 samples of PSR 1046-58; the integration lasting 10-20
minutes depending on the frequency of observation. Fig. 1 shows
samples of integrated pulse profiles for these pulsars.
<P>
Manual checks of phase drift are done after most of the observations
for any evidence of a glitch. The first glitch in PSR 1046-58 was observed
in this monitoring programme. It occurred on December 9, 1997, about 200 days
after the commencement of intensive observations on this pulsar. The spin-up
was rather large in magnitude and was accompanied by an increase in spin-down
rate.  We also observed a large glitch on PSR 1737-30 (the largest so far
on this pulsar) on May 3, 1998. Figs. 2 and 3 show the rotational behaviour
of the two pulsars over a 400-day period.
<P>
<b><A NAME="SECTION00023000000000000000">
2.3 Glitch Recovery Model</b>
<P>
Early observations of pulsar response to a glitch on time scales of
weeks provided strong evidence for a superfluid interior(<I>e.g.</I>
Baym <I>et al.</I>, 1969) - all forms of ``normal'' matter would respond to
such a disturbance on unobservably short time-scales. A neutron star, as
shown in Fig. 4, is believed to have a rigid outer crust which consists
primarily of iron nuclei; an inner crust of neutron-rich nuclei,
electrons and <SUP>1</SUP><I>S</I><SUB>0</SUB> neutron superfluids; an outer 
(fluid) core containing <SUP>3</SUP><I>P</I><SUB>2</SUB> superfluid neutrons 
with an admixture of <SUP>1</SUP><I>S</I><SUB>0</SUB> super-conducting protons 
and normal electrons; and an inner (plasma) core that may possibly contain 
pion condensate, condensed kaons or matter in some exotic state (Alpar, 1989).
<P>
<center><img src="urama1.jpg">
<p>
<small><b>Figure 3:</b>
The rotational parameters of PRS <I>B</I>1737-58 showing the glitch
that occured on MJD 50936.8 &#177 0.1.<br>
(a) shows the residual frequency after subtracting the pre-glitch model 
while (b) is the frequency derivative. <br>
In panel (a) the uncertainties are much smaller than the
size of the symbols.</small></center>
<p>
Pulsar glitch recovery is now known to occur over all observed time scales,
from a few hours to years. These recovery time scales have been
explained in terms of the superfluid components of the neutron star
interior and their coupling to the lattice crust. As the quality and
number of observations has increased, so too has the number of proposed
possible mechanisms of inter-component coupling. Some of the currently
surviving models are, briefly: the vortex creep model (Alpar <I>et al.</I>,
1984, 1993, 1996); the core ``shell'' model (Sedrakian <I>et al.</I>, 1995); the
three-component model (Takatsuka &amp; Tamagaki, 1989) and the vortex
corotating model (Jones, 1990). In both the vortex creep and the vortex
corotating models, glitch behaviours are attributed to the interaction
between the crustal superfluid and the rest of the star while the core
``shell'' model is based on the core superfluid coupling/decoupling from
the rest of the star. In the three-component model both the core and
crust superfluids contribute to glitch behaviours. These and other
models are yet to provide adequate explanations for all the observed
glitch characteristics - glitch sizes, recovery, interval, etc.
<P>
In a modified three-component model we extended the original
three-component model to accomodate all the observed glitch recovery
time scales (Urama & Okeke, 1999). The model is based on the principle
of the corotation of both the core and crust superfluids which are not
pinned. Here, the shortest recovery timescale has been attributed to the
core superfluid while the interglitch period is seen to result from the
crust superfluid.
<P>
<center><img src="fig5.gif">
<p>
<small><b>Figure 4:</b>
A schematic representation of a possible cross-section of a
<!-- MATH: $1.4M_{\odot}$ -->
<IMG WIDTH="51" HEIGHT="29" ALIGN="MIDDLE" BORDER="0"
 SRC="img25.gif"
 ALT="$1.4M_{\odot}$"> neutron star.</small></center>
<P>
<H2><A NAME="SECTION00030000000000000000">
3 Discussion and Conclusion</A>
</H2>
<P>
The Vela and Crab pulsars are the only pulsars whose glitch recoveries
have been monitored in sufficient detail. While the Vela pulsar has
generally shown giant glitches, all the Crab pulsar glitches have been
small. For the Crab pulsar, it is shown that the spin-down rate,
<!-- MATH: $\dot{\Omega}_{c},$ -->
<IMG  WIDTH="27" HEIGHT="37" ALIGN="MIDDLE" BORDER="0"
 SRC="img27.gif"
 ALT="$\dot{\Omega}_{c},$"> never returns to its preglitch value; there is a
permanent offset in </EM>
<!-- MATH: $\dot{\Omega}_{c}.$ -->
<IMG
 WIDTH="27" HEIGHT="37" ALIGN="MIDDLE" BORDER="0"
 SRC="img28.gif"
 ALT="$\dot{\Omega}_{c}.$"> The issue of whether the same
mechanisms are involved in Vela-type and Crab-type glitches continues to
be contentious. The differences in their glitch characteristics have
been too easily attributed to the age difference between them.
Unfortunately all the data on real-time glitch detection have been
obtained only on these two. The HartRAO pulsar project has contributed to
the debate on component interaction with the first observation of short
(~ hours) recovery following a glitch in the Vela pulsar (Flanagan,
1990), and of a pair of glitches with atypical characteristics. We hope
that by extending this project to include other young pulsars we can
also shed some light on the structural evolution of the pulsar interior.
<P>
<B>Acknowledgements</B>
<P>
One of us, J.O.U., acknowledges the support of HartRAO and an IAU Commission
38 travel grant towards his visit to South Africa. He also wishes to
thank the staff of HartRAO and Claire Flanagan for hospitality and
wonderful assistance with this work.
<P>
<b>References</b>
<OL>
<LI>Alpar, M.A., Chau, H.F., Cheng, K.S., and Pines, D. 1996, <EM>Ap.J.,</EM> 
459, 706
<LI>Alpar, M.A., Chau, H.F., Cheng, K.S., and Pines, D. 1993, <EM>Ap.J.,</EM> 
409, 345
<LI>Alpar, M.A.  1989, <I>Inside Neutron Stars</I>: H. &#214;gelman and
E.&nbsp;P.&nbsp;J. van den Heuvel (eds.), Timing Neutron Stars (Kluwer Academic
Publishers), 431
<LI>Alpar, M.A., Langer, S., and Sauls, J.A. 1984, <EM>Ap.J.,</EM> 282, 533
<LI>Baym, G., Pethick, C.J., and Pines, D. 1969, <EM>Nature,</EM> 224, 673
<LI>Flanagan, C.S. 1990, <EM>Nature,</EM> 345,416
<LI>Flanagan, C.S. 1995, PhD Thesis, Rhodes University
<LI>Hewish, A., Bell, S.J., Pilkington, J.D.H., Scott, P.F., and Collins,
R.A.  1968, <EM>Nature,</EM> 217, 709
<LI>Jones, P.B.  1990, <EM>M.N.R.A.S.,</EM> 246, 315
<LI>Sedrakian, A.D., Sedrakian, D.M., Cordes, J.M., and Terzian Yervant
1995, <EM>Ap.J.,</EM> 447, 324
<LI>Takatsuka, T., and Tamagaki, R. 1989, <EM>Prog. of Theor. Phys.,</EM> 82, 945
<LI>Taylor, J.H., Manchester, R.N., and Lyne, A.G. 1993, <EM>Ap.J. Suppl.
Series,</EM> 88, 529
<LI>Urama, J.O., and Okeke, P.N. 1999, <EM>Astrophys. Lett. &amp; Comm.,</EM> submitted
</ol>
<p><a href="as4.html"><img src="../backarr.gif" border=0></a>
<P>
<ADDRESS>
<I>WGSSA</I>
<BR><I>2000-03-07</I>
</ADDRESS>
</BODY>
</HTML>
