Abstract. 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 highmass 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. PACS

We investigate the possibility that GW170817 was not the merger of two conventional neutron stars (NS), but involved at least one if not two hybrid stars with a quark matter core that might even belong to a third family of compact stars. To this end, we develop a Bayesian analysis method for selecting the most probable equation of state (EoS) under a set of constraints from compact star physics, which now also include the tidal deformability from GW170817 and the first result for the mass and radius determination for PSR J0030+0451 by the NICER Collaboration. We apply this method for the first time to a two-parameter family of hybrid EoS based on the DD2 model with nucleonic excluded volume for hadronic matter and the color superconducting generalized nlNJL model for quark matter. The model has a variable onset density for deconfinement and can mimic the effects of pasta phases with the possibility of producing a third family of hybrid stars in the mass-radius diagram. The main findings of this study are that: (1) the presence of multiple configurations for a given mass (twins or even triples) corresponds to a set of disconnected lines in the Λ 1 – Λ 2 diagram of tidal deformabilities for binary mergers, so that merger events from the same mass range may result in a probability landscape with different peak positions; (2) the Bayesian analysis with the above observational constraints favors an early onset of the deconfinement transition, at masses of M onset ≤ 0.8 M ⊙ with an M–R relationship that in the range of observed neutron star masses is almost indistinguishable from that of a soft hadronic Akmal, Pandharipande, and Ravenhall (APR) EoS; (3) a few, yet fictitious measurements of the NICER experiment two times more accurate than the present value and a different mass and radius that would change the posterior likelihood so that hybrid EoS with a phase transition onset in the range M onset = 1.1–1.6 M ⊙ would be favored.

We investigate the robustness of third family solutions for hybrid compact stars with a quark matter core that correspond to the occurrence of high-mass twin stars against a softening of the phase transition by means of a construction that mimics the effects of pasta structures in the mixed phase. We consider a class of hybrid equations of state that exploits a relativistic mean-field model for the hadronic as well as for the quark matter phase. We present parametrizations that correspond to branches of high-mass twin star pairs with maximum masses between 2.05 M and 1.48 M having radius differences between 3.2 and 1.5 km, respectively. When compared to a Maxwell construction with a fixed value of critical pressure Pc, the effect of the mixed phase construction consists in the occurrence of a region of pressures around Pc belonging to the coexistence of hadronic and quark matter phases between the onset pressure at PH and the end of the transition at PQ. The maximum broadening which would still allow mass twin compact stars is found to be (PQ − PH )max ≈ Pc for all parametrizations within the present class of models. At least the heavier of the neutron stars of the binary merger GW170817 could have been a member of the third family of hybrid stars. We present the example of another class of hybrid star equations of state for which the appearance of the third family branch is not as robust against mixed phase effects as that of the present work.PACS numbers:

The effect of pasta phases on the quark-hadron phase transition is investigated for a set of relativistic mean-field equations of state for both hadron and quark matter. The results of the full numerical solution with pasta phases are compared with those of an interpolating construction used in previous works, for which we demonstrate an adequate description of the numerical results. A one-to-one mapping of the free parameter of the construction to the physical surface tension of the quark-hadron interface is obtained for which a fit formula is given. For each pair of quark and hadron matter models the critical value of the surface tension is determined, above which the phase transition becomes close to the Maxwell construction. This result agrees well with earlier theoretical estimates. The study is extended to neutron star matter in beta equilibrium with electrons and muons and is applied to investigate the effect of pasta phases on the structure of hybrid compact stars and the robustness of a possible third family solution.