Slowly rotating two-fluid neutron star model

We study stationary axisymmetric configurations of a star model consisting of two barotropic fluids, which are uniformly rotating at two different rotation rates. Analytic approximate solutions in the limit of slow rotation are obtained with the classical method of Chandrasekhar, which consists of an expansion of the solution in terms of the rotation rate, and which is generalized to the case of two fluids in order to apply it to the present problem. This work has a direct application to neutron star models, in which the neutron superfluid can rotate at a different speed than the fluid of charged components. Two cases are considered, the case of two non-interacting fluids, and the case of an interaction of a special type, corresponding to the vortices of the neutron superfluid being completely pinned to the second fluid. The special case of the equation of state P proportional to rho(2) is solved explicitly as an illustration of the foregoing results.