Rapidly rotating superfluid neutron stars in Newtonian dynamics

We develop a formulation for constructing and examining rapidly rotating Newtonian neutron star models that contain two superfluids, taking account of the effect of the rotation velocity difference between two superfluids. We assume neutron stars to be composed of the superfluid neutrons and a mixture of the superfluid protons and the normal fluid electrons. To describe Newtonian dynamics of the two superfluids, the Newtonian version of the so-called two-fluid formalism is employed. The effect of the rotation velocity difference on the structure of equilibrium state is treated as a small perturbation to rapidly rotating superfluid stars whose angular velocities of two superfluids are assumed to be exactly the same. We derive basic equations for the perturbed structures of rapidly rotating superfluid stars due to the rotation velocity difference between two superfluids. Assuming the superfluids to obey a simple analytical equation of state proposed by Prix, Comer and Andersson, we obtain numerical solutions for the perturbations and find that the density distributions of the superfluids are strongly dependent on the parameter σ, which appears in the analytical equation of state and characterizes the so-called symmetry energy. It is also found that if the analytical equation of state of Prix et al. is assumed, the perturbations can be represented in terms of the universal functions that are independent of the parameters of the equation of state.