We formulate hydrodynamic equations for nonsuperfluid multicomponent magnetized charged relativistic mixtures, taking into account chemical reactions as well as viscosity, diffusion, thermodiffusion, and thermal conductivity effects. The resulting equations have a rather simple form and can be readily applied, e.g., for studying magnetothermal evolution of neutron stars. We also establish a link between our formalism and the results known in the literature, and express the phenomenological diffusion coefficients through momentum transfer rates which are calculated from microscopic theory.

We study equation of state (EOS) of an accreting neutron star crust. Usually, such EOS is obtained assuming (implicitly) that the free (unbound) neutrons and nuclei in the inner crust move together. We argue, that this assumption violates the condition $\mu_n^\infty={\rm const}$, required for hydrostatic (and diffusion) equilibrium of unbound neutrons ($\mu^\infty_n$ is the redshifted neutron chemical potential). We construct a new EOS respecting this condition, working in the compressible liquid-drop approximation. We demonstrate that it is close to the catalyzed EOS in most part of the inner crust, being very different from EOSs of accreted crust discussed in the literature. In particular, the pressure at the outer-inner crust interface does not coincide with the neutron drip pressure, usually calculated in the literature, and is determined by hydrostatic (and diffusion) equilibrium conditions within the star. We also find an instability at the bottom of fully accreted crust that transforms nuclei into homogeneous nuclear matter. It guarantees that the structure of fully accreted crust remains self-similar during accretion.

It is well-known that r-mode oscillations of rotating neutron stars may be unstable with respect to the gravitational wave emission. It is highly unlikely to observe a neutron star with the parameters within the instability window, a domain where this instability is not suppressed. But if one adopts the `minimal' (nucleonic) composition of the stellar interior, a lot of observed stars appear to be within the r-mode instability window. One of the possible solutions to this problem is to account for hyperons in the neutron star core. The presence of hyperons allows for a set of powerful (lepton-free) non-equilibrium weak processes, which increase the bulk viscosity, and thus suppress the r-mode instability. Existing calculations of the instability windows for hyperon NSs generally use reaction rates calculated for the $\Sigma^-\Lambda$ hyperonic composition via the contact $W$ boson exchange interaction. In contrast, here we employ hyperonic equations of state where the $\Lambda$ and $\Xi^-$ are the first hyperons to appear (the $\Sigma^-$'s, if they are present, appear at much larger densities), and consider the meson exchange channel, which is more effective for the lepton-free weak processes. We calculate the bulk viscosity for the non-paired $npe\mu\Lambda\Xi^-$ matter using the meson exchange weak interaction. A number of viscosity-generating non-equilibrium processes is considered (some of them for the first time in the neutron-star context). The calculated reaction rates and bulk viscosity are approximated by simple analytic formulas, easy-to-use in applications. Applying our results to calculation of the instability window, we argue that accounting for hyperons may be a viable solution to the r-mode problem.

Force on proton vortices in superfluid and superconducting matter of neutron stars is calculated at vanishing stellar temperature. Both longitudinal (dissipative) and transverse (Lorentz-type) components of the force are derived in a coherent way and compared in detail with the corresponding expressions available in the literature. This allows us to resolve a controversy about the form of the Lorentz-type force component acting on proton vortices. The calculated force is a key ingredient in magnetohydrodynamics of superconducting neutron stars and is important for modeling the evolution of stellar magnetic field.

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 study the bulk viscosity in two (anti-spaghetti and Swiss cheese) phases of non- spherical nuclei in the mantle of an oscillating neutron star near the boundary with the stellar core. The bulk viscosity is produced by non-equilibrium Urca neutrino emis- sion processes. In the mantle, the direct Urca process may be open (Gusakov et al., 2004) if neutrons and protons move in a periodic potential created by a lattice of non- spherical nuclei (which allows the nucleons to have large quasi-momenta and satisfy direct Urca momentum-conservation). This bulk viscosity can dominate over that due to the modified Urca process in the outer stellar core and over the shear viscosity. The bulk viscosity depends strongly on temperature, oscillation frequency and nucleon superfluidity. The enhanced bulk viscosity in the mantle can control propagation and damping of neutron star oscillations.

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.

In [Gusakov et al.\ PRD, 96, 103012, (2017)], we proposed a self-consistent method to study the quasistationary evolution of the magnetic field in neutron-star cores. Here we apply it to calculate the instantaneous particle velocities and other parameters of interest, which are fixed by specifying the magnetic field configuration. Interestingly, we found that the magnetic field can lead to generation of a macroscopic fluid motion with the velocity, significantly exceeding the diffusion particle velocities. This result calls into question the standard view on the magnetic field evolution in neutron stars and suggests a new, shorter timescale for such evolution.

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].

Buoyancy of proton vortices is considered as one of the important mechanisms of magnetic field expulsion from the superconducting interiors of neutron stars. Here we show that the generally accepted expression for the buoyancy force is not correct and should be modified. The correct expression is derived for both neutron and proton vortices. It is argued that this force is already contained in the coarse-grained hydrodynamics of Bekarevich & Khalatnikov and its various multifluid extensions, but is absent in the hydrodynamics of Hall. Some potentially interesting buoyancy-related effects are briefly discussed.

The relativistic dynamic equations are derived for a superfluid-superconducting mixture coupled to the electromagnetic field. For definiteness, and bearing in mind possible applications of our results to neutron stars, it is assumed that the mixture is composed of superfluid neutrons, superconducting protons, and normal electrons. Proton superconductivity of both I and II types is analysed, and possible presence of neutron and proton vortices (or magnetic domains in the case of type-I proton superconductivity) is allowed for. The derived equations neglect all dissipative effects except for the mutual friction dissipation and are valid for arbitrary temperatures (i.e. they do not imply that all nucleons are paired), which is especially important for magnetar conditions. It is demonstrated that these general equations can be substantially simplified for typical neutron stars, for which a kind of magnetohydrodynamic approximation is justified. Our results are compared to the nonrelativistic formulations existing in the literature and a number of discrepancies are found. In particular, it is shown that, generally, the electric displacement ${\pmb D}$ does not coincide with the electric field ${\pmb E}$, contrary to what is stated in the previous works. The relativistic framework developed here is easily extendable to account for more sophisticated microphysics models and it provides the necessary basis for realistic modelling of neutron stars.

The relativistic analogue of the Hall-Vinen-Bekarevich-Khalatnikov (HVBK) hydrodynamics is derived making use of the phenomenological method similar to that used by Bekarevich and Khalatnikov [1] in their derivation of HVBK-hydrodynamics. The resulting equations describe a finite-temperature superfluid liquid with the distributed vorticity. The main dissipative effects, including mutual friction, are taken into account. The proposed hydrodynamics is needed for reliable modeling of the dynamical properties of superfluid neutron stars.

We analyse the oscillations of general relativistic superfluid hyperon stars, following the approach suggested by Gusakov & Kantor and Gusakov et al. and generalizing it to the nucleon-hyperon matter. We show that the equations governing the oscillations can be split into two weakly coupled systems with the coupling parameters $s_{\rm e}$, $s_{\rm \mu}$, and $s_{\rm str}$. The approximation $s_{\rm e} = s_{\rm \mu} = s_{\rm str} = 0$ (decoupling approximation) allows one to drastically simplify the calculations of stellar oscillation spectra. An efficiency of the presented scheme is illustrated by the calculation of sound speeds in the nucleon-hyperon matter composed of neutrons (n), protons (p), electrons (e), muons ($\mu$), as well as $\rm \Lambda$, ${\rm \Xi}^-$, and ${\rm \Xi}^0$-hyperons. However, the gravity oscillation modes (g-modes) cannot be treated within this approach, and we discuss them separately. For the first time we study the composition g-modes in superfluid hyperon stars with the $\rm npe\mu\Lambda$ core and show that there are two types of g-modes (`muonic' and `$\Lambda$--hyperonic') in such stars. We also calculate the g-mode spectrum and find out that the eigenfrequencies $\nu$ of the superfluid g-modes can be exceptionally large (up to $\nu \approx 742~{\rm Hz}$ for a considered stellar model).

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.

For the first time nonradial oscillations of superfluid nonrotating stars are self-consistently studied at finite stellar temperatures. We apply a realistic equation of state and realistic density dependent model of critical temperature of neutron and proton superfluidity. In particular, we discuss three-layer configurations of a star with no neutron superfluidity at the centre and in the outer region of the core but with superfluid intermediate region. We show, that oscillation spectra contain a set of modes whose frequencies can be very sensitive to temperature variations. Fast temporal evolution of the pulsation spectrum in the course of neutron star cooling is also analysed.

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 study transport properties of a strongly interacting superfluid mixture of two Fermi-liquids. A typical example of such matter is the neutron-proton liquid in the cores of neutron stars. To describe the mixture, we employ the Landau theory of Fermi-liquids, generalized to allow for the effects of superfluidity. We formulate the kinetic equation and analyze linear response of the system to vector (e.g., electromagnetic) perturbation. In particular, we calculate the transverse and longitudinal polarization functions for both liquid components. We demonstrate, that they can be expressed through the Landau parameters of the mixture and polarization functions of noninteracting matter (when the Landau quasiparticle interaction is neglected). Our results can be used, e.g., for studies of the kinetic coefficients and low-frequency long-wavelength collective modes in superfluid Fermi-mixtures.

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.

The hydrodynamics, describing dynamical effects in superfluid neutron stars, essentially differs from the standard one-fluid hydrodynamics. In particular, we have four bulk viscosity coefficients in the theory instead of one. In this paper we calculate these coefficients, for the first time, assuming they are due to non-equilibrium beta-processes (such as modified or direct Urca process). The results of our analysis are used to estimate characteristic damping times of sound waves in superfluid neutron stars. It is demonstrated that all four bulk viscosity coefficients lead to comparable dissipation of sound waves and should be considered on the same footing.

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.

We examine radial oscillations of superfluid neutron stars at finite internal temperatures. For this purpose we generalize the description of relativistic superfluid hydrodynamics to the case of superfluid mixtures. We show that in a neutron star at hydrostatic and beta-equilibrium the red-shifted temperature gradient is smoothed out by neutron superfluidity (but not by proton superfluidity). We calculate radial oscillation modes of neutron stars assuming "frozen" nuclear composition in the pulsating matter. The resulting pulsation frequencies show a strong temperature dependence in the temperature range (0.1-1) T_cn, where T_cn is the critical temperature of neutron superfluidity. Combining our results with thermal evolution, we obtain a significant evolution of the pulsation spectrum, associated with highly efficient Cooper pairing neutrino emission, for 20 years after superfluidity onset.

We study the "minimal" cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and the enhanced cooling is produced by the neutrino emission due to Cooper pairing of neutrons. Extending our previous consideration (Gusakov et al. 2004a), we include the effects of accreted envelopes of light elements. We employ phenomenological density-dependent critical temperatures T_{cp}(\rho) and T_{cnt}(\rho) of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature T_{cns}(\rho) of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of Gusakov et al. (2004a) and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations.

We have derived a set of equations to describe the thermal evolution of a neutron star which undergoes small-amplitude radial pulsations. We have taken into account, in the frame of the General Theory of Relativity, the pulsation damping due to the bulk and shear viscosity and the accompanying heating of the star. The neutrino emission of a pulsating non-superfluid star and its heating due to the bulk viscosity are calculated assuming that both processes are determined by the non-equilibrium modified Urca process. Analytical and numerical solutions to the set of equations of the stellar evolution are obtained for linear and strongly non-linear deviations from beta-equilibrium. It is shown that a pulsating star may be heated to very high temperatures, while the pulsations damp very slowly with time (a power law damping for 100-1000 years), as long as the damping is determined by the bulk viscosity. The contribution of the shear viscosity to the damping becomes important in a rather cool star with a low pulsation energy.

The entrainment matrix (also termed the Andreev-Bashkin matrix or the mass-density matrix) for a neutron-proton mixture is derived at a finite temperature in a neutron star core. The calculation is performed in the frame of the Landau Fermi-liquid theory generalized to account for superfluidity of nucleons. It is shown, that the temperature dependence of the entrainment matrix is described by a universal function independent on an actual model of nucleon-nucleon interaction employed. The results are presented in the form convenient for their practical use. The entrainment matrix is important, e.g., in kinetics of superfluid nucleon mixtures or in studies of the dynamical evolution of neutron stars (in particular, in the studies of star pulsations and pulsar glitches).

We study the cooling of superfluid neutron stars whose cores consist of nucleon matter with the Akmal-Pandharipande-Ravenhall equation of state. This equation of state opens the powerful direct Urca process of neutrino emission in the interior of most massive neutron stars. Extending our previous studies (Gusakov et al. 2004a, Kaminker et al. 2005), we employ phenomenological density-dependent critical temperatures T_{cp}(\rho) of strong singlet-state proton pairing (with the maximum T_{cp}^{max} \sim 7e9 K in the outer stellar core) and T_{cnt}(\rho) of moderate triplet-state neutron pairing (with the maximum T_{cnt}^{max} \sim 6e8 K in the inner core). Choosing properly the position of T_{cnt}^{max} we can obtain a representative class of massive neutron stars whose cooling is intermediate between the cooling enhanced by the neutrino emission due to Cooper pairing of neutrons in the absence of the direct Urca process and the very fast cooling provided by the direct Urca process non-suppressed by superfluidity.

We investigate the combined effect of neutron and proton superfluidities on the cooling of neutron stars whose cores consist of nucleons and electrons. We consider singlet-state pairing of protons and triplet-state pairing of neutrons in the cores of neutron stars. The critical superfluid temperatures T_c are assumed to depend on the density of matter. We study two types of neutron pairing with different components of the total angular momentum of Cooper pairs along the quantization axis (|m_J| =0 or 2). Our calculations are compared with observations of thermal emission from isolated neutron stars. We show that the observations can be interpreted by using two classes of superfluidity models: (1) strong proton superfluidity with a maximum critical temperature in the stellar core T_c^{max} > 4 \times 10^9 K and weak neutron superfluidity of any type (T_c^{max} < 2 \times 10^8 K); (2) strong neutron superfluidity (pairing with |m_J|=0) and weak proton superfluidity. The two types of models reflect an approximate symmetry with respect to an interchange of the critical temperatures of neutron and proton pairing.

The impact of nuclear physics theories on cooling of isolated neutron stars is analyzed. Physical properties of neutron star matter important for cooling are reviewed such as composition, the equation of state, superfluidity of various baryon species, neutrino emission mechanisms. Theoretical results are compared with observations of thermal radiation from neutron stars. Current constraints on theoretical models of dense matter, derived from such a comparison, are formulated.

We show that the direct Urca process of neutrino emission is allowed in two possible phases of nonspherical nuclei (inverse cylinders and inverse spheres) in the mantle of a neutron star near the crust-core interface. The process is open because neutrons and protons move in a periodic potential created by inhomogeneous nuclear structures. In this way the nucleons acquire large quasimomenta needed to satisfy momentum-conservation in the neutrino reaction. The appropriate neutrino emissivity in a nonsuperfluid matter is about 2--3 orders of magnitude higher than the emissivity of the modified Urca process in the stellar core. The process may noticeably accelerate the cooling of low-mass neutron stars.

We simulate cooling of superfluid neutron stars with nucleon cores where direct Urca process is forbidden. We adopt density dependent critical temperatures $T_{cp}(\rho)$ and $T_{cn}(\rho)$ of singlet-state proton and triplet-state neutron pairing in a stellar core and consider a strong proton pairing (with maximum $T_{cp}^{max} \ga 5 \times 10^9$ K) and a moderate neutron pairing ($T_{cn}^{max} \sim 6 \times 10^8$ K). When the internal stellar temperature $T$ falls below $T_{cn}^{max}$, the neutrino luminosity $L_{CP}$ due to Cooper pairing of neutrons behaves $\propto T^8$, just as that produced by modified Urca process (in a non-superfluid star) but is higher by about two orders of magnitude. In this case the Cooper-pairing neutrino emission acts like an enhanced cooling agent. By tuning the density dependence $T_{cn}(\rho)$ we can explain observations of cooling isolated neutron stars in the scenario in which direct Urca process or similar process in kaon/pion condensed or quark matter are absent.

We consider cooling of neutron stars (NSs) with superfluid cores composed of neutrons, protons, and electrons (assuming singlet-state pairing of protons, and triplet-state pairing of neutrons). We mainly focus on (nonstandard) triplet-state pairing of neutrons with the $|m_J| = 2$ projection of the total angular momentum of Cooper pairs onto quantization axis. The specific feature of this pairing is that it leads to a power-law (nonexponential) reduction of the emissivity of the main neutrino processes by neutron superfluidity. For a wide range of neutron critical temperatures $T_{cn}$, the cooling of NSs with the $|m_J| = 2$ superfluidity is either the same as the cooling with the $m_J = 0$ superfluidity, considered in the majority of papers, or much faster. The cooling of NSs with density dependent critical temperatures $T_{cn}(\rho)$ and $T_{cp}(\rho)$ can be imitated by the cooling of the NSs with some effective critical temperatures $T_{cn}$ and $T_{cp}$ constant over NS cores. The hypothesis of strong neutron superfluidity with $|m_J| = 2$ is inconsistent with current observations of thermal emission from NSs, but the hypothesis of weak neutron superfluidity of any type does not contradict to observations.

We calculate and provide analytic fits of the factors which describe the reduction of the neutrino emissivity of modified Urca and nucleon-nucleon bremsstrahlung processes by superfluidity of neutrons and protons in neutron-star cores. We consider $^1$S$_0$ pairing of protons and either $^1$S$_0$ or $^3$P$_2$ pairing of neutrons. We analyze two types of $^3$P$_2$ pairing: the familiar pairing with zero projection of the total angular momentum of neutron pairs onto quantization axis, $m_J=0$; and the pairing with $|m_J|=2$ which leads to the gap with nodes at the neutron Fermi surface. Combining the new data with those available in the literature we fully describe neutrino emission by nucleons from neutron star cores to be used in simulations of cooling of superfluid neutron stars.