Hydrodynamic instabilities and turbulence in neutron stars have been suggested to be related to observable spin variations in pulsars, such as spin glitches, timing noise, and precession (nutation). Accounting for the stabilizing effects of the stellar magnetic field, we revisit the issue of whether the inertial modes of a neutron star can become unstable when the neutron and proton condensates flow with respect to one another. The neutron and proton condensates are coupled through the motion of imperfectly pinned vorticity (vortex slippage) and vortex-mediated scattering (mutual friction). Two-stream instabilities that occur when the two condensates rotate with respect to one another in the outer core are stabilized by the toroidal component of the magnetic field. This stabilization occurs when the Alfv\'en speed of the toroidal component of the magnetic field becomes larger than the relative rotational velocity of the condensates, corresponding to toroidal field strengths in excess of $\simeq 10^{10}\,{\rm G}$. In contrast with previous studies, we find that spin down of a neutron star under a steady torque is stable. The Donnelly--Glaberson instability is not stabilized by the magnetic field, and could play an important role if neutron stars undergo precession.

Motivated by stability considerations and observational evidence, we argue that magnetars possess highly-tangled internal magnetic fields. We propose that the quasi-periodic oscillations (QPOs) seen to accompany giant flares can be explained as torsional modes supported by a tangled magnetic field, and we present a simple model that supports this hypothesis for SGR 1900+14. Taking the strength of the tangle as a free parameter, we find that the magnetic energy in the tangle must dominate that in the dipolar component by a factor of $\sim 14$ to accommodate the observed 28 Hz QPO. Our simple model provides useful scaling relations for how the QPO spectrum depends on the bulk properties of the neutron star and the tangle strength. The energy density in the tangled field inferred for SGR 1900+14 renders the crust nearly dynamically irrelevant, a significant simplification for study of the QPO problem. The predicted spectrum is about three times denser than observed, which could be explained by preferential mode excitation or beamed emission. We emphasize that field tangling is needed to stabilize the magnetic field, so should not be ignored in treatment of the QPO problem.

We propose a scenario for the quasi-periodic oscillations observed in magnetar flares wherein a tangled component of the stellar magnetic field introduces nearly isotropic stress that gives the fluid core of the star an effective shear modulus. In a simple, illustrative model of constant density, the tangled field eliminates the problematic Alfv\'en continuum that would exist in the stellar core for an organized field. For a tangled field energy density comparable to that inferred from the measured dipole fields of $\sim 10^{15}$ G in SGRs 1806-20 and 1900+14, torsional modes exist with fundamental frequencies of about 20 Hz, and mode spacings of $\sim 10$ Hz. For fixed stellar mass and radius, the model has only one free parameter, and can account for {\em every} observed QPO under 160 Hz to within 3 Hz for both SGRs 1806-20 and 1900+14. The combined effects of stratification and crust stresses generally decrease the frequencies of torsional oscillations by $<10$% for overtones and increase the lowest-frequency fundamentals by up to 50%, and so the star can be treated as having constant density to a generally good first approximation. We address the issue of mode excitation by sudden readjustment of the stellar magnetosphere. While the total energy in excited modes is well within the energy budget of giant flares, the surface amplitude is $< 10^{-3}$ of the stellar radius for global oscillations, and decreases strongly with mode frequency. The 626 Hz QPO reported for SGR 1806-20 is particularly problematic to excite beyond a surface amplitude of $10^{-6}$ of the stellar radius.

Pinning of superfluid vortices is predicted to prevail throughout much of a neutron star. Based on the idea of Alpar et al., I develop a description of the coupling between the solid and liquid components of a neutron star through {\em thermally-activated vortex slippage}, and calculate the the response to a spin glitch. The treatment begins with a derivation of the vortex velocity from the vorticity equations of motion. The activation energy for vortex slippage is obtained from a detailed study of the mechanics and energetics of vortex motion. I show that the "linear creep" regime introduced by Alpar et al. and invoked in fits to post-glitch response is not realized for physically reasonable parameters, a conclusion that strongly constrains the physics of post-glitch response through thermal activation. Moreover, a regime of "superweak pinning", crucial to the theory of Alpar et al. and its extensions, is probably precluded by thermal fluctuations. The theory given here has a robust conclusion that can be tested by observations: {\em for a glitch in spin rate of magnitude $\Delta\nu$, pinning introduces a delay in the post-glitch response time}. The delay time is $t_d=7 (t_{sd}/10^4\mbox{yr})((\Delta\nu/\nu)/10^{-6})$ d where $t_{sd}$ is the spin-down age; $t_d$ is typically weeks for the Vela pulsar and months in older pulsars, and is independent of the details of vortex pinning. Post-glitch response through thermal activation cannot occur more quickly than this timescale. Quicker components of post-glitch response as have been observed in some pulsars, notably, the Vela pulsar, cannot be due to thermally-activated vortex motion but must represent a different process, such as drag on vortices in regions where there is no pinning. I also derive the mutual friction force for a pinned superfluid at finite temperature for use in other studies of neutron star hydrodynamics.

Sudden relaxation of the magnetic field in the core of a magnetar produces mechanical energy primarily in the form of shear waves which propagate to the surface and enter the magnetosphere as relativistic Alfv\'en waves. Due to a strong impedance mismatch, shear waves excited in the star suffer many reflections before exiting the star. If mechanical energy is deposited in the core and is converted {\em directly} to radiation upon propagation to the surface, the rise time of the emission is at least seconds to minutes, and probably minutes to hours for a realistic magnetic field geometry, at odds with observed rise times of $\lap 10$ ms for both and giant flares. Mechanisms for both small and giant flares that rely on the sudden relaxation of the magnetic field of the core are rendered unviable by the impedance mismatch, requiring the energy that drives these events to be stored in the magnetosphere just before the flare. ends, unless the waves are quickly damped.

The superfluid core of a neutron star is usually assumed to corotate with the crust over timescales longer than minutes. I show that the interaction between the neutron superfluid and the type II superconductor of the outer core increases the coupling time to weeks or longer. I suggest that observed post-glitch response over timescales of weeks to years represents recovery of the outer core. Spin glitches could originate in either the inner crust or the outer core.

Shear-driven turbulence in the superfluid interior of a neutron star exerts a fluctuating torque on the rigid crust, causing the rotational phase to walk randomly. The phase fluctuation spectrum is calculated analytically for incompressible Kolmogorov turbulence and is found to be red; the half-power point is set by the observed spin-down rate, the crust-superfluid lag, and the dynamical response time of the superfluid. Preliminary limits are placed on the latter quantities using selected time- and frequency-domain data. It is found that measurements of the normalization and slope of the power spectrum are reproduced for reasonable choices of the turbulence parameters. The results point preferentially to the neutron star interior containing a turbulent superfluid rather than a turbulent Navier-Stokes fluid. The implications for gravitational wave detection by pulsar timing arrays are discussed briefly.

We study statistical properties of stochastic variations in pulse arrival times, timing noise, in radio pulsars using a new analysis method applied in the time domain. The method proceeds in two steps. First, we subtract low-frequency wander using a high-pass filter. Second, we calculate the discrete correlation function of the filtered data. As a complementary method for measuring correlations, we introduce a statistic that measures the dispersion of the data with respect to the data translated in time. The analysis methods presented here are robust and of general usefulness for studying arrival time variations over timescales approaching the average sampling interval. We apply these methods to timing data for 32 pulsars. In two radio pulsars, PSRs B1133+16 and B1933+16, we find that fluctuations in arrival times are correlated over timescales of 10 - 20 d with the distinct signature of a relaxation process. Though this relaxation response could be magnetospheric in origin, we argue that damping between the neutron star crust and interior liquid is a more likely explanation. Under this interpretation, our results provide the first evidence independent from pulsar spin glitches of differential rotation in neutron stars. PSR B0950+08, shows evidence for quasi-periodic oscillations that could be related to mode switching.

Pinning of superfluid vortices to the nuclear lattice of the inner crust of a neutron star supports a velocity difference between the superfluid and the solid as the star spins down. Under the Magnus force that arises on the vortex lattice, vortices undergo {\em vortex creep} through thermal activation or quantum tunneling. We examine the hydrodynamic stability of this situation. Vortex creep introduces two low-frequency modes, one of which is unstable above a critical wavenumber for any non-zero flow velocity of the superfluid with respect to the solid. For typical pinning parameters of the inner crust, the superfluid flow is unstable over length scales $\lap 10$ m and over timescales as fast as months. The vortex lattice could degenerate into a tangle, and the superfluid flow could become turbulent. Unexpectedly large dissipation would suppress this instability.

Pinning of superfluid vortices to magnetic flux tubes in the outer core of a neutron star supports a velocity difference of $\sim 10^5$ \cms\ between the neutron superfluid and the proton-electron fluid as the star spins down. Under the Magnus force that arises on the vortex array, vortices undergo {\em vortex creep} through thermal activation or quantum tunneling. We examine the hydrodynamic stability of this situation. Vortex creep introduces two low-frequency modes, one of which is unstable above a critical wavenumber for any non-zero flow velocity of the neutron superfluid with respect to the charged fluid. For typical pinning parameters of the outer core, the superfluid flow is unstable over wavelengths $\lambda\lap 10$ m and over timescales of $\sim (\lambda/{1 m})^{1/2}$ yr down to $\sim 1$ d. The vortex lattice could degenerate into a tangle, and the superfluid flow would become turbulent. We suggest that superfluid turbulence could be responsible for the red timing noise seen in many neutron stars, and find a predicted spectrum that is generally consistent with observations.

We investigate a thermo-resistive instability in the outer crusts of magnetars wherein a perturbation in temperature increases ohmic heating. We show that magnetars of characteristic age {\tau}_{age} ~ 10^4 yr are unstable over timescales as short as days if strong current sheets are present in the outer crust. This instability could play an important role in the thermal and magnetic field evolution of magnetars, and may be related to bursting activity in magnetars.

I study the dynamics of a superfluid vortex in a random potential, as in the inner crust of a neutron star. Below a critical flow velocity of the ambient superfluid, a vortex is effectively immobilized by lattice forces even in the limit of zero dissipation. Low-velocity, translatory motion is not dynamically possible, a result with important implications for understanding neutron star precession and the dynamical properties of superfluid nuclear matter.

The characteristic physical timescales near stellar-mass compact objects are measured in milliseconds. These timescales -- the free-fall time, the fastest stable orbital period, and stellar spin periods -- encode the fundamental physical properties of compact objects: mass, radius, and angular momentum. The characteristic temperature of matter in the vicinity of neutron stars is such that the principal electromagnetic window into their realms is the X-ray band. Because of these connections to the fundamental properties of neutron stars, X-ray timing studies remain today the most direct means of probing their structure and dynamics. While current X-ray observatories have revealed many relevant and fascinating phenomena, they lack the sensitivity to fully exploit them to uncover the fundamental properties of compact objects and their extreme physics. With this white paper, we summarize and highlight the science opportunities that will accompany an order-of-magnitude improvement in X-ray timing sensitivity, a goal attainable in the coming decade.

We show the existence of a strong trend between neutron star surface temperature and the dipolar component of the magnetic field extending through three orders of field magnitude, a range that includes magnetars, radio-quiet isolated neutron stars, and many ordinary radio pulsars. We suggest that this trend can be explained by the decay of currents in the crust over a time scale of few Myr. We estimate the minimum temperature that a NS with a given magnetic field can reach in this interpretation.

Aims: To determine whether ``vortex creep'' in neutron stars, the slow motion of neutron vortices with respect to pinning sites in the core or inner crust, is consistent with observations of long-period precession. Methods: Using the concept of vortex drag, I discuss the precession dynamics of a star with imperfectly-pinned (i.e., "creeping'') vortices. Results: The precession frequency is far too high to be consistent with observations, indicating that the standard picture of the outer core (superfluid neutrons in co-existence with type II, superconducting protons) should be reconsidered. There is a slow precession mode, but it is highly over-damped and cannot complete even a single cycle. Moreover, the vortices of the inner crust must be able to move with little dissipation with respect to the solid.

Young, rapidly rotating neutron stars could accelerate protons to energies of $\sim 1$ PeV close to the stellar surface, which scatter with x-rays from the stellar surface through the $\Delta$ resonance and produce pions. The pions subsequently decay to produce muon neutrinos. We find that the energy spectrum of muon neutrinos consists of a sharp rise at $\sim 50$ TeV, corresponding to the onset of the resonance, above which the flux drops as $\epsilon_\nu^{-2}$ up to an upper-energy cut-off that is determined by either kinematics or by the maximum energy to which protons are accelerated. We predict event rates as high as 10-50 km$^{-2}$ yr$^{-1}$ from relatively young, close neutron stars. Such fluxes would be detectable by IceCube.

Young, rapidly rotating neutron stars could accelerate ions from their surface to energies of $\sim 1$ PeV. If protons reach such energies, they will produce pions (with low probability) through resonant scattering with x-rays from the stellar surface. The pions subsequently decay to produce muon neutrinos. Here we calculate the energy spectrum of muon neutrinos, and estimate the event rates at Earth. The spectrum consists of a sharp rise at $\sim 50$ TeV, corresponding to the onset of the resonance, above which the flux drops with neutrino energy as $\epsilon_\nu^{-2}$ up to an upper-energy cut-off that is determined by either kinematics or by the maximum energy to which protons are accelerated. We estimate event rates as high as 10-100 km^${-2}$ yr$^{-1}$ from some candidates, a flux that would be easily detected by IceCube. Lack of detection would allow constraints on the energetics of the poorly-understood pulsar magnetosphere.

We present an analysis of the spin-down parameters for 131 radio pulsars for which $\ddot\nu$ has been well determined. These pulsars have characteristic ages ranging from $10^{3} - 10^{8}$ yr and spin periods in the range 0.4--30 s; nearly equal numbers of pulsars have $\ddot\nu>0$ as $\ddot\nu<0$. We find a strong correlation of $\ddot\nu$ with $\dot{\nu}$, {\em independent of the sign of} $\ddot\nu$. We suggest that this trend can be accounted for by small, stochastic deviations in the spin-down torque that are directly proportional (in magnitude) to the spin-down torque.

Stairs, Lyne & Shemar have found that arrival time residuals from PSR B1828-11 vary periodically with a period of 500 days. This behavior can be accounted for by precession of the radiopulsar, an interpretation that is reinforced by the detection of variations in its pulse profile on the same timescale. Here, we model the period residuals from PSR B1828-11 in terms of precession of a triaxial rigid body. We include two contributions to the residuals: (i) the geometric effect, which arises because the times at which the pulsar emission beam points toward the observer varies with precession phase; (ii) the spindown contribution, which arises from any dependence of the spindown torque acting on the pulsar on the angle between its spin and magnetic axes. We use the data to probe numerous properties of the pulsar, most notably its shape, and the dependence of its spindown torque on the angle between its spin and magnetic axes, for which we assume a sum of a spin-aligned component (with a weight 1-a) and a dipolar component perpendicular to the magnetic beam axis (weight a), rather than the vacuum dipole torque (a=1). We find that a variety of shapes are consistent with the residuals, with a slight statistical preference for a prolate star. Moreover, a range of torque possibilities fit the data equally well, with no strong preference for the vacuum model. In the case of a prolate star we find evidence for an angle-dependent spindown torque. Our results show that the combination of geometrical and spin-down effects associated with precession can account for the principal features of PSR B1828-11's timing behavior, without fine tuning of the parameters.

Neutron stars are efficient accelerators for bringing charges up to relativistic energies. We show that if positive ions are accelerated to ~1 PeV near the surface of a young neutron star (t_age < about 10^5 yr), protons interacting with the star's radiation field will produce beamed mu neutrinos with energies of ~50 TeV that could produce the brightest neutrino sources at these energies yet proposed. These neutrinos would be coincident with the radio beam, so that if the star is detected as a radio pulsar, the neutrino beam will sweep the Earth; the star would be a ``neutrino pulsar''. Looking for muon neutrino emission from young neutron stars will provide a valuable probe of the energetics of the neutron star magnetosphere.

I show that the standard picture of the neutron star core containing coexisting neutron and proton superfluids, with the proton component forming a type II superconductor threaded by flux tubes, is inconsistent with observations of long-period (~1 yr) precession in isolated pulsars. I conclude that either the two superfluids coexist nowhere in the stellar core, or the core is a type I superconductor rather than type II. Either possibility would have interesting implications for neutron star cooling and theories of spin jumps (glitches).

I summarize the evidence for precession of isolated neutron stars and theoretical effort to understand the observations. I discuss factors that might set the precession period, describe constraints on the material properties of the crust, and conclude with a brief discussion of possible sources of stress that would deform a neutron star to the extent required.

We calculate the crustal rigidity parameter, b, of a neutron star (NS), and show that b is a factor 40 smaller than the standard estimate due to Baym & Pines (1971). For a NS with a relaxed crust, the NS's free-precession frequency is directly proportional to b. We apply our result for b to PSR 1828-11, a 2.5 Hz pulsar that appears to be precessing with period 511 d. Assuming this 511-d period is set by crustal rigidity, we show that this NS's crust is not relaxed, and that its reference spin (roughly, the spin for which the crust is most relaxed) is 40 Hz, and that the average spindown strain in the crust is 5 \times 10^{-5}. We also briefly describe the implications of our b calculation for other well-known precession candidates.

The neutron vortices thought to exist in the inner crust of a neutron star interact with nuclei and are expected to pin to the nuclear lattice. Evidence for long-period precession in pulsars, however, requires that pinning be negligible. We estimate the strength of vortex pinning and show that hydrodynamic forces present in a precessing star are likely sufficient to unpin all of the vortices of the inner crust. In the absence of precession, however, vortices could pin to the lattice with sufficient strength to explain the giant glitches observed in many radio pulsars.

Many radio pulsars exhibit glitches wherein the star's spin rate increases fractionally by $\sim 10^{-10} - 10^{-6}$. Glitches are ascribed to variable coupling between the neutron star crust and its superfluid interior. With the aim of distinguishing among different theoretical explanations for the glitch phenomenon, we study the response of a neutron star to two types of perturbations to the vortex array that exists in the superfluid interior: 1) thermal motion of vortices pinned to inner crust nuclei, initiated by sudden heating of the crust, (e.g., a starquake), and 2) mechanical motion of vortices, (e.g., from crust cracking by superfluid stresses). Both mechanisms produce acceptable fits to glitch observations in four pulsars, with the exception of the 1989 glitch in the Crab pulsar, which is best fit by the thermal excitation model. The two models make different predictions for the generation of internal heat and subsequent enhancement of surface emission. The mechanical glitch model predicts a negligible temperature increase. For a pure and highly-conductive crust, the thermal glitch model predicts a surface temperature increase of as much as $\sim$ 2%, occurring several weeks after the glitch. If the thermal conductivity of the crust is lowered by a high concentration of impurities, however, the surface temperature increases by $\sim$ 10% about a decade after a thermal glitch. A thermal glitch in an impure crust is consistent with the surface emission limits following the January 2000 glitch in the Vela pulsar. Future surface emission measurements coordinated with radio observations will constrain glitch mechanisms and the conductivity of the crust.

Pulse timing of the isolated pulsar PSR B1882-11 shows strong Fourier power at periods ~1000, 500 and 250 d, correlated with changes in the pulse profile (Stairs, Lyne & Shemar 2000). We study the extent to which these data can be explained by precession of the star's rigid crust coupled to the magnetic dipole torque. We find that the correlated changes in the pulse duration and spin period derivative can be explained as precession at a period of ~500 d with a wobble angle of ~3 deg if the star's dipole moment is nearly orthogonal to its symmetry axis. The dipole torque produces a harmonic at ~250 d. Comparison of the predicted spin dynamics with the observed pulse durations requires the radio beam to have a non-standard ``hour-glass'' shape. We make predictions of variations in beam polarization and pulse profile with which to test this interpretation. The precession interpretation of PSR B1828-11 seriously challenges the current understanding of the liquid interior of the neutron star. In particular, if the internal liquid is in a superfluid state, its rotational vortices cannot be significantly pinned to the crust.

The neutron star crust is rigid material floating on a neutron-proton liquid core. As the star's spin rate slows, the changing stellar shape stresses the crust and causes fractures. These starquakes may trigger pulsar glitches as well as the jumps in spin-down rate that are observed to persist after some glitches. Earlier studies found that starquakes in spinning-down neutron stars push matter toward the magnetic poles, causing temporary misalignment of the star's spin and angular momentum. After the star relaxes to a new equilibrium orientation, the magnetic poles are closer to the equator, and the magnetic braking torque is increased. The magnitude and sign of the predicted torque changes are in agreement with the observed persistent spin-down offsets. Here we examine the relaxation processes by which the new equilibrium orientation is reached. We find that the neutron superfluid in the inner crust slows as the star's spin realigns with the angular momentum, causing the crust to spin more rapidly. For plausible parameters the time scale and the magnitude of the crust's spin up agree with the giant glitches in the Vela and other pulsars.

With the aim of constraining the structural properties of neutron stars and the equation of state of dense matter, we study sudden spin-ups, glitches, occurring in the Vela pulsar and in six other pulsars. We present evidence that glitches represent a self-regulating instability for which the star prepares over a waiting time. The angular momentum requirements of glitches in Vela indicate that at least 1.4% of the star's moment of inertia drives these events. If glitches originate in the liquid of the inner crust, Vela's `radiation radius' $R_\infty$ must exceed ~12 km for a mass of 1.4 solar masses. The isolated neutron star RX J18563-3754 is a promising candidate for a definitive radius measurement, and offers to further our understanding of dense matter and the origin of glitches.

Gravitational, magnetic and superfluid forces can stress the crust of an evolving neutron star. Fracture of the crust under these stresses could affect the star's spin evolution and generate high-energy emission. We study the growth of strain in the crust of a spinning down, magnetized neutron star and examine the initiation of crust cracking (a {\em starquake}). In preliminary work (Link, Franco & Epstein 1998), we studied a homogeneous model of a neutron star. Here we extend this work by considering a more realistic model of a solid, homogeneous crust afloat on a liquid core. In the limits of astrophysical interest, our new results qualitatively agree with those from the simpler model: the stellar crust fractures under shear stress at the rotational equator, matter moves to higher latitudes and the star's oblateness is reduced. Magnetic stresses favor faults directed toward the magnetic poles. Thus our previous conclusions concerning the star's spin response still hold; namely, asymmetric redistribution of matter excites damped precession which could ultimately lead to an increase in the spin-down torque. Starquakes associated with glitches could explain the permanent {\em offsets} in period derivative observed to follow glitches in at least three pulsars.

With the aim of constraining the structural properties of neutron stars and the equation of state of dense matter, we study sudden spin-ups, glitches, occurring in the Vela pulsar and in six other pulsars. We present evidence that glitches represent a self-regulating instability for which the star prepares over a waiting time. The angular momentum requirements of glitches in Vela indicate that at least 1.4% of the star's moment of inertia drives these events. If glitches originate in the liquid of the inner crust, Vela's `radiation radius' must exceed ~12 km for a mass of 1.4 solar masses. Observational tests of whether other neutron stars obey this constraint will be possible in the near future.

Based on the study of temporal asymmetry of 631 gamma ray bursts from the BATSE 3B catalog by Link and Epstein [Ap J 466, 764 (1996)], we identify the population of bursts whose rising times are longer than their decays, thus showing atypical profiles. We analyse their sky distribution, morphology, time-space clustering and other average properties and compare them with those associated with the bulk of the bursts. We show how most of the peculiar bursts analysed are consistent with recent fireball models, but a fraction of bursts ($\sim 4$% of the total sample) appear to be inconsistent.

The recent temperature measurements of the two older isolated neutron stars PSR 1929+10 and PSR 0950+08 (ages of $3\times 10^6$ and $2\times 10^7$ yr, respectively) indicate that these objects are heated. A promising candidate heat source is friction between the neutron star crust and the superfluid it is thought to contain. We study the effects of superfluid friction on the long-term thermal and rotational evolution of a neutron star. Differential rotation velocities between the superfluid and the crust (averaged over the inner crust moment of inertia) of $\bar\omega\sim 0.6$ rad s$^{-1}$ for PSR 1929+10 and $\sim 0.02$ rad s$^{-1}$ for PSR 0950+08 would account for their observed temperatures. These differential velocities could be sustained by pinning of superfluid vortices to the inner crust lattice with strengths of $\sim$ 1 MeV per nucleus. Pinned vortices can creep outward through thermal fluctuations or quantum tunneling. For thermally-activated creep, the coupling between the superfluid and crust is highly sensitive to temperature. If pinning maintains large differential rotation ($\sim 10$ rad s$^{-1}$), a feedback instability could occur in stars younger than $\sim 10^5$ yr causing oscillations of the temperature and spin-down rate over a period of $\sim 0.3 t_{\rm age}$. For stars older than $\sim 10^6$ yr, however, vortex creep occurs through quantum tunneling, and the creep velocity is too insensitive to temperature for a thermal-rotational instability to occur. These older stars could be heated through a steady process of superfluid friction.

The persistent increases in spin-down rate offsets seen to accompany glitches in the Crab and other pulsars suggest increases in the spin-down torque. We interpret these offsets as due to starquakes occurring as the star spins down and the rigid crust becomes less oblate. We study the evolution of strain in the crust, the initiation of starquakes, and possible consequences for magnetic field and torque evolution. Crust cracking occurs as equatorial material shears under the compressive forces arising from the star's decreasing circumference, and matter moves to higher latitudes along a fault inclined to the equator. A starquake is most likely to originate near one of the two points on the rotational equator farthest from the magnetic poles. The material breaks along a fault approximately aligned with the magnetic poles. We suggest that the observed offsets come about when a starquake perturbs the star's mass distribution, producing a misalignment of the angular momentum and spin axes. Subsequently, damped precession to a new rotational state increases the angle alpha between the rotation and magnetic axes. The resulting increase in external torque appears as a permanent increase in the spin-down rate. Repeated starquakes would continue to increase alpha, making the pulsar more of an orthogonal rotator.

We show that the brighter gamma-ray bursts (GRBs) in the BATSE 3B catalog concentrate slightly toward the Galactic plane and center, suggesting that at least some bursts originate within the Galaxy. To develop an interpretation of this brightness-dependent anisotropy, we consider GRBs distributed in a thick disk centered on the Galaxy. As an approximation to the true luminosity distribution, we divide bursts into two distinct luminosity classes. Most bursts originate from low-luminosity, nearby sources while the relatively few high-luminosity sources trace the shape of the thick disk. We find that characteristic disk dimensions as small as 20 kpc in thickness and 30 kpc in radial extent can match the observed brightness-dependent anisotropy as well as the number-flux distribution; these dimensions are smaller than previously considered viable. For bursts distributed in a disk of these compact dimensions, the low-luminosity sources must have peak powers of $\lapprox 5\times 10^{40}$ erg s$^{-1}$ and rate densities of $\gapprox 5\times 10^{-3}$ yr$^{-1}$ kpc$^{-3}$. The high-luminosity sources must have a peak luminosity of $\simeq 10^{42}$ erg s$^{-1}$ but their rate density need be only 5% or less than that of the low-luminosity sources. If the Andromeda Galaxy contained a similar population of high-luminosity sources, their peak gamma-ray fluxes would be less than $\sim 0.1$ s$^{-1}$ cm$^{-2}$, making their detection problematic. Statistical analyses of the gamma-ray burst events obtained since the compilation of the 3B catalog will determine whether the thick-disk interpretation is preferred over the isotropic one.

Variation in the angle $\alpha$ between a pulsar's rotational and magnetic axes would change the torque and spin-down rate. We show that sudden increases in $\alpha$, coincident with glitches, could be responsible for the persistent increases in spin-down rate that follow glitches in the Crab pulsar. Moreover, changes in $\alpha$ at a rate similar to that inferred for the Crab pulsar account naturally for the very low braking index of the Vela pulsar. If $\alpha$ increases with time, all pulsar ages obtained from the conventional braking model are underestimates. Decoupling of the neutron star liquid interior from the external torque cannot account for Vela's low braking index. Variations in the Crab's pulse profile due to changes in $\alpha$ might be measurable.

We study the temporal asymmetry of over 600 bursts from the BATSE 3B catalog, encompassing a 200-fold range in peak flux. By comparing the rates of rise and fall of the flux near the highest burst peak, we find that about two-thirds of the bursts exhibit a preferred asymmetry in the sense that the flux rises more rapidly than it falls, confirming the conclusions of previous studies employing smaller databases. The statistical significance of the average time asymmetry of the sample is $>99.999$%; therefore, models that predict time symmetry of the burst profile are ruled out. We find no statistically significant correlation between burst temporal asymmetry and peak. This result is consistent with both cosmological and local interpretations of the gamma ray burst phenomenon.

We examine the thermal and dynamical response of a neutron star to a sudden perturbation of the inner crust temperature. During the star's evolution, starquakes and other processes may deposit $\gap 10^{42}$ ergs, causing significant internal heating and increased frictional coupling between the crust and the more rapidly rotating neutron superfluid the star is expected to contain. Through numerical simulation we study the propagation of the thermal wave created by the energy deposition, the induced motion of the interior superfluid, and the resulting spin evolution of the crust. We find that energy depositions of $\sim 10^{40}$ ergs produce gradual spin-ups above the timing noise level, while larger energy depositions produce sudden spin jumps resembling pulsar glitches. For a star with a temperature in the observed range of the Vela pulsar, an energy deposition of $\sim 10^{42}$ ergs produces a large spin-up taking place over minutes, similar to the Vela ``Christmas'' glitch. Comparable energy deposition in a younger and hotter ``Crab-like'' star produces a smaller spin-up taking place over $\sim 1$ day, similar to that seen during the partially time-resolved Crab glitch of 1989.

Differential rotation between the neutron star crust and a more rapidly rotating interior superfluid leads to frictional heating that affects the star's long-term thermal evolution and resulting surface emission. Here we present the results of thermal evolution simulations based on two models of the vortex pinning forces that bracket a range of plausible pinning strengths. We include the effects of superfluidity, magnetic fields, and temperature gradients. As representative standard and accelerated neutrino emission processes taking place in the core, we consider the modified Urca process in normal baryonic matter, and the much faster quark Urca process. Comparison of our results with neutron star surface temperature data, including the recent temperature measurement of the Geminga pulsar, shows that stars with soft equations of state and modest frictional heating are in closest agreement with the data; stars with stronger frictional heating have temperatures inconsistent with the upper limit of PSR 1929+10. Stiffer stars undergoing standard cooling generally have temperatures lying above the Vela detection, a situation worsened by the inclusion of frictional heating. Stars undergoing accelerated cooling without frictional heating have temperatures that fall far below most temperature measurements; the Vela and Geminga detections being the most compelling examples. Only in stiff stars, which have thick crusts, can the inclusion of strong frictional heating raise the temperature at late stages in the evolution to a level consistent with the data. However, such a large amount of heating leads to a temperature at $\sim 1000$ yr in excess of the Crab upper limit. Suppression of accelerated neutrino emission processes, perhaps by superfluid pairing in the core, may yield acceptable cooling models.