We calculate the spectrum of inertial oscillation modes in a slowly rotating superfluid neutron star, including, for the first time, both the effects of finite temperatures and entrainment between superfluid neutrons and protons. We work in the Newtonian limit and assume minimal core composition (neutrons, protons and electrons). We also developed an approximate method that allows one to calculate the superfluid r-mode analytically. Finally, we derive and analyze dispersion relations for inertial modes in the superfluid NS matter in the short wavelength limit.

We propose a general method to self-consistently study the quasistationary evolution of the magnetic field in the cores of neutron stars. The traditional approach to this problem is critically revised. Our results are illustrated by calculation of the typical timescales for the magnetic field dissipation as functions of temperature and the magnetic field strength.

We address the problem of magnetic field dissipation in the neutron star cores, focusing on the role of neutron superfluidity. Contrary to the results in the literature, we show that in the finite-temperature superfluid matter composed of neutrons, protons, and electrons, magnetic field dissipates exclusively due to Ohmic losses and non-equilibrium beta-processes, and only an admixture of muons restores (to some extent) the role of particle relative motion for the field dissipation. The reason for this discrepancy is discussed.

We calculate the finite-temperature r-mode spectrum of a slowly rotating superfluid Newtonian neutron star neglecting the entrainment between neutron and proton liquid components (i.e., neglecting the off-diagonal element of the entrainment matrix). We show that for `minimal' NS core composition (neutrons, protons, and electrons) only two $m=2$ r-modes exist --- normal mode, which is similar to ordinary r-mode in a nonsuperfluid star, and a superfluid temperature-dependent mode. Accounting for muons in the core dramatically modifies the oscillation spectrum, resulting in an infinite set of superfluid r-modes, whose frequencies vary with temperature. We demonstrate that the normal r-mode can exhibit avoided crossings with superfluid modes at certain `resonance' temperatures, where it dissipates strongly, which leads to substantial suppression of the r-mode instability near these temperatures. The corresponding instability windows are calculated and discussed.

To put new constraints on the r-mode instability window, we analyse the formation of millisecond pulsars (MSPs) within the recycling scenario, making use of three sets of observations: (a) X-ray observations of neutron stars (NSs) in low-mass X-ray binaries; (b) timing of millisecond pulsars; and (c) X-ray and UV observations of MSPs. As shown in previous works, r-mode dissipation by shear viscosity is not sufficient to explain observational set (a), and enhanced r-mode dissipation at the red-shifted internal temperatures $T^\infty\sim 10^8$ K is required to stabilize the observed NSs. Here, we argue that models with enhanced bulk viscosity can hardly lead to a self-consistent explanation of observational set (a) due to strong neutrino emission, which is typical for these models (unrealistically powerful energy source is required to keep NSs at the observed temperatures). We also demonstrate that the observational set (b), combined with the theory of internal heating and NS cooling, provides evidence of enhanced r-mode dissipation at low temperatures, $T^\infty\sim 2\times 10^7$ K. Observational set (c) allows us to set an upper limit on the internal temperatures of MSPs, $T^\infty<2\times 10^7$ K (assuming a canonical NS with the accreted crust). Recycling scenario can produce MSPs at these temperatures only if r-mode instability is suppressed in the whole MSP spin frequency range ($\nu\lesssim 750$ Hz) at temperatures $2\times 10^7\lesssim T^\infty\lesssim 3 \times 10^7$ K, providing thus a new constraint on the r-mode instability window. These observational constraints are analysed in more details in application to the resonance uplift scenario of Gusakov et al. [Phys. Rev. Lett., 112 (2014), 151101].

In the recent papers by Gusakov, Chugunov, and Kantor (2014) a new scenario describing evolution of rapidly rotating neutron stars in low-mass X-ray binaries was proposed. The scenario accounts for a resonant interaction of normal r modes with superfluid inertial modes at some specific internal stellar temperatures ("resonance temperatures"). This interaction results in an enhanced damping of r mode and appearance of the "stability peaks" in the temperature -- spin frequency plane, which split the r-mode instability window in the vicinity of the resonance temperatures. The scenario suggests that the hot and rapidly rotating NSs spend most of their life climbing up these peaks and, in particular, are observed there at the moment. We analyze in detail possible observational signatures of this suggestion. In particular, we show that these objects may exhibit `anti-glitches' -- sudden frequency jumps on a time scale of hours-months.

The spin-down of a neutron star, e.g. due to magneto-dipole losses, results in compression of the stellar matter and induces nuclear reactions at phase transitions between different nuclear species in the crust. We show that this mechanism is effective in heating recycled pulsars, in which the previous accretion process has already been compressing the crust, so it is not in nuclear equilibrium. We calculate the corresponding emissivity and confront it with available observations, showing that it might account for the likely thermal ultraviolet emission of PSR J0437-4715.

Recent observation of a sudden spin down, occurring on a timescale not exceeding two weeks, of the magnetar 1E2259+586 (see Archibald et al. 2013, where this event was dubbed an `anti-glitch') has not still received any interpretation in terms of the standard scenario of pulsar glitches proposed by Anderson and Itoh (1975). Motivated by this observation, here we present a toy model that allows, under certain conditions, for anti-glitches in neutron stars within the standard approach.

A new class of neutron stars (NSs) -- hot rapidly rotating non-accreting NSs, which we propose to call HOFNARs (HOt and Fast Non-Accreting Rotators) or "hot widows" (in analogy with "black widow" pulsars) -- is suggested. We argue that such stars should originate from the low-mass X-ray binaries (LMXBs) provided that they were unstable with respect to excitation of $r$-modes at the end of accretion epoch (when their low-mass companions ceased to fill the Roche lobe). High temperature of "hot widows"/HOFNARs is maintained by $r$-mode dissipation rather than by accretion. We analyse observational properties of "hot widows"/HOFNARs and demonstrate that these objects form a specific separate class of neutron stars. In particular, some of the most stable X-ray sources among the candidates to quiescent LMXB systems (qLMXBs), can, in fact, belong to that new class. We formulate observational criteria which allow to distinguish "hot widows"/HOFNARs from qLMXB systems, and argue that available observations of X-ray sources 47 Tuc X5 and X7 satisfy (or, at least, do not contradict) these criteria. In addition, we discuss indirect evidences in favor of "hot widows"/HOFNARs existence, following from the analysis of observations and predictions of population synthesis theories. If that new class of NSs does exist, it would prove the possibility to emit gravitational waves by mass-current multipole. Various applications of our results, such as prospects for constraining superdense matter properties with hot widows"/HOFNARs, are analyzed.

In a previous paper [M. E. Gusakov, A. I. Chugunov, and E. M. Kantor, Phys. Rev. Lett. 112, 151101 (2014)], we introduced a new scenario that explains the existence of rapidly rotating warm neutron stars (NSs) observed in low-mass X-ray binaries (LMXBs). Here it is described in more detail. The scenario takes into account the interaction between superfluid inertial modes and the normal (quadrupole) $m=2$ $r$-mode, which can be driven unstable by Chandrasekhar-Friedman-Schutz (CFS) mechanism. This interaction can only occur at some fixed "resonance" stellar temperatures; it leads to formation of the "stability peaks" which stabilize a star in the vicinity of these temperatures. We demonstrate that a NS in LMXB spends a substantial fraction of time on the stability peak, that is, in the region of stellar temperatures and spin frequencies, that has been previously thought to be CFS unstable with respect to excitation of $r$-modes. We also find that the spin frequencies of NSs are limited by the CFS instability of normal (octupole) $m=3$ $r$-mode rather than by $m=2$ $r$-mode. This result agrees with the predicted value of the cutoff spin frequency $\sim 730$ Hz in the spin distribution of accreting millisecond X-ray pulsars. In addition, we analyze evolution of a NS after the end of the accretion phase and demonstrate that millisecond pulsars can be born in LMXBs within our scenario. Besides millisecond pulsars, our scenario also predicts a new class of LMXB descendants - hot and rapidly rotating nonaccreting NSs ("hot widows"/HOFNARs). Further comparison of the proposed theory with observations of rotating NSs can impose new important constraints on the properties of superdense matter.

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 solve 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 eigenfrequencies. 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 consider an instability of rapidly rotating neutron stars in low-mass X-ray binaries (LMXBs) with respect to excitation of r-modes (which are analogous to Earth's Rossby waves controlled by the Coriolis force). We argue that finite temperature effects in the superfluid core of a neutron star lead to a resonance coupling and enhanced damping (and hence stability) of oscillation modes at certain stellar temperatures. Using a simple phenomenological model we demonstrate that neutron stars with high spin frequency may spend a substantial amount of time at these `resonance' temperatures. This finding allows us to explain puzzling observations of hot rapidly rotating neutron stars in LMXBs and to predict a new class of hot, non-accreting, rapidly rotating neutron stars, some of which may have already been observed and tentatively identified as quiescent LMXB (qLMXB) candidates. We also impose a new theoretical limit on the neutron star spin frequency, explaining the cut-off spin frequency ~730 Hz, following from the statistical analysis of accreting millisecond X-ray pulsars. Besides explaining the observations, our model provides a new tool to constrain superdense matter properties comparing measured and theoretically predicted resonance temperatures.

Observations of massive ($M \approx 2.0~M_\odot$) neutron stars (NSs), PSRs J1614-2230 and J0348+0432, rule out most of the models of nucleon-hyperon matter employed in NS simulations. Here we construct three possible models of nucleon-hyperon matter consistent with the existence of $2~M_\odot$ pulsars as well as with semi-empirical nuclear matter parameters at saturation, and semi-empirical hypernuclear data. Our aim is to calculate for these models all the parameters necessary for modelling dynamics of hyperon stars (such as equation of state, adiabatic indices, thermodynamic derivatives, relativistic entrainment matrix, etc.), making them available for a potential user. To this aim a general non-linear hadronic Lagrangian involving $\sigma\omega\rho\phi\sigma^\ast$ meson fields, as well as quartic terms in vector-meson fields, is considered. A universal scheme for calculation of the $\ell=0,1$ Landau Fermi-liquid parameters and relativistic entrainment matrix is formulated in the mean-field approximation. Use of this scheme allow us to obtain numerical tables with the equation of state, Landau quasiparticle effective masses, adiabatic indices, the $\ell=0,1$ Landau Fermi-liquid parameters, and the relativistic entrainment matrix for the selected models of nucleon-hyperon matter. These data are available on-line and suitable for numerical implementation in computer codes modelling various dynamical processes in NSs, in particular, oscillations of superfluid NSs and their cooling.

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 onset (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 calculations 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.

We study the effects of finite stellar temperatures on the oscillations of superfluid neutron stars. The importance of these effects is illustrated with a simple example of a radially pulsating general relativistic star. Two main effects are taken into account: (i) temperature dependence of the entrainment matrix and (ii) the variation of the size of superfluid region with temperature. Four models are considered, which include either one or both of these two effects. Pulsation spectra are calculated for these models, and asymptotes for eigenfrequencies at temperatures close to critical temperature of neutron superfluidity, are derived. It is demonstrated that models that allow for the temperature effect (ii) but disregard the effect (i), yield unrealistic results. Eigenfunctions for the normal- and superfluid-type pulsations are analyzed. It is shown that superfluid pulsation modes practically do not appear at the neutron-star surface and, therefore, can hardly be observed by measuring the modulation of the electromagnetic radiation from the star. The e-folding times for damping of pulsations due to the shear viscosity and nonequilibrium modified Urca processes are calculated and their asymptotes at temperatures close to the neutron critical temperature, are obtained. It is demonstrated that superfluid pulsation modes are damped by 1--3 orders of magnitude faster than normal modes.

We calculate the important quantity of superfluid hydrodynamics, the relativistic entrainment matrix for a nucleon-hyperon mixture at arbitrary temperature. In the nonrelativistic limit this matrix is also termed the Andreev-Bashkin or mass-density matrix. Our results can be useful for modeling the pulsations of massive neutron stars with superfluid nucleon-hyperon cores and for studies of the kinetic properties of superfluid baryon matter.

We calculate the relativistic entrainment matrix Y_ik at zero temperature for nucleon-hyperon mixture composed of neutrons, protons, Lambda and Sigma^- hyperons, as well as of electrons and muons. This matrix is analogous to the entrainment matrix (also termed mass-density matrix or Andreev-Bashkin matrix) of non-relativistic theory. It is an important ingredient for modelling the pulsations of massive neutron stars with superfluid nucleon-hyperon cores. The calculation is done in the frame of the relativistic Landau Fermi-liquid theory generalized to the case of superfluid mixtures; the matrix Y_ik is expressed through the Landau parameters of nucleon-hyperon matter. The results are illustrated with a particular example of the sigma-omega-rho mean-field model with scalar self-interactions. Using this model we calculate the matrix Y_ik and the Landau parameters. We also analyze stability of the ground state of nucleon-hyperon matter with respect to small perturbations.

We consider sound waves in superfluid nucleon-hyperon matter of massive neutron-star cores. We calculate and analyze the speeds of sound modes and their damping times due to the shear viscosity and non-equilibrium weak processes of particle transformations. For that, we employ the dissipative relativistic hydrodynamics of a superfluid nucleon-hyperon mixture, formulated recently [M.E. Gusakov and E.M. Kantor, Phys. Rev. D78, 083006 (2008)]. We demonstrate that the damping times of sound modes calculated using this hydrodynamics and the ordinary (nonsuperfluid) one, can differ from each other by several orders of magnitude.

We calculated bulk viscosity due to non-equilibrium weak processes in superfluid nucleon-hyperon matter of neutron stars. For that, the dissipative relativistic hydrodynamics, formulated in paper [1] for superfluid mixtures, was extended to the case when both nucleons and hyperons are superfluid. It was demonstrated that in the most general case (when neutrons, protons, Lambda, and Sigma^{-} hyperons are superfluid), non-equilibrium weak processes generate sixteen bulk viscosity coefficients, with only three of them being independent. In addition, we corrected an inaccuracy in a widely used formula for the bulk viscosity of non-superfluid nucleon-hyperon matter.

One of the effective mechanisms of neutrino energy losses in red giants, presupernovae and in the cores of white dwarfs is the emission of neutrino-antineutrino pairs in the process of plasmon decay. In this paper, we numerically calculate the emissivity due to plasmon decay in a wide range of temperatures (10^7-10^11) K and densities (200-10^14) g cm^-3. Numerical results are approximated by convenient analytical expressions. We also calculate and approximate by analytical expressions the neutrino luminosity of white dwarfs due to plasmon decay, as a function of their mass and internal temperature. This neutrino luminosity depends on the chemical composition of white dwarfs only through the parameter mu_e (the net number of baryons per electron) and is the dominant neutrino luminosity in all white dwarfs at the neutrino cooling stage.

The electrostatic potential that keeps approximate charge neutrality in neutron star matter is self-consistently introduced into the formalism for rotochemical heating presented in a previous paper by Fernandez and Reisenegger. Although the new formalism is more rigorous, we show that its observable consequences are indistinguishable from those of the previous one, leaving the conclusions of the previous paper unchanged.

Comments on the article "Pulsar dynamics: magnetic dipole model revisited".