At second order in perturbation theory, the $r$-modes of uniformly rotating stars include an axisymmetric part that can be identified with differential rotation of the background star. If one does not include radiation-reaction, the differential rotation is constant in time and has been computed by S\'a. It has a gauge dependence associated with the family of time-independent perturbations that add differential rotation to the unperturbed equilibrium star: For stars with a barotropic equation of state, one can add to the time-independent second-order solution arbitrary differential rotation that is stratified on cylinders (that is a function of distance $\varpi$ to the axis of rotation). We show here that the gravitational radiation-reaction force that drives the $r$-mode instability removes this gauge freedom: The expontially growing differential rotation of the unstable second-order $r$-mode is unique. We derive a general expression for this rotation law for Newtonian models and evaluate it explicitly for slowly rotating models with polytropic equations of state.

We study the inertial modes of slowly rotating, fully relativistic compact stars. The equations that govern perturbations of both barotropic and non-barotropic models are discussed, but we present numerical results only for the barotropic case. For barotropic stars all inertial modes are a hybrid mixture of axial and polar perturbations. We use a spectral method to solve for such modes of various polytropic models. Our main attention is on modes that can be driven unstable by the emission of gravitational waves. Hence, we calculate the gravitational-wave growth timescale for these unstable modes and compare the results to previous estimates obtained in Newtonian gravity (i.e. using post-Newtonian radiation formulas). We find that the inertial modes are slightly stabilized by relativistic effects, but that previous conclusions concerning eg. the unstable r-modes remain essentially unaltered when the problem is studied in full general relativity.

Several recent surprises appear dramatically to have improved the likelihood that the spin of rapidly rotating, newly formed neutron stars (and, possibly, of old stars spun up by accretion) is limited by a nonaxisymmetric instability driven by gravitational waves. Except for the earliest part of the spin-down, the axial l=m=2 mode (an r-mode) dominates the instability, and the emitted waves may be observable by detectors with the sensitivity of LIGO II. A review of these hopeful results is followed by a discussion of constraints on the instability set by dissipative mechanisms, including viscosity, nonlinear saturation, and energy loss to a magnetic field driven by differential rotation.

We study the r-modes and rotational ``hybrid'' modes of relativistic stars. As in Newtonian gravity, the spectrum of low-frequency rotational modes is highly sensitive to the stellar equation of state. If the star and its perturbations obey the same one-parameter equation of state (as with isentropic stars), there exist {\it no pure r-modes at all} - no modes whose limit, for a star with zero angular velocity, is an axial-parity perturbation. Rotating stars of this kind similarly have no pure g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed pressure and density. We compute the post-Newtonian corrections to the $l=m$ r-modes of isentropic and non-isentropic uniform density stars.

Despite the recent excitement over the r-mode instability of rotating stars, these modes are not yet well-understood for stellar models appropriate to neutron stars - perfect fluid models in which both the equilibrium and perturbed configurations obey the same one-parameter equation of state. In spherical stars of this kind, r-modes and g-modes form a degenerate zero-frequency subspace. Rotation splits the degeneracy to zeroth order in the star's angular velocity $\Omega$, and the resulting modes are generically hybrids, whose limit as $\Omega\to 0$ is a stationary current with axial and polar parts. Lindblom and Ipser have recently found these hybrid rotational modes in an analytic study of the Maclaurin spheroids. This dissertation studies them in Newtonian stars (both uniform density models and $n=1$ polytropes) and reports a first computation of the hybrid rotational modes of relativistic stars.

A brief review of the stability of rotating relativistic stars is followed by a more detailed discussion of recent work on an instability of r-modes, modes of rotating stars that have axial parity in the slow-rotation limit. These modes may dominate the spin-down of neutron stars that are rapidly rotating at birth, and the gravitational waves they emit may be detectable.

Almost none of the r-modes ordinarily found in rotating stars exist, if the star and its perturbations obey the same one-parameter equation of state; and rotating relativistic stars with one-parameter equations of state have no pure r-modes at all, no modes whose limit, for a star with zero angular velocity, is a perturbation with axial parity. Similarly (as we show here) rotating stars of this kind have no pure g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed pressure and density. Where have these modes gone? In spherical stars of this kind, r-modes and g-modes form a degenerate zero-frequency subspace. We find that rotation splits the degeneracy to zeroth order in the star's angular velocity $\Omega$, and the resulting modes are generically hybrids, whose limit as $\Omega\to 0$ is a stationary current with axial and polar parts. Because each mode has definite parity, its axial and polar parts have alternating values of $l$. We show that each mode belongs to one of two classes, axial-led or polar-led, depending on whether the spherical harmonic with lowest value of $l$ that contributes to its velocity field is axial or polar. We numerically compute these modes for slowly rotating polytropes and for Maclaurin spheroids, using a straightforward method that appears to be novel and robust. Timescales for the gravitational-wave driven instability and for viscous damping are computed using assumptions appropriate to neutron stars.