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

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.

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