With the first detection of gravitational waves from a binary system of neutron stars, GW170817, a new window was opened to study the properties of matter at and above nuclear-saturation density. Reaching densities a few times that of nuclear matter and temperatures up to $100\,\rm{MeV}$, such mergers also represent potential sites for a phase transition (PT) from confined hadronic matter to deconfined quark matter. While the lack of a postmerger signal in GW170817 has prevented us from assessing experimentally this scenario, two theoretical studies have explored the postmerger gravitational-wave signatures of PTs in mergers of binary systems of neutron stars. We here extend and complete the picture by presenting a novel signature of the occurrence of a PT. More specifically, using fully general-relativistic hydrodynamic simulations and employing a suitably constructed equation of state that includes a PT, we present the occurrence of a "delayed PT", i.e. a PT that develops only some time after the merger and produces a metastable object with a quark-matter core, i.e. a hypermassive hybrid star. Because in this scenario, the postmerger signal exhibits two distinct fundamental gravitational-wave frequencies -- before and after the PT -- the associated signature promises to be the strongest and cleanest among those considered so far, and one of the best signatures of the production of quark matter in the present Universe.

We study in detail the nuclear aspects of a neutron-star merger in which deconfinement to quark matter takes place. For this purpose, we make use of the Chiral Mean Field (CMF) model, an effective relativistic model that includes self-consistent chiral symmetry restoration and deconfinement to quark matter and, for this reason, predicts the existence of different degrees of freedom depending on the local density/chemical potential and temperature. We then use the out-of-chemical-equilibrium finite-temperature CMF equation of state in full general-relativistic simulations to analyze which regions of different QCD phase diagrams are probed and which conditions, such as strangeness and entropy, are generated when a strong first-order phase transition appears. We also investigate the amount of electrons present in different stages of the merger and discuss how far from chemical equilibrium they can be and, finally, draw some comparisons with matter created in supernova explosions and heavy-ion collisions.

If a phase transition is allowed to take place in the core of a compact star, a new stable branch of equilibrium configurations can appear, providing solutions with the same mass as the purely hadronic branch and hence giving rise to twin-star configurations. We perform an extensive analysis of the features of the phase transition leading twin-star configurations and, at the same time, fulfilling the constraints coming from the maximum mass of $2M_\odot$ and the information following gravitational-wave event GW170817. In particular, we use a general equation of state for the neutron-star matter that parametrizes the hadron-quark phase transition between the model describing the hadronic phase and a constant speed of sound for the quark phase. We find that the largest number of twin-star solutions has masses in the neutron-star branch in the range $1-2M_\odot$ and twin-branch masses $\gtrsim 2M_\odot$. The analysis of the masses, radii and tidal deformabilities also reveals that when twin stars appear, the tidal deformability shows two distinct branches with the same mass, thus differing considerably from the behaviour expected for neutron stars. In addition, we find that the data from GW170817 is compatible with the existence of hybrid stars and could also be interpreted as produced by the merger of a binary system of hybrid stars or of a hybrid star with a neutron star. The presence of a hybrid star in the inspiral phase can be established clearly if future gravitational-wave detections measure chirp masses $\mathcal{M}\lesssim 1.2M_\odot$ and tidal deformabilities of $\Lambda_{1.4}\lesssim 400$ for $1.4M_\odot$ stars. Finally, combining all observational information available, we set constraints on the parameters that characterise the phase transition, the maximum masses, and the radii of $1.4M_\odot$ stars described by equations of state leading to twin-star configurations.

In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals.

Merging binaries of neutron stars are not only strong sources of gravitational waves, but also have the potential of revealing states of matter at densities and temperatures not accessible in laboratories. A crucial and long-standing question in this context is whether quarks are deconfined as a result of the dramatic increase in density and temperature following the merger. We present the first fully general-relativistic simulations of merging neutron stars including quarks at finite temperatures that can be switched off consistently in the equation of state. Within our approach, we can determine clearly what signatures a quark-hadron phase transition would leave in the gravitational-wave signal. In particular, we show that if the conditions are met for a phase transition to take place at several times nuclear saturation density, they would lead to a post-merger signal considerably different from the one expected from the inspiral, that can only probe the hadronic part of the equations of state, and to an anticipated collapse of the merged object. We also show that the phase transition leads to a very hot and dense quark core that, when it collapses to a black hole, produces a ringdown signal different from the hadronic one. Finally, in analogy with what is done in heavy-ion collisions, we use the evolution of the temperature and density in the merger remnant to illustrate the properties of the phase transition in a QCD phase diagram.

Inferring the properties of dense matter is one of the most exciting prospects from the measurement of gravitational waves from neutron star mergers. However, it will require reliable numerical simulations that incorporate viscous dissipation and energy transport if these can play a significant role within the survival time of the post-merger object. We calculate timescales for typical forms of dissipation and find that thermal transport and shear viscosity will not be important unless neutrino trapping occurs, which requires temperatures above about 10 MeV and gradients over lengthscales of 0.1 km or less. On the other hand, if direct-Urca processes remain suppressed, leaving modified-Urca processes to establish flavor equilibrium, then bulk viscous dissipation could provide significant damping to density oscillations observed right after the merger. When comparing with data from a state-of-the-art merger simulation we find that the bulk viscosity takes values close to its resonant maximum in a typical neutron-star merger, motivating a more careful assessment of the role of bulk viscous dissipation in the gravitational-wave signal from merging neutron stars.

Determining the differential-rotation law of compact stellar objects produced in binary neutron stars mergers or core-collapse supernovae is an old problem in relativistic astrophysics. Addressing this problem is important because it impacts directly on the maximum mass these objects can attain and hence on the threshold to black-hole formation under realistic conditions. Using the results from a large number of numerical simulations in full general relativity of binary neutron star mergers described with various equations of state and masses, we study the rotational properties of the resulting hypermassive neutron stars. We find that the angular-velocity distribution shows only a modest dependence on the equation of state, thus exhibiting the traits of "quasi-universality" found in other aspects of compact stars, both isolated and in binary systems. The distributions are characterized by an almost uniformly rotating core and a "disk". Such a configuration is significantly different from the $j-{\rm constant}$ differential-rotation law that is commonly adopted in equilibrium models of differentially rotating stars. Furthermore, the rest-mass contained in such a disk can be quite large, ranging from $\simeq 0.03\,M_{\odot}$ in the case of high-mass binaries with stiff equations of state, up to $\simeq 0.2\,M_{\odot}$ for low-mass binaries with soft equations of state. We comment on the astrophysical implications of our findings and on the long-term evolutionary scenarios that can be conjectured on the basis of our simulations.

The inner regions of the most massive compact stellar objects might be occupied by a phase of quarks. Since the observations of the massive pulsars PSR J1614-2230 and of PSR J0348+0432 with about two solar masses, the equations of state constructing relativistic stellar models have to be constrained respecting these new limits. We discuss stable hybrid stars, i.e. compact objects with an outer layer composed of nuclear matter and with a core consisting of quark matter (QM). For the outer nuclear layer we utilize a density dependent nuclear equation of state and we use a chiral SU(3) Quark-Meson model with a vacuum energy pressure to describe the objects core. The appearance of a disconnected mass-radius branch emerging from the hybrid star branch implies the existence of a third family of compact stars, so called twin stars. Twin stars did not emerge as the transition pressure has to be relatively small with a large jump in energy density, which could not be satisfied within our approach. This is, among other reasons, due to the fact that the speed of sound in QM has to be relatively high, which can be accomplished by an increase of the repulsive coupling. This increase on the other hand yields too high transition pressures for twins stars to appear.

We have constructed models of uniformly rotating compact stars including hyperons, Bose-Einstein condensates of antikaons and quarks. First order phase transitions from hadronic to antikaon condensed matter and then to quark matter are considered here. For the equation of state undergoing phase transitions to antikaon condensates, the third family of compact stars are found to exist in the fixed angular velocity sequences. However, the third family solution disappears when the compact stars rotate very fast. For this equation of state, the fixed baryon number supramassive sequence shows a second stable part after the unstable region but no back bending phenomenon. On the other hand, we observe that the rotation gives rise to a second maximum beyond the neutron star maximum for the equation of state involving phase transitions to both antikaon condensed and quark matter. In this case, the back bending phenomenon has been observed in the supramassive sequence as a consequence of the first order phase transition from $K^-$ condensed to quark matter. And the back bending segment contains stable configurations of neutron stars.

The quark hadron phase transition in a rotating compact star has been studied. The NLZ-model has been used for the hadronic sector and the MIT Bag model has been used for the quark sector. It has been found that rotating twin star (third family) solutions are obtained upto $\Omega \approx 4000 s^{-1}$. Stars which are rotating faster than this limit do not show any twin star solution. A backbending in moment of inertia is also observed in the supermassive rest mass sequences. The braking index is found to diverge for a star having pure quark core.

Baryonic matter at high density and low temperature is a color superconductor. In real world, this state of matter may naturally appear inside compact stars. A construction of a hybrid compact star with two flavor color superconducting quark matter in its interior is presented.

Realistic nonstrange hybrid compact stars with color superconducting quark matter in their interior are constructed. It is shown that a positively charged two-flavor color superconducting phase could naturally appear in the core of a hybrid star as one of the components of a globally neutral mixed phase. The negatively charged normal quark phase is the other component of the mixed phase. The quark core of the star is surrounded by another mixed phase made of hadronic and normal quark matter. The two mixed phases are separated by a sharp interface. Finally, the lowest density regions of the star are made of pure hadronic matter and nuclear crust.

We study properties of compact stars with the deconfinement phase transition in their interiors. The equation of state of cold baryon-rich matter is constructed by combining a relativistic mean-field model for the hadronic phase and the MIT Bag model for the deconfined phase. In a narrow parameter range two sequences of compact stars (twin stars), which differ by the size of the quark core, have been found. We demonstrate the possibility of a rapid transition between the twin stars with the energy release of about $10^{52}$ ergs. This transition should be accompanied by the prompt neutrino burst and the delayed gamma-ray burst.

Recent progress in the understanding of the high density phase of neutron stars advances the view that a substantial fraction of the matter consists of hyperons. The possible impacts of a highly attractive interaction between hyperons on the properties of compact stars are investigated. We find that a hadronic equation of state with hyperons allows for a first order phase transition to hyperonic matter. The corresponding hyperon stars can have rather small radii of R=8 km.

We investigate the properties of charge-neutral $\beta$-equilibrium cold quark matter within the Nambu-Jona-Lasinio model. The calculations are carried out for different ratios of coupling constants characterizing the vector and scalar 4-fermion interaction, $\xi\equiv G_V/G_S$. It is shown that for $\xi < 0.4$ matter is self-bound and for $\xi < 0.65$ it has a first order phase transition of the liquid-gas type. The Gibbs conditions in the mixed phase are applied for the case of two chemical potentials associated with the baryon number and electric charge. The characteristics of the quark stars are calculated for $\xi = 0, 0.5$ and 1. It is found that the phase transition leads to a strong density variation at the surface of these stars. For $\xi = 1$ the properties of quark stars show behaviors typical for neutron stars. At $\xi\goo 0.4$ the stars near to the maximum mass have a large admixture of strange quarks in their interiors.

We investigate various properties of neutron star matter within an effective chiral $SU(3)_L \times SU(3)_R$ model. The predictions of this model are compared with a Walecka-type model. It is demonstrated that the importance of hyperon degrees are strongly depending on the interaction used, even if the equation of state near saturation density is nearly the same in both models. While the Walecka-type model predicts a strange star core with strangeness fraction $f_S \approx 4/3$, the chiral model allows only for $f_S \approx 1/3$ and predicts that $\Sigma^0$, $\Sigma^+$ and $\Xi^0$ will not exist in star, in contrast to the Walecka-type model.

The effects of internal quark structure of baryons on the composition and structure of neutron star matter with hyperons are investigated in the quark-meson coupling (QMC) model. The QMC model is based on mean-field description of nonoverlapping spherical bags bound by self-consistent exchange of scalar and vector mesons. The predictions of this model are compared with quantum hadrodynamic (QHD) model calibrated to reproduce identical nuclear matter saturation properties. By employing a density dependent bag constant through direct coupling to the scalar field, the QMC model is found to exhibit identical properties as QHD near saturation density. Furthermore, this modified QMC model provides well-behaved and continuous solutions at high densities relevant to the core of neutron stars. Two additional strange mesons are introduced which couple only to the strange quark in the QMC model and to the hyperons in the QHD model. The constitution and structure of stars with hyperons in the QMC and QHD models reveal interesting differences. This suggests the importance of quark structure effects in the baryons at high densities.