Gravitational wave observations of GW170817 placed bounds on the tidal deformabilities of compact stars, allowing one to probe equations of state for matter at supranuclear densities. Here we design new parametrizations for hybrid hadron-quark equations of state, which give rise to low-mass twin stars, and test them against GW170817. We find that GW170817 is consistent with the coalescence of a binary hybrid star-neutron star. We also test and find that the I-Love-Q relations for hybrid stars in the third family agree with those for purely hadronic and quark stars within $\sim 3\%$ for both slowly and rapidly rotating configurations, implying that these relations can be used to perform equation-of-state independent tests of general relativity and to break degeneracies in gravitational waveforms for hybrid stars in the third family as well.

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 (~2 M ⊙) 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.

A class of hybrid compact star equations of state is investigated that joins by a Maxwell construction a low-density phase of hadronic matter, modeled by a relativistic mean-field approach with excluded nucleon volume, with a high-density phase of color superconducting two-flavor quark matter, described within a nonlocal covariant chiral quark model. It is found that the occurrence of a stable branch of hybrid compact stars requires a nonvanishing vector meson coupling in the quark model that exceeds a minimal value which depends on the presence of a diquark condensate. It is shown that these hybrid stars do not form a third family disconnected from the second family of ordinary neutron stars unless additional (de)confining effects are introduced with a density-dependent bag pressure. A suitably chosen density dependence of the vector meson coupling assures that at the same time the $2M_{\odot }$ maximum mass constraint is fulfilled on the hybrid star branch. A twofold interpolation method is realized which implements both the density dependence of a confining bag pressure at the onset of the hadron-to-quark matter transition and the stiffening of quark matter at higher densities by a density-dependent vector meson coupling. For three parametrizations of this class of hybrid equation of state the properties of corresponding compact star sequences are presented, including mass twins of neutron and hybrid stars at 2.00, 1.39 and $1.20M_{\odot }$, respectively, and the hybrid compact star (third) families. The sensitivity of the hybrid equation of state and the corresponding compact star sequences to variations of the interpolation parameters at the 10% level is investigated and it is found that the feature of third family solutions for compact stars is robust against such a variation. This advanced description of hybrid star matter allows us to interpret GW170817 as a merger not only of two neutron stars but also of a neutron star with a hybrid star or of two hybrid stars.

We show that in the three-polytrope model of Hebeler et al. [ Astrophys. J. 773, 11 (2013)] for the neutron star equation of state (EoS) at supersaturation densities a third family of compact stars can be obtained which confirms the possibility of high-mass twin stars that have coincident masses $M_1=M_2\approx 2M_{\odot }$ and significantly different radii $|R_1-R_2|>2–3$ km. We show that the causality constraint puts severe limitations on the maximum mass of the third family sequence which can be relaxed when this scheme is extended to four polytropes thus mimicking a realistic high-density matter EoS.

We briefly review the five possible real polynomial solutions of hypergeometric differential equations. Three of them are the well known classical orthogonal polynomials, but the other two are different with respect to their orthogonality properties. We then focus on the family of polynomials which exhibits a finite orthogonality. This family, to be referred to as the Romanovski polynomials, is required in exact solutions of several physics problems ranging from quantum mechanics and quark physics to random matrix theory. It appears timely to draw attention to it by the present study. Our survey also includes several new observations on the orthogonality properties of the Romanovski polynomials and new developments from their Rodrigues formula.