Massive cold hybrid stars in a modified Polyakov-Nambu-Jona-Lasinio model

We propose a modified Polyakov-loop Nambu--Jona-Lasinio (mPNJL) model in which the Polyakov potential is given by an explicit dependence on the quark chemical potential, allowing it to remain finite at zero temperature and thus to describe the confinement-deconfinement transition in cold dense matter. Combining this modified quark sector with hadronic equations of state via a Maxwell construction, we find that, depending on the model parameters, the equation of state can exhibit either two phase transitions, from hadronic matter to confined (quarkyonic) quark matter and subsequently to deconfined quark matter, or a single transition directly from hadronic to deconfined quark matter or from hadronic to quarkyonic quark matter. Stable massive cold hybrid stars with only quarkyonic and/or deconfined quark phase are obtained. We systematically examine how the parameters of the modified Polyakov potential and the quark vector interactions control the location of these transitions, and find that repulsive vector interactions are essential to obtain a stable quark core. Hybrid stars with quarkyonic and/or a deconfined core can reach maximum masses above $2M_\odot$, provided a sufficiently stiff hadronic equation of state is used at low density. In the core of the maximum-mass configurations, the speed of sound exceeds the conformal limit, $c_s^2 = 1/3$, for the quarkyonic core stars. This work establishes the qualitative role of each model parameter in shaping hybrid-star structure.