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

@misc{imam2026massive,
  title = {{Massive cold hybrid stars in a modified Polyakov-Nambu-Jona-Lasinio model}},
  author = {Sk Md Adil Imam and Pedro Costa and Mariana Dutra and Odilon Lourenço and Renan Pereira and Constança Providência},
  howpublished = {arXiv},
  year = {2026},
  month = sep,
  doi = {10.48550/arxiv.2608.12653},
  eprint = {2608.12653},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2608.12653},
  note = {CC BY 4.0},
}

TY  - UNPB
TI  - Massive cold hybrid stars in a modified Polyakov-Nambu-Jona-Lasinio model
AU  - Sk Md Adil Imam
AU  - Pedro Costa
AU  - Mariana Dutra
AU  - Odilon Lourenço
AU  - Renan Pereira
AU  - Constança Providência
PY  - 2026
DA  - 2026/09/07
PB  - arXiv
DO  - 10.48550/arxiv.2608.12653
UR  - https://arxiv.org/abs/2608.12653
AB  - 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 quark matter and subsequently to deconfined quark matter, or a single transition directly from hadronic to deconfined quark matter or from hadronic to confined quark matter. Stable massive cold hybrid stars with only confined 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 confined and/or a deconfined core can reach maximum masses above 2M_⊙, 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 shows variations at finite baryon densities, with c_s^2(μ) departing from the asymptotic conformal value c_s^2 = 1/3 in confined core stars. These variations serve as a diagnostic of the equation-of-state stiffness while remaining fully consistent with causality and thermodynamic stability. This work establishes the qualitative role of each model parameter in shaping hybrid-star structure.
ER  -