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.

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 long-held notion that the highest-energy cosmic rays are of distant extragalactic origin is challenged by observations that events above $\sim 10^{20}$ eV do not exhibit the expected high-energy cutoff from photopion production off the cosmic microwave background. We suggest that these unexpected ultra-high-energy events are due to iron nuclei accelerated from young strongly magnetized neutron stars through relativistic MHD winds. We find that neutron stars whose initial spin periods are shorter than $\sim 4 (B_S/10^{13}{\rm G})^{1/2}$ ms, where $B_S$ is the surface magnetic field, can accelerate iron cosmic rays to greater than $\sim 10^{20}$ eV. These ions can pass through the remnant of the supernova explosion that produced the neutron star without suffering significant spallation reactions. For plausible models of the Galactic magnetic field, the trajectories of the iron ions curve sufficiently to be consistent with the observed arrival directions of the highest energy events.

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.

The absence of the expected GZK cutoff strongly challenges the notion that the highest-energy cosmic rays are of distant extragalactic origin. We discuss the possibility that these ultra-high-energy events originate in our Galaxy and propose that they may be due to iron nuclei accelerated from young, strongly magnetic neutron stars. Newly formed pulsars accelerate ions from their surface through relativistic MHD winds. We find that pulsars whose initial spin periods are shorter than $\sim 4 (B_S/10^{13}{\rm G})$ ms, where $B_S$ is the surface magnetic field, can accelerate iron ions to greater than $10^{20} eV$. These ions can pass through the remnant of the supernova explosion that produced the pulsar without suffering significant spallation reactions. Depending on the structure of the galactic magnetic field, the trajectories of the iron ions from galactic sources can be consistent with the observed arrival directions of the highest energy events.

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 show that the short spin-up time observed for the Vela pulsar during the 1988 ``Christmas'' glitch implies that the coupling time of the pulsar core to its crust is less than $\sim$ 10 seconds. Ekman pumping cannot explain the fast core-crust coupling and a more effective coupling is necessary. The internal magnetic field of the Vela pulsar can provide the necessary coupling if the field threads the core with a magnitude that exceeds $ 10^{13}$ Gauss for a normal interior and $ 10^{11}$ Gauss for a superconducting interior. These lower bounds favor the hypothesis that the interior of neutron stars contains superfluid neutrons and protons and challenge the notion that pulsar magnetic fields decay over million year time scales or that magnetic flux is expelled from the core as the star slows.

We examine the dynamics of a rotating viscous fluid following an abrupt change in the angular velocity of the solid bounding surface. We include the effects of a density stratification and compressibility which are important in astrophysical objects such as neutron stars. We confirm and extend the conclusions of previous studies that stratification restricts the Ekman pumping process to a relatively thin layer near the boundary, leaving much of the interior fluid unaffected. We find that finite compressibility further inhibits Ekman pumping by decreasing the extent of the pumped layer and by increasing the time for spin-up. The results of this paper are important for interpreting the spin period discontinuities (``glitches'') observed in rotating neutron stars.

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.