Tkachenko modes in rotating neutron stars: the effect of compressibility and implications for pulsar timing noise

Long wavelength oscillations (Tkachenko waves) of the triangular lattice of quantized vortices in superfluid neutron stars have been suggested as one of the possible explanations for the timing noise observed in many radio pulsars, in particular for the 100-1000 day variations in the spin of PSR B1828-11. Most studies to date have, however, been based on the hydrodynamics developed for superfluid Helium. In this paper we extend the formulation to a two fluid neutron and proton system, relevant for neutron star interiors and include the effect of chemical coupling, compressibility and mutual friction between the components. In particular we find that chemical coupling and compressibility can have a drastic effect on the mode structure. However, for the slower pulsars rotating at 1-10 Hz (such as PSR 1828-11), most choices of parameters in the equation of state lead to Tkachenko oscillations with frequencies in the correct range to explain the timing noise. We also investigate the case of more rapidly rotating pulsars (above 100 Hz) for which we find that there is a vast portion of parameter space in which there are no Tkachenko modes, but only modified sound waves at much higher frequencies

Are gravitational waves spinning down PSR J1023+0038?

The pulsar J1203+0038 rotates with a frequency $\nu\approx 592$ Hz and has been observed to transition between a radio state, during which it is visible as a millisecond radio pulsar, and and a Low Mass X-ray Binary state, during which accretion powered X-ray pulsations are visible. Timing during the two phases reveals that during the LMXB phase the neutron star is spinning down at a rate of $\dot{\nu}\approx -3 \times 10^{-15}$ Hz/s, which is approximately 27\% faster than the rate measured during the radio phase, $\dot{\nu}\approx -2.4 \times 10^{-15}$ Hz/s, and at odds with the predictions of accretion models. In this letter we suggest that the increase in spin-down rate is compatible with gravitational wave emission, and in particular to the creation of a `mountain' during the accretion phase. We show that asymmetries in pycno-nuclear reaction rates in the crust can lead to a large enough mass quadrupole to explain the observed spin-down rate, which so far has no other self-consistent explanation, and that radio timing at the onset of the next millisecond radio pulsar phase can test this scenario. Another possibility is that an unstable $r$-mode with amplitude $\alpha\approx 5\times10^{-8}$ may be present in the system.Comment: 5 pages, submitted to PR

A new mechanism for saturating unstable r-modes in neutron stars

We consider a new mechanism for damping the oscillations of a mature neutron star. The new dissipation channel arises if superfluid vortices are forced to cut through superconducting fluxtubes. This mechanism is interesting because the oscillation modes need to exceed a critical amplitude in order for it to operate. Once it acts the effect is very strong (and nonlinear) leading to efficient damping. The upshot of this is that modes are unlikely to ever evolve far beyond the critical amplitude. We consider the effect of this new dissipation channel on the r-modes, that may be driven unstable by the emission of gravitational waves. Our estimates show that the fluxtube cutting leads to a saturation threshold for the instability that can be smaller than that of other proposed mechanisms. This suggests that the idea may be of direct astrophysical relevance

Mesoscopic pinning forces in neutron star crusts

The crust of a neutron star is thought to be comprised of a lattice of nuclei immersed in a sea of free electrons and neutrons. As the neutrons are superfluid their angular momentum is carried by an array of quantized vortices. These vortices can pin to the nuclear lattice and prevent the neutron superfluid from spinning down, allowing it to store angular momentum which can then be released catastrophically, giving rise to a pulsar glitch. A crucial ingredient for this model is the maximum pinning force that the lattice can exert on the vortices, as this allows us to estimate the angular momentum that can be exchanged during a glitch. In this paper we perform, for the first time, a detailed and quantitative calculation of the pinning force \emph{per unit length} acting on a vortex immersed in the crust and resulting from the mesoscopic vortex-lattice interaction. We consider realistic vortex tensions, allow for displacement of the nuclei and average over all possible orientation of the crystal with respect to the vortex. We find that, as expected, the mesoscopic pinning force becomes weaker for longer vortices and is generally much smaller than previous estimates, based on vortices aligned with the crystal. Nevertheless the forces we obtain still have maximum values of order $f_{\rm{pin}}\approx 10^{15}$ dyn/cm, which would still allow for enough angular momentum to be stored in the crust to explain large Vela glitches, if part of the star is decoupled during the event.Comment: 17 pages, 16 figures, 5 table