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.

(Abridged). The type-I X-ray bursting low mass X-ray binary KS 1731-260 was recently detected for the first time in quiescence by Wijnands et al., following an approximately 13 yr outburst which ended in Feb 2001. Unlike all other known transient neutron stars, the duration of this recent outburst is as long as the thermal diffusion time of the crust. The large amount of heat deposited by reactions in the crust will have heated the crust to temperatures much higher than the equilibrium core temperature. As a result, the thermal luminosity currently observed from the neutron star is dominated not by the core, but by the crust. Moreover, the level and the time evolution of quiescent luminosity is determined mostly by the amount of heat deposited in the crust during the most recent outburst. Using estimates of the outburst mass accretion rate, our calculations of the quiescent flux immediately following the end of the outburst agree with the observed quiescent flux to within a factor of a few. We present simulations of the evolution of the quiescent lightcurve for different scenarios of the crust microphysics, and demonstrate that monitoring observations (with currently flying instruments) spanning from 1--30 yr can measure the crust cooling timescale and the total amount of heat stored in the crust. These quantities have not been directly measured for any neutron star.

We analyze the global hydrodynamic flow in the ocean of an accreting, rapidly rotating, non-magnetic neutron star in an LMXB during a type I X-ray burst. Our analysis takes into account the rapid rotation of the star and the lift-up of the burning ocean during the burst. We find a new regime for spreading of a nuclear burning front, where the flame is carried along a coherent shear flow across the front. If turbulent viscosity is weak, the speed of flame propagation is ~20 km/s, while, if turbulent viscosity is dynamically important, the flame speed increases, and reaches the maximum value, ~300 km/s, when the eddy overturn frequency is comparable to the Coriolis parameter. We show that, due to rotationally reduced gravity, the thermonuclear runaway is likely to begin on the equator. The equatorial belt is ignited first, and the flame then propagates from the equator to the poles. Inhomogeneous cooling (equator first, poles second) drives strong zonal currents which may be unstable to formation of Jupiter-type vortices; we conjecture that these vortices are responsible for modulation of X-ray flux in the tails of some bursts. We consider the effect of strong zonal currents on the frequency of modulation of the X-ray flux and show that the large values of the frequency drifts observed in some bursts can be accounted for within our model combined with the model of homogeneous radial expansion. Also, if inhomogeneities are trapped in the forward zonal flows around the propagating burning front, chirps with large frequency ranges (~25-500 Hz) may be detectable during the burst rise. We argue that an MHD dynamo within the burning front can generate a small-scale magnetic field, which may enforce vertically rigid flow in the front's wake and can explain the coherence of oscillations in the burst tail.

Transiently accreting neutron stars in quiescence (Lx<10^34 erg/s) have been observed to vary in intensity by factors of few, over timescales of days to years. If the quiescent luminosity is powered by a hot NS core, the core cooling timescale is much longer than the recurrence time, and cannot explain the observed, more rapid variability. However, the non-equilibrium reactions which occur in the crust during outbursts deposit energy in iso-density shells, from which the thermal diffusion timescale to the photosphere is days to years. The predicted magnitude of variability is too low to explain the observed variability unless - as is widely believed - the neutrons beyond the neutron-drip density are superfluid. Even then, variability due to this mechanism in models with standard core neutrino cooling processes is less than 50 per cent - still too low to explain the reported variability. However, models with rapid core neutrino cooling can produce variability by a factor as great as 20, on timescales of days to years following an outburst. Thus, the factors of few intensity variability observed from transiently accreting neutron stars can be accounted for by this mechanism only if rapid core cooling processes are active.

Motivated by the narrow range of spin frequencies of nearly 20 accreting neutron stars, Bildsten (1998) conjectured that their spin-up had been halted by the emission of gravitational waves. He also pointed out that small nonaxisymmetric temperature variations in the accreted crust will lead to "wavy" electron capture layers, whose horizontal density variations naturally create a mass quadrupole moment. We present a full calculation of the crust's elastic adjustment to these density perturbations and find that the elastic response of the crust reduces Bildsten's original estimate of the quadrupole moment in the thin outer crust by a factor of 20-50. However, this basic picture, when applied to capture layers in the deep inner crust, can generate quadrupoles in the necessary range as long as there are ~5% lateral temperature variations in the inner crust. By calculating the thermal flow throughout the core and the crust, we find that temperature gradients this large are easily maintained by asymmetric heat sources or lateral composition gradients in the crust. We also derive a general relation between the stresses and strains in the crust and the maximum quadrupole moment they can generate. We show under quite general conditions that maintaining a quadrupole of the magnitude necessary to balance the accretion torque requires dimensionless strains close to 0.01 at near-Eddington accretion rates, of order the breaking strain of conventional materials.

The presence of a viscous boundary layer under the solid crust of a neutron star dramatically increases the viscous damping rate of the fluid r-modes. We improve previous estimates of this damping rate by including the effect of the Coriolis force on the boundary-layer eigenfunction and by using more realistic neutron-star models. If the crust is assumed to be perfectly rigid, the gravitational radiation driven instability in the r-modes is completely suppressed in neutron stars colder than about 1.5 x 10^8 K. Energy generation in the boundary layer will heat the star, and will even melt the crust if the amplitude of the r-mode is large enough. We solve the heat equation explicitly (including the effects of thermal conduction and neutrino emission) and find that the r-mode amplitude needed to melt the crust is approximately a_c = 5 x 10^{-3} for maximally rotating neutron stars. If the r-mode saturates at an amplitude larger than a_c, the heat generated is sufficient to maintain the outer layers of the star in a mixed fluid-solid state analogous to the pack ice on the fringes of the Arctic Ocean. We argue that in young, rapidly rotating neutron stars this effect considerably delays the formation of the crust. By considering the dissipation in the ice flow, we show that the final spin frequency of stars with r-mode amplitude of order unity is close to the value estimated for fluid stars without a crust.

R-modes in neutron stars with crusts are damped by viscous friction at the crust-core boundary. The magnitude of this damping, evaluated by Bildsten and Ushomirsky (BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency for a neutron star spun up by accretion in a Low-Mass X-ray binary (LMXB). In this paper we explore the mechanical coupling between the core r-modes and the elastic crust, using a toy model of a constant density neutron star with a constant shear modulus crust. We find that, at spin frequencies in excess of ~50 Hz, the r-modes strongly penetrate the crust. This reduces the relative motion (slippage) between the crust and the core compared to the rigid crust limit. We therefore revise down, by as much as a factor of 10^2-10^3, the damping rate computed by BU, significantly reducing the maximal possible spin frequency of neutron star with a solid crust. The dependence of the crust-core slippage on the spin frequency is complicated, and is very sensitive to the physical thickness of the crust. If the crust is sufficiently thick, the curve of the critical spin frequency for the onset of the r-mode instability becomes multi-valued for some temperatures; this is related to the avoided crossings between the r-mode and the higher-order torsional modes in the crust. The critical frequencies are comparable to the observed spins of neutron stars in LMXBs and millisecond pulsars.

We summarize the observations of the spin periods of rapidly accreting neutron stars. If gravitational radiation is responsible for balancing the accretion torque at the observed spin frequencies of ~300 Hz, then the brightest of these systems make excellent gravitational wave sources for LIGO-II and beyond. We review the recent theoretical progress on two mechanisms for gravitational wave emission: mass quadrupole radiation from deformed neutron star crusts and current quadrupole radiation from r-mode pulsations in neutron star cores.

Recent observations suggest that neutron stars in low-mass X-ray binaries rotate within a narrow range of spin frequencies clustered around 300 Hz. A proposed explanation for this remarkable fact is that gravitational radiation from a steady-state r-mode oscillation in the neutron star's core halts the spin-up due to accretion. For the neutron star transients, balancing the time-averaged accretion torque with gravitational wave emission from steady-state, constant amplitude r-mode pulsations implies a quiescent luminosity too bright to be consistent with observations (in particular of Aql X-1). The viscous dissipation (roughly 10 MeV per accreted nucleon for a spin of 300 Hz) from such an r-mode makes the core sufficiently hot to power a thermal luminosity of order 1e34 erg/s when accretion halts. This is the minimum quiescent luminosity that the neutron star must emit when viscous heating in the core is balanced by radiative cooling from the surface, as is the case when the core of the star is superfluid. We therefore conclude that either the accretion torque is much less than (dM/dt)(GMR)^{1/2}, or that a steady-state r-mode does not limit the spin rate of the neutron star transients. Future observations with AXAF and XMM promise to further constrain the amount of viscous dissipation in the neutron star core.

We use a simple physical model to study the nonlinear behaviour of the r-mode instability. We assume that r-modes (Rossby waves) are excited in a thin spherical shell of rotating incompressible fluid. For this case, exact Rossby wave solutions of arbitrary amplitude are known. We find that: (a) These nonlinear Rossby waves carry ZERO physical angular momentum and positive physical energy, which is contrary to the folklore belief that the r-mode angular momentum and energy are negative. (b) Within our model, we confirm the differential drift reported by Rezzolla, Lamb and Shapiro (1999). Radiation reaction is introduced into the model by assuming that the fluid is electrically charged; r-modes are coupled to electromagnetic radiation through current (magnetic) multipole moments. We find that: (c) To linear order in the mode amplitude, r-modes are subject to the CFS instability, as expected. (d) Radiation reaction decreases the angular velocity of the shell and causes differential rotation (which is distinct from but similar in magnitude to the differential drift reported by Rezzolla et al.) prior to saturation of the r-mode growth. This is contrary to the phenomenological treatments to date, which assume that the loss of stellar angular momentum is accounted for by the r-mode growth. We demonstrate, for the first time, that r-mode radiation reaction leads to differential rotation. (e) We show that for l=2 r-mode electromagnetic radiation reaction is equivalent to gravitational radiation reaction in the lowest post-Newtonian order.

Recent work has raised the exciting possibility that r-modes (Rossby waves) in rotating neutron star cores might be strong gravitational wave sources. We estimate the effect of a solid crust on their viscous damping rate and show that the dissipation rate in the viscous boundary layer between the oscillating fluid and the nearly static crust is >10^5 times higher than that from the shear throughout the interior. This increases the minimum frequency for the onset of the gravitational r-mode instability to at least 500 Hz when the core temperature is less than 10^10 K. It eliminates the conflict of the r-mode instability with the accretion-driven spin-up scenario for millisecond radio pulsars and makes it unlikely that the r-mode instability is active in accreting neutron stars. For newborn neutron stars, the formation of a solid crust shortly after birth affects their gravitational wave spin-down and hence detectability by ground-based interferometric gravitational wave detectors.

Several classes of stars (most notably O and B main-sequence stars, as well as accreting white dwarfs and neutron stars) rotate quite rapidly, at spin frequencies greater than the typical g-mode frequencies. We discuss how rapid rotation modifies the $\kappa$-mechanism excitation and observability of g-mode oscillations. We find that, by affecting the timescale match between the mode period and the thermal time at the driving zone, rapid rotation stabilizes some of the g-modes that are excited in a non-rotating star, and, conversely, excites g-modes that are damped in absence of rotation. The fluid velocities and temperature perturbations are strongly concentrated near the equator for most g-modes in rapidly rotating stars, which means that a favorable viewing angle may be required to observe the pulsations. Moreover, the stability of modes of the same $l$ but different $m$ is affected differently by rotation. We illustrate this by considering g-modes in Slowly Pulsating B-type stars as a function of the rotation rate.

Observations of quasi-periodic oscillations (QPOs) in the luminosity from many accreting neutron stars (NS) have led us to investigate a source of periodicity prevalent in other stars: non-radial oscillations. After summarizing the structure of the atmosphere and ocean of an accreting NS, we discuss the various low l g-modes with frequencies in the 1-100 Hz range. Successful identification of a non-radial mode with an observed frequency would yield new information about the thermal and compositional makeup of the NS, as well as its radius. We close by discussing how rapid rotation changes the g-mode frequencies.

We present an analytic calculation of the thermonuclear depletion of the light elements lithium, beryllium, and boron in fully convective, low-mass stars. Under the presumption that the pre--main-sequence star is always fully mixed during contraction, we find that the burning of these rare light elements can be computed analytically, even when the star is degenerate. Using the effective temperature as a free parameter, we constrain the properties of low-mass stars from observational data, independently of the uncertainties associated with modeling their atmospheres and convection. Our analytic solution explains the dependence of the age at a given level of elemental depletion on the stellar effective temperature, nuclear cross sections, and chemical composition. Most importantly, our results allow observers to translate lithium non-detections in young cluster members into a model-independent minimum age for that cluster. Using this procedure, we have found lower limits to the ages of the Pleiades (100 Myr) and Alpha Persei (60 Myr) clusters. Recent experimental work on the low energy resonance in the ^10B(p,\alpha)^7Be reaction has greatly enhanced estimates of the destruction rate of ^10B, making it possible for stars with M>0.1 M_sun to deplete both ^10B and ^11B before reaching the main sequence. Moreover, there is an interesting range of masses, 0.085 M_sun < M < 0.13 M_sun, where boron depletion occurs on the main sequence in less than a Hubble time, providing a potential ``clock'' for dating low-mass stars.

We present an analytic calculation of the thermonuclear depletion of lithium in contracting, fully convective, pre-main sequence stars of mass M < 0.5 M_sun. Previous numerical work relies on still-uncertain physics (atmospheric opacities and convection, in particular) to calculate the effective temperature as a unique function of stellar mass. We assume that the star's effective temperature, T_eff, is fixed during Hayashi contraction and allow its actual value to be a free parameter constrained by observation. Using this approximation, we compute lithium burning analytically and explore the dependence of lithium depletion on T_eff, M, and composition. Our calculations yield the radius, age, and luminosity of a pre-main sequence star as a function of lithium depletion. This allows for more direct comparisons to observations of lithium depleted stars. Our results agree with those numerical calculations that explicitly determine stellar structure during Hayashi contraction. In agreement with Basri, Marcy, and Graham (1996), we show that the absence of lithium in the Pleiades star HHJ 3 implies that it is older than 100 Myr. We also suggest a generalized method for dating galactic clusters younger than 100 Myr (i.e., those with contracting stars of M > 0.08 M_sun) and for constraining the masses of lithium depleted stars.