We study non-radial oscillations of neutron stars with superfluid baryons, in a general relativistic framework, including finite temperature effects. Using a perturbative approach, we derive the equations describing stellar oscillations, which we sol
ve by numerical integration, employing different models of nucleon superfluidity, and determining frequencies and gravitational damping times of the quasi-normal modes. As expected by previous results, we find two classes of modes, associated to superfluid and non-superfluid degrees of freedom, respectively. We study the temperature dependence of the modes, finding that at specific values of the temperature, the frequencies of the two classes of quasi-normal modes show avoided crossings, and their damping times become comparable. We also show that, when the temperature is not close to the avoided crossings, the frequencies of the modes can be accurately computed by neglecting the coupling between normal and superfluid degrees of freedom. Our results have potential implications on the gravitational wave emission from neutron stars.
We demonstrate a possibility of existence of a peculiar temperature-dependent composition $g$-modes in superfluid neutron stars. We calculate the Brunt-V$ddot{rm a}$is$ddot{rm a}$l$ddot{rm a}$ frequency for these modes, as well as their eigenfrequenc
ies. The latter turn out to be rather large, up to $sim 500$ Hz for a chosen model of a neutron star. This result indicates, in particular, that use of the barotropic equation of state may be not a good approximation for calculation of inertial modes even in most rapidly rotating superfluid neutron stars.
We suggest a specific new class of low-frequency g-modes in superfluid neutron stars. We determine the Brunt-Vaisala frequency for these modes and demonstrate that they can be unstable with respect to convection. The criterion for the instability ons
et (analogue of the well known Schwarzschild criterion) is derived. It is very sensitive to equation of state and a model of nucleon superfluidity. In particular, convection may occur for both positive and negative temperature gradients. Our results have interesting implications for neutron star cooling and seismology.
We analyze damping of oscillations of general relativistic superfluid neutron stars. To this aim we extend the method of decoupling of superfluid and normal oscillation modes first suggested in [Gusakov & Kantor PRD 83, 081304(R) (2011)]. All calcula
tions are made self-consistently within the finite temperature superfluid hydrodynamics. The general analytic formulas are derived for damping times due to the shear and bulk viscosities. These formulas describe both normal and superfluid neutron stars and are valid for oscillation modes of arbitrary multipolarity. We show that: (i) use of the ordinary one-fluid hydrodynamics is a good approximation, for most of the stellar temperatures, if one is interested in calculation of the damping times of normal f-modes; (ii) for radial and p-modes such an approximation is poor; (iii) the temperature dependence of damping times undergoes a set of rapid changes associated with resonance coupling of neighboring oscillation modes. The latter effect can substantially accelerate viscous damping of normal modes in certain stages of neutron-star thermal evolution.
We show that suppression of the baryon energy gaps, caused by the relative motion of superfluid and normal liquid components, can substantially influence dynamical properties and evolution of neutron stars. This effect has been previously ignored in the neutron-star literature.
We show that equations governing pulsations of superfluid neutron stars can be splitted into two sets of weakly coupled equations, one describing the superfluid modes and another one -- the normal modes. The coupling parameter s is small, |s| ~ 0.01-
0.05, for realistic equations of state. Already an approximation s=0 is sufficient to calculate the pulsation spectrum within the accuracy of a few percents. Our results indicate, in particular, that emission of gravitational waves from superfluid pulsation modes is suppressed in comparison to that from normal modes. The proposed approach allows to drastically simplify modeling of pulsations of superfluid neutron stars.
