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.

The linear response of a superfluid, rotating uniformly in a cylindrical container and threaded with a large number of vortex lines, to an impulsive increase in the angular velocity of the container is investigated. At zero temperature and with perfect pinning of vortices to the top and bottom of the container, we demonstrate that the system oscillates persistently with a frequency proportional to the vortex line tension parameter to the quarter power. This low-frequency mode is generated by a secondary flow analogous to classical Ekman pumping that is periodically reversed by the vortex tension in the boundary layers. We compare analytic solutions to the two-fluid equations of Chandler & Baym (1986) with the spin-up experiments of Tsakadze & Tsakadze (1980) in helium II and find the frequency agrees within a factor of four, although the experiment is not perfectly suited to the application of the linear theory. We argue that this oscillatory Ekman pumping mode, and not Tkachenko modes provide a natural explanation for the observed oscillation. In neutron stars, the oscillation period depends on the pinning interaction between neutron vortices and flux tubes in the outer core. Using a simplified pinning model, we demonstrate that strong pinning can accommodate modes with periods of days to years, which are only weakly damped by mutual friction over longer timescales.

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.

Following a glitch, the crust and magnetized plasma in the outer core of a neutron star are believed to rapidly establish a state of co-rotation within a few seconds by process analogous to classical Ekman pumping. However, in ideal magnetohydrodynamics, a final state of co-rotation is inconsistent with conservation of energy of the system. We demonstrate that, after the Ekman-like spin up is completed, magneto-inertial waves continue to propagate throughout the star, exciting torsional oscillations in the crust and plasma. The crust oscillation is irregular and quasi-periodic, with a dominant frequency of the order of seconds. Crust oscillations commence after an Alfv\'en crossing time, approximately half a minute at the magnetic pole, and are subsequently damped by the electron viscosity over approximately an hour. In rapidly rotating stars, the magneto-inertial spectrum in the core approaches a continuum, and crust oscillations are damped by resonant absorption analogous to quasi-periodic oscillations in magnetars. The oscillations predicted are unlikely to be observed in timing data from existing radio telescopes, but may be visible to next generation telescope arrays.

A recipe is presented for interpreting non-invasively the transport processes at work during relaxation of a cylindrical, superfluid-filled vessel, after it is accelerated impulsively and then allowed to respond to the viscous torque exerted by the contained fluid. The recipe exploits a recently published analytic solution for Ekman pumping in a two-component superfluid, which treats the back-reaction self-consistently in arbitrary geometry for the first time. The applicability of the recipe to He II, 3He, 3He-4He mixtures and Bose-Einstein condensates is assessed, and the effects of turbulence discussed.

The spin up of helium II is studied by calculating the spin-down recovery of a superfluid-filled container after an impulsive acceleration and comparing with experiments. The calculation takes advantage of a recently published analytic solution for the spin up of a Hall-Vinen-Bekharevich-Khalatnikov superfluid that treats the back-reaction torque exerted by the viscous component self-consistently in arbitrary geometry for the first time. Excellent agreement at the 0.5% level is obtained for experiments at $T=1.57\,{\rm K}$, after correcting for the non-uniform rotation in the initial state, confirming that vortex tension and pinning (which are omitted from the theory) play a minimal role under certain conditions (small Rossby number, smooth walls). The dependence of the spin-down time on temperature and the mass fraction of the viscous component are also investigated. Closer to the lambda point, the predicted onset of turbulence invalidates the linear Ekman theory.

Nonaxisymmetric, meridional circulation inside a neutron star, excited by a glitch and persisting throughout the post-glitch relaxation phase, emits gravitational radiation. Here, it is shown that the current quadrupole contributes more strongly to the gravitational wave signal than the mass quadrupole evaluated in previous work. We calculate the signal-to-noise ratio for a coherent search and conclude that a large glitch may be detectable by second-generation interferometers like the Laser Interferometer Gravitational-Wave Observatory. It is shown that the viscosity and compressibility of bulk nuclear matter, as well as the stratification length-scale and inclination angle of the star, can be inferred from a gravitational wave detection in principle.

A two-component hydrodynamic model is constructed of the global superfluid flow induced by two-component Ekman pumping during the recovery stage of a glitch. The model successfully accounts for the quasi-exponential recovery observed in pulsars like Vela and the "overshoot" observed in pulsars like the Crab. By fitting the model to high-resolution timing data, three important constitutive coefficients in bulk nuclear matter can be extracted: the shear viscosity, the mutual friction parameter, and the charged fluid fraction. The fitted coefficients for the Crab and Vela are compared with theoretical predictions for several equations of state, including the color-flavor locked and two-flavor color superconductor phases of quark matter.