We present a Bayesian analysis for a new class of realistic models of two-phase equations of state (EoS) for hybrid stars and demonstrate that the observation of a pair of high-mass twin stars would have a sufficient discriminating power to favor hybrid EoS with a strong first order phase transition over alternative EoS. Such a measurement would provide evidence for the existence of a critical endpoint in the QCD phase diagram.

We point out that the very existence of a "horizontal branch" in the mass-radius characteristics for neutron stars indicates a strong first-order phase transition and thus supports the existence of a critical endpoint (CEP) of first order phase transitions in the QCD phase diagram. This branch would sample a sequence of hybrid stars with quark matter core, leading to the endpoint of stable compact star configurations with the highest possible baryon densities. Since we know of the existence of compact stars with $2M_\odot$, this hypothetical branch has to lie in the vicinity of this mass value, if it exists. We report here a correlation between the maximal radius of the horizontal branch and the pressure at the onset of hadron-to-quark matter phase transition which is likely to be a universal quantity of utmost relevance to the upcoming experiments with heavy-ion collisions at NICA and FAIR.

We explore systematically a new class of two-phase equations of state (EoS) for hybrid stars that is characterized by three main features : (1) stiffening of the nuclear EoS at supersaturation densities due to quark exchange effects (Pauli blocking) between hadrons, modelled by an excluded volume correction, (2) stiffening of the quark matter EoS at high densities due to multiquark interactions and (3) possibility for a strong first order phase transition with an early onset and large density jump. The third feature results from a Maxwell construction for the possible transition from the nuclear to a quark matter phase and its properties depend on the two parameters used for (1) and (2), respectively. Varying these two parameters one obtains a class of hybrid EoS that yields solutions of the Tolman-Oppenheimer-Volkoff (TOV) equations for sequences of hadronic and hybrid stars in the mass-radius diagram which cover the full range of patterns according to the Alford-Han-Prakash classification following which a hybrid star branch can be either absent, connected or disconnected with the hadronic one. The latter case often includes a tiny connected branch. The disconnected hybrid star branch, also called "third family", corresponds to high-mass twin stars characterized by the same gravitational mass but different radii. We perform a Bayesian analysis and demonstrate that the observation of such a pair of high-mass twin stars would have a sufficient discriminating power to favor hybrid EoS with a strong first order phase transition over alternative EoS.

Aims: We present a new microscopic hadron-quark hybrid equation of state model for astrophysical applications, from which compact hybrid star configurations are constructed. These are composed of a quark core and a hadronic shell with a first-order phase transition at their interface. The resulting mass-radius relations are in accordance with the latest astrophysical constraints. Methods: The quark matter description is based on a quantum chromodynamics (QCD) motivated chiral approach with higher-order quark interactions in the Dirac scalar and vector coupling channels. For hadronic matter we select a relativistic mean-field equation of state with density-dependent couplings. Since the nucleons are treated in the quasi-particle framework, an excluded volume correction has been included for the nuclear equation of state at suprasaturation density which takes into account the finite size of the nucleons. Results: These novel aspects, excluded volume in the hadronic phase and the higher-order repulsive interactions in the quark phase, lead to a strong first-order phase transition with large latent heat, i.e. the energy-density jump at the phase transition, which fulfils a criterion for a disconnected third-family branch of compact stars in the mass-radius relationship. These twin stars appear at high masses ($\sim$ 2 M$_\odot$) that are relevant for current observations of high-mass pulsars. Conclusions: This analysis offers a unique possibility by radius observations of compact stars to probe the QCD phase diagram at zero temperature and large chemical potential and even to support the existence of a critical point in the QCD phase diagram.

We study the quark-hadron phase transition with the finite-size effects in neutron stars. The finite-size effects should be, generally, taken into account in the phase transition of multi-component system. The behavior of the phase transition, however, strongly depends on the models for quark and hadron matter, surface tension, neutrino fraction, and temperature. We find that, if the surface tension is strong, the EOS becomes similar to the case of a Maxwell construction for any hadron and/or quark model, though we adopt the Gibbs conditions. We also find that the mass-radius relations for that EOS are consistent with the observations, and our model is then applicable to realistic astrophysical phenomena such as the thermal evolution of compact stars.

We explore the possibility of obtaining heavy hybrid stars within the framework of the two flavor Nambu--Jona-Lasinio model with 8-quark interactions in the scalar and in the vector channel. The main impact of the 8-quark scalar channel is to reduce the onset of quark matter, while the 8-quark vector channel acts to stiff the equation of state at high densities. Within the parameter space where the 4-quark vector channel is small, and the 8-quark vector channel sizeable, stable stars with $2 M_\odot$ and above are found to hold quark matter in their cores.

We study the role of finite-size effects at the hadron-quark phase transition in a new hybrid equation of state constructed from an ab-initio Br\"uckner-Hartree-Fock equation of state with the realistic Bonn-B potential for the hadronic phase and a covariant non-local Nambu--Jona-Lasinio model for the quark phase. We construct static hybrid star sequences and find that our model can support stable hybrid stars with an onset of quark matter below $2 M_\odot$ and a maximum mass above $2.17 M_\odot$ in agreement with recent observations. If the finite-size effects are taken into account the core is composed of pure quark matter. Provided that the quark vector channel interaction is small, and the finite size effects are taken into account, quark matter appears at densities 2-3 times the nuclear saturation density. In that case the proton fraction in the hadronic phase remains below the value required by the onset of the direct URCA process, so that the early onset of quark matter shall affect on the rapid cooling of the star.

We study the possibility that the transition from hadron matter to quark matter at vanishing temperatures proceeds via crossover, similar to the crossover behavior found with lattice QCD studies at high temperatures. The purpose is to examine astrophysical consequences of this postulate by constructing hybrid star sequences fulfilling current experimental data.

We present two types of models for hybrid compact stars composed of a quark core and a hadronic mantle with an abrupt first order phase transition at the interface which are in accordance with the latest astrophysical measurements of two 2 M_sun pulsars. While the first is a schematic one, the second one is based on a QCD motivated nonlocal PNJL model with density-dependent vector coupling strength. Both models support the possibility of so called twin compact stars which have the same mass but different radius and internal structure at high mass (~2 M_sun), provided they exhibit a large jump \Delta \epsilon in the energy density of the first order phase transition fulfilling \Delta \epsilon/\epsilon_crit > 0.6. We conclude that the measurement of high-mass twin stars would support the existence of a first order phase transition in symmetric matter at zero temperature entailing the existence of a critical end point in the QCD phase diagram.

Nonlocal PNJL models allow for a detailed description of chiral quark dynamics with running quark masses and wave function renormalization in accordance with lattice QCD (LQCD) in vacuum. Their generalization to finite temperature T and chemical potential \mu{} allows to reproduce the \mu-dependence of the pseudocritical temperature from LQCD when a nonvanishing vector meson coupling is adjusted. This restricts the region for the critical endpoint in the QCD phase diagram and stiffens the quark matter equation of state (EoS). It is demonstrated that the construction of a hybrid EoS for compact star applications within a two-phase approach employing the nonlocal PNJL EoS and an advanced hadronic EoS leads to the masquerade problem. A density dependence of the vector meson coupling is suggested as a possible solution which can be adjusted in a suitable way to describe hybrid stars with a maximum mass in excess of 2 M_sun with a possible early onset of quark deconfinement even in the cores of typical (M ~ 1.4 M_sun) neutron stars.