We investigate the possibility that GW170817 has not been the merger of two conventional neutron stars (NS) but involved at least one if not two hybrid stars with a quark matter core which 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 NICER. 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 of deconfinement and can mimic the effects of pasta phases with the possibility of a third family of hybrid stars in the mass-radius ($M-R$) diagram. The main findings of this study are that: 1) the presence of multiple configurations for a given mass (twins) corresponds to a set of disconnected lines in the 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_{rm onset}le 0.8~M_odot$ with a $M-R$ relationship that in the range of observed neutron star masses is almost indistinguishable from that of a soft hadronic APR EoS; 3) a few yet fictitious measurements of the NICER experiment with a $1sigma$ range that is half of the present value and different mass and radius would change the posterior likelihood so that hybrid EoS with a phase transition onset in the range $M_{rm onset} = 1.1 - 1.6~M_odot $ would be favored.
We present a first step in developing a benchmark equation-of-state (EoS) model for multi-messenger astronomy that unifies the thermodynamics of quark and hadronic degrees of freedom. A Lagrangian approach to the thermodynamic potential of quark-meso
n-nucleon (QMN) matter was used. In this approach, dynamical chiral-symmetry breaking is described by the scalar mean-field dynamics coupled to quarks and nucleons and their chiral partners, whereby its restoration occurs in the hadronic phase by parity doubling, as well as in the quark phase. Quark confinement was achieved by an auxiliary scalar field that parametrizes a dynamical infrared cutoff in the quark sector, serving as an ultraviolet cutoff for the nucleonic phase space. The gap equations were solved for the isospin-symmetric case, as well as for neutron star (NS) conditions. We also calculated the mass-radius (MR) relation of NSs and their tidal deformability parameter. The obtained EoS is in accordance with nuclear matter properties at saturation density and with the flow constraint from heavy ion collision experiments. For isospin-asymmetric matter, a sequential occurrence of light quark flavors is obtained, allowing for a mixed phase of chirally-symmetric nucleonic matter with deconfined down quarks. The MR relations and TDs for compact stars fulfill the constraints from the latest astrophysical observations for PSR J0740+6620, PSR J0030+0451, and the NS merger GW170817, whereby the tension between the maximum mass and compactness constraints rather uniquely fixes the model parameters. The model predicts the existence of stars with a core of chirally restored but purely hadronic (confined) matter for masses beyond $1.8~M_odot$. Stars with pure-quark matter cores are found to be unstable against the gravitational collapse. This instability is shifted to even higher densities if repulsive interactions between quarks are included.
In this chapter we will introduce an effective equation of state (EoS) model based on polytropes that serves to study the so called mass twins scenario, where two compact stars have approximately the same mass but (significant for observation) quite
different radii. Stellar mass twin configurations are obtained if a strong first-order phase transition occurs in the interior of a compact star. In the mass-radius diagram of compact stars, this will lead to a third branch of gravitationally stable stars with features that are very distinctive from those of white dwarfs and neutron stars. We discuss rotating hybrid star sequences in the slow rotation approximation and in full general relativity and draw conclusions for an upper limit on the maximum mass of nonrotating compact stars that has recently be deduced from the observation of the merger event GW170817.
We extend the recently developed hybrid quark-meson-nucleon model by augmenting a six-point scalar interaction and investigate the consequences for neutron-star sequences in the mass-radius diagram. The model has the characteristic feature that, at i
ncreasing baryon density, the chiral symmetry is restored within the hadronic phase by lifting the mass splitting between chiral partner states (parity doubling), before quark deconfinement takes place. At low temperature and finite baryon density, the model predicts a first-, second-order chiral phase transition, or a crossover, depending on the expectation value of the scalar field, and a first-order deconfinement phase transition. We discuss two sets of free parameters, which result in compact-star mass-radius relations that are at tension with the combined constraints for maximum-mass ($2~M_odot$) and the compactness (GW170817). We find that the most preferable mass-radius relations result in isospin-symmetric phase diagram with rather low temperature for the critical point of the chiral phase transition.
We describe a multi-messenger interpretation of GW170817, which yields a robust lower limit on NS radii. This excludes NSs with radii smaller than about 10.7 km and thus rules out very soft nuclear matter. We stress the potential of this type of cons
traints when future detections become available. A very similar argumentation may yield an upper bound on the maximum mass of nonrotating NSs. We also discuss simulations of NS mergers, which undergo a first-order phase transition to quark matter. We point out a different dynamical behavior. Considering the gravitational-wave signal, we identify an unambiguous signature of the QCD phase transition in NS mergers. The occurrence of quark matter through a strong first-order phase transition during merging leads to a characteristic shift of the dominant postmerger frequency. The frequency shift is indicative for a phase transition if it is compared to the postmerger frequency which is expected for purely hadronic EoS models. A very strong deviation of several 100 Hz is observed for hybrid EoSs in an otherwise tight relation between the tidal deformability and the postmerger frequency. We address the potential impact of a first-order phase transition on the electromagnetic counterpart of NS mergers. Our simulations suggest that there would be no significant qualitative differences between a system undergoing a phase transition to quark matter and purely hadronic mergers. The quantitative differences are within the spread which is found between different hadronic EoS models. This implies on the one hand that GW170817 is compatible with a possible transition to quark matter. On the other hand these considerations show that it may not be easy to identify quantitative differences between purely hadronic mergers and events in which quark matter occurs considering solely their electromagnetic counterpart or their nucleosynthesis products. (abridged)
Recent lattice QCD studies at vanishing density exhibit the parity-doubling structure for the low-lying baryons around the chiral crossover temperature. This finding is likely an imprint of the chiral symmetry restoration in the baryonic sector of QC
D, and is expected to occur also in cold dense matter, which makes it of major relevance for compact stars. By contrast, typical effective models for compact star matter embody chiral physics solely in the deconfined sector, with quarks as degrees of freedom. In this contribution, we present a description of QCD matter based on the effective hybrid quark-meson-nucleon model. Its characteristic feature is that, under neutron-star conditions, the chiral symmetry is restored in a first-order phase transition deep in the hadronic phase, before the deconfinement of quarks takes place. We discuss the implications of the parity doubling of baryons on the mass-radius relation for compact stars obtained in accordance with the modern constraints on the mass from PSR J0348+0432, the compactness from GW170817, as well as the direct URCA process threshold. We show that the existence of high-mass stars might not necessarily signal the deconfinement of quarks.
We investigate the equation of state for a recently developed hybrid quark-meson-nucleon model under neutron star conditions of $beta-$equilibrium and charge neutrality. The model has the characteristic feature that at increasing baryon density chira
l symmetry is restored in a first order transition within the hadronic phase by lifting the mass splitting between chiral partner states, before quark deconfinement takes place. Most important for this study are the nucleon (neutron, proton) and $N(1535)$ states. We present three sets for the two free parameters which result in compact star mass-radius relations in accordance with modern constraints on the mass from PSR~J0437-4715 and on the compactness from GW170817. We also consider the threshold for the direct URCA process for which a new relationship is given and suggest as an additional constraint on the parameter choice of the model that this process shall become operative at best for stars with masses above the range for binary radio pulsars, $M>1.4~M_odot$.
In this review article we discuss selected developments regarding the role of the equation of state (EOS) in simulations of core-collapse supernovae. There are no first-principle calculations of the state of matter under supernova conditions since a
wide range of conditions is covered, in terms of density, temperature and isospin asymmetry. Instead, model EOS are commonly employed in supernova studies. These can be divided into regimes with intrinsically different degrees of freedom: heavy nuclei at low temperatures, inhomogeneous nuclear matter where light and heavy nuclei coexist together with unbound nucleons, and the transition to homogeneous matter at high densities and temperatures. In this article we discuss each of these phases with particular view on their role in supernova simulations.
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 p
hase 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 discuss an idea for how accreting millisecond pulsars could contribute to the understanding of the QCD phase transition in the high-density nuclear matter equation of state (EoS). It is based on two ingredients, the first one being a ``phase diagr
am of rapidly rotating compact star configurations in the plane of spin frequency and mass, determined with state-of-the-art hybrid equations of state, allowing for a transition to color superconducting quark matter. The second is the study of spin-up and accretion evolution in this phase diagram. We show that the quark matter phase transition leads to a characteristic line in the Omega-M plane, the phase border between neutron stars and hybrid stars with a quark matter core. Along this line a change in the pulsars moment of inertia entails a waiting point phenomenon in the accreting millisecond X-ray pulsar (AMXP) evolution: most of these objects should therefore be found along the phase border in the Omega-M plane, which may be viewed as the AMXP analog of the main sequence in the Hertzsprung-Russell diagram for normal stars. In order to prove the existence of a high-density phase transition in the cores of compact stars we need population statistics for AMXPs with sufficiently accurate determination of their masses and spin frequencies.
