In this paper, we discuss the impact of rotation on the particle composition of rotating neutron stars (pulsars). Particular emphasis is put on the formation of quark matter during stellar spin-down, driven by continuous gravitational compression. Our study is based on modern models for the nuclear equation of state whose parameters are tightly constrained by nuclear data, neutron star masses, and the latest estimates of neutron star radii.

In this paper, we use a three flavor non-local Nambu--Jona-Lasinio (NJL) model, an~improved effective model of Quantum Chromodynamics (QCD) at low energies, to investigate the existence of deconfined quarks in the cores of neutron stars. Particular emphasis is put on the possible existence of quark matter in the cores of rotating neutron stars (pulsars). In contrast to non-rotating neutron stars, whose particle compositions do not change with time (are frozen in), the type and structure of the matter in the cores of rotating neutron stars depends on the spin frequencies of these stars, which opens up a possible new window on the nature of matter deep in the cores of neutron stars. Our study shows that, depending on mass and rotational frequency, up to around 8% of the mass of a massive neutron star may be in the mixed quark-hadron phase, if the phase transition is treated as a Gibbs transition. We also find that the gravitational mass at which quark deconfinement occurs in rotating neutron stars varies quadratically with spin frequency, which can be fitted by a simple formula.

Non-rotating neutron stars are generally treated in theoretical studies as perfect spheres. Such a treatment, however, may not be correct if strong magnetic fields are present (such as for magnetars) and/or the pressure of the matter in the cores of neutron stars is non-isotropic (e.g., color superconducting). In this paper, we investigate the structure of non-spherical neutron stars in the framework of general relativity. Using a parameterized metric to model non-spherical mass distributions, we first derive a stellar structure equation for deformed neutron stars. Numerical investigations of this model equation show that the gravitational masses of deformed neutron stars depend rather strongly on the degree and type (oblate or prolate) of stellar deformation. In particular, we find that the mass of a neutron star increases with increasing oblateness but decreases with increasing prolateness. If this feature carries over to a full two-dimensional treatment of deformed neutron stars, this opens up the possibility that, depending on the type of stellar deformation, there may exist multiple maximum-mass neutron stars for one and for the same model for the nuclear equation of state.